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and  C.15  of  Appendix  C. 


Abbreviations 

ASCE 

ATFC 
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CK  U 
o 

DSS 
IBM 
ICES-LEASE 
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JHRP 
MP 

NYSTAB 

OCR 

PI 

PSA 

PSP 

SHANSEP 


American  Society  of  Civil  Engineers . 

American  Testing  and  Engineering  Corporation. 

Control  Data  Corporation  (Computer  System  at  Purdue  Univ. ) 

K  Consolidated  -  Undrained  Shear  Test. 
o 

Direct  Simple  Shear. 

International  Business  Machines  Corporation  (Computer). 

Integrated  Civil  Engineering  System  -  Lease  Electronic 

Accounting  System. 

Indiana  State  Highway  Commission. 

Joint  Highway  Research  Project. 

Morgenstern-Price  method  of  slices  for  slope  stability 

analysis. 

New  York  State,  Department  of  Transportation,  Slope 

Stability  Analysis. 

Over  Consolidation  Ratio. 

Plasticity  Index. 

Plane  Strain  Active. 

Plane  Strain  Passive. 

Stress  History  and  Normalized  Soil  Engineering  Properties. 
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Special  Symbols 

A  prime   ('"")   indicates   that  the  variable  is   in  terms  of  effective  stresses 

A     (Delta)   indicates   a  finite  change  or  variation  in  a  variable. 

L     (Sigma)   used  as  summation  symbol. 

£    Symbol  for  Center  Line. 

Symbols 

a         -  Moment  arm  for  body  forces  in  simple  slope  stability 

analysis. 
a.  to  a    -  Geometry  dependent  coefficients  in  standardized  equation 

for  the  factor  of  safety. 
A1  ,A  ,A    -  Coefficients  introduced  in  the  general  equation  for  the 

factor  of  safety. 
A,B,C,D.E  -  Zones  of  constant  stress. 

-  Different  types  of  analysis  (see  Table  B.M. 
b         -  Base  of  a  slope. 
B         -  Crown  width  of  an  embankment. 
c,  c'      -  Strength  intercept  in  terms  of  total  (c)  and  effective  (c') 

stresses,  respectively. 

c,  to  e    -  Strength  intercepts  of  soils  1  to  n. 
1     n 


I      w  \ 

c  -  Normalized  strength  of  weak  seam  (=— )  • 
n  c£ 

c  -  Required  strength  intercept  to  make  slope  in  equilibrium, 
r 

c  -  Undrained  strength  intercept. 

C,  C  -  Stress  independent  strength  (=c  dl). 

C  =  c  dl 
r  r 

A(-er)  -  The   change  in  —  due  to  a  change     Ac'    in  c'      (=  ^jj"~)' 
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A(-^-)Y     -  The  change  in  —  due  to  a  change  Ly   in  y 

(-    c'         C   V 
(Y  +^y)  h     yh    ' 

Yh^h     "  The  chanSe  in  ~—  due  to  a  change  Ah  in  h 

(-     c'  C'v 

v"  Y  <h  +  Ah)  "  Yh 

d         -  Maximum  distance  between  a  shear  surface  and  the  chord 

of  the  shear  surface,  measured  normal  to  the  chord. 
D         -  Depth  to  a  weak  layer 
E         -  Interslice  normal  forces 

-  As  subscript,  indicates  embankment. 

f         -  Used  as  subscript,  indicates  at  failure  or  on  the  failure 

plane. 

ff  -     Used  as  subscript,   indicates  on   failure  plane  at   failure. 

f  -     Partial  coefficient  for  skin  friction  on  piles . 

f  -     Partial  coefficient  with  respect  to  strength  intercept   c. 

f  -     Partial  coefficient    for  dead  load  and  water  pressure, 

g 

f  -     Partial   coefficient  for  live  loads. 

P 

f  -     Partial  coefficient   for  wind  loads. 

f,  -     Partial  coefficient  with  respect  to  strength  angle  $. 

f(x)  -     A   function  relating  the  horizontal  and  vertical  interslice 

forces   along  a  shear  surface  as  proposed  by  Morgenstern 

and  Price   (1965) . 
F  -     Factor  of  safety. 

-  Used  as   subscript,   indicates   foundation. 

F  -     Modified  Bishop  factor  of  safety  for  circular  failure 

surface . 

-  Janbu's   factor  of  safety   (Figure  6.20). 
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F^       -  The  value  of  F  ,  when  moment  is  taken  at  infinity 

(=  l^m  F  )  =  Janbu's  factor  of  safety. 

F        -  Factor  of  safety  with  respect  to  strength  intercept  c, 

F        -  Normalized  factor  of  safety  (=  -^-  ). 

CE 
-  Factor  of  safety  with  respect  to  strength  angle  <j). 


n 


AF(Ax  )   -  The  change  in  the  factor  of  safety  due  to  a  change  in  the 
independent  variable  x, ,  other  parameters  kept  constant . 


AF  AF(Axi) 


The  tangents  or  secants  of  the  F  vs.  x.  curves, 


Ax.  Ax.  i 

1  1 

other  variables  held  constant. 

AF   ,  =  AF(Ac') 

c' 

AFU  =  AF(Ah) 

n 

AFg  =  AF(A6) 

AF^  =  AF(Ay) 

AF,,  =  AF(A4>') 

EAF  =  EAF(Ax.),  i  «  1,  n. 

1  J 

h,  H  -  Height  of  a  slope  or  embankment. 

H  -  Critical  height  of  a  slope. 

i  -  Subscript  indicating  an  array  of  variables. 

-  Inclination  of  ground  surface. 

k,    to  k       -  Geometry  dependent  coefficients   in  standardized  equation 
In 

for  the   factor  of  safety. 

k,  -  Horizontal  earthquake  coefficient. 

k  -  Vertical  earthquake  coefficient. 

K  -  Coefficient  of  earth  pressure  at  rest. 

o 

1  -  Length  of  a  failure  surface. 

dl  -  Infinitesimal  increment  of  length  along  the   failure  surface. 


M 
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L  -     Length  of  the   chord  of  a  shear  surface. 

-  Moment   around  the  center  0  of  a  circular  or  log.   spiral 
failure  surface. 

n  -     Number  of  slices  in  a  method  of  slices  slope  stability 

analysis . 

-  Depth  factor  (Taylor's  stability  chart). 

N»  N'  -     Total  and  effective  normal  forces  acting  on  a  shear  surface, 

0  -     Center  of  a  circular  or  log.   spiral  failure  surface. 

q.»  -     Maximum  shear  stress   at  failure   (=  ;r  (a,   -  a„)    ). 

f  2       1         3  f 

Q  -  Resultant  of  uniform  surcharge,   acting  on  top  of  a  alice. 

r  -  Radius  of  a  circular  or  log.   spiral  failure  surface. 

rQ  -  Initial  radius  of  log.   spiral  failure  surface. 

rM  -  Maximum  radius   of  log.   spiral  failure  surface, 

r  -  Pore  pressure  parameter  (=   ^—) . 

R  -  Resultant  of  normal  stress    (force)   and  normal  stress   depen- 

dent limiting  shear  stress    (force)   on  a  shear  surface. 

-  Ranf(0.). 

Ranf(x)       -     Standard  random  function,   taking  any  value  between  0 
and  l(for  x  =  0. ). 

R  -     Resultant  of  strength  intercept  resistance   forces. 

te 
R.  -     i       number  in  a  random  sequence  of  Ranf(0.). 

s    __  -     Shear  strength  based  on  effective  stresses. 

s  -     Undrained  shear  strength,  based  on  total  stresses. 

S,  S  -     Shear  resistance  required  along  a  potential  shear  surface 

for  equilibrium  of  slope. 

S  -     Shear  resistance  available  along  an  assumed  shear  surface. 


XXX 

u  -  Pore  water  pressure. 

u  -  Pore  pressure  constant, 
c 

U  -  Degree  of  consolidation. 

U  -  Water  force  acting  on  base  of  a  slice. 

Ug  -  Water  force  acting  on  top  of  a  slice. 

W  -  Used  as  subscript,  indicates  weak  soil  layer. 

W,  W'  -  Total  and  effective  weight  of  a  sliding  mass,  or  weight  of 

an  individual  slice. 

x.  ,  x  -  Independent  parameters  for  stability  analysis. 

Ax  -  Width  of  a  slice. 

X  -  Interslice  shear  force. 

y  .  -  Minimum  elevation  to  which  a  failure  surface  may  extend, 
nun 


Greek  Alphabet 

a        -  Inclination  of  shear  surface  with  respect  to  the  horizontal. 

-  Inclination  of  weak  soil  layer. 

ct-       -  Inclination  of  shear  surface  at  intersection  with  top  of 

s lope . 
6        -  Slope  inclination. 

Bj.,  SR    -  Inclination  of  left  and  right  side  of  an  embankment. 
Y        -  Unit  weight  of  a  soil. 
6        -  Inclination  of  uniform  surcharge,  acting  on  top  of  a  slice, 

measured  positive  counterclockwise  from  vertical. 
9        -  One  half  the  center  angle  of  a  circular  failure  surface. 

-  Angle  between  r  and  r  for  a  log.  spiral  surface. 

-  Initial  inclination  of  a  failure  surface  at  bottom  of  slope. 

-  Inclination  of  a  line  segment  of  a  trial  failure  surface. 
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Q-r  -     Lower  angle  limitation  for  line  segment  of  irregular 

failure  surface  (=  clockwise  direction  limit). 
9U        -  Upper  angle  limitation  for  line  segment  of  irregular 

failure  surface  (=  counterclockwise  direction  limit). 
9q        -  Central  angle  of  circular  failure  surface,  measured  from 

horizontal  radius  above  top  of  slope  to  radius  through 

intersection  of  failure  surface  with  ground  surface  at 

the  top  of  the  slope.   See  Figure  C.lU. 
A9        -  Angle  deflection  between  consecutive  line  segments  of  a 

circular  failure  surface. 
A0        -  Maximum  angle  deflection  between  consecutive  line  segments 

of  a  circular  failure  surface. 

AG  .      -  Minimum  angle  deflection  between  consecutive  line  segments 
mm  &  ■' 

of  a  circular  failure  surface. 

X        -  Stability  parameter  used  by  Janbu  (195^+b),  (=   t~  ) . 

a  ,  o  -  Maximum  and  minimum  principal  stresses . 

0,0'        -     Total  and  effective  normal  stress  on  failure  plane. 

o,        -  Horizontal  normal  stress, 
h 

o  -     Vertical  normal  stress, 

v 

a  '       -   Initial  vertical  effective  normal  stress, 
vo 

0    '       -  Vertical  consolidation  stress, 
vc 

a  '       -  Maximum  past  vertical  stress, 
vm 

T         -  Shear  stress  on  a  potential  failure  plane. 
T        -  Shear  stress  on  a  failure  plane  at  failure. 

t        -  Shear  stress  on  a  horizontal  plane  (direct-simple  shear 

h 


test) 
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$>   *'     -  strength  angle  in  terms  of  total  (<f>)  and  effective  (<J>') 

stresses,  respectively. 
*1  t0  ^n  ~  "tren&tn  angles  of  soils  1  to  n. 

<J>r        _  Strength  angle  required  for  equilibrium  of  a  slope. 
<J>         -  Undrained  strength  angle. 
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HIGHLIGHT  SUMMARY 

This  study  is  a  continuation  of  a  previous  one  by  Siegel  (1975) 
and  deals  with  computerized  slope  stability  analysis. 

All  investigations  are  related  to  a  computer  program,  STABL,  which 
performs  analyses  of  general  slope  stability  problems. 

A  parametric  study  was  conducted  on  the  STABL  program.  The  study 
showed  that  variations  in  the  soil  strength  parameters  c  and  <J>  have  the 
largest  influence  on  the  stability  of  a  slope. 

The  STABL  program  was  modified  to  provide  additional  options  in 
terms  of  a  new  sliding  block  surface  generator,  several  piezometric 
surfaces,  and  the  modified  Bishop  factor  of  safety  for  circular  failure 
surfaces. 

Results  of  analysis  performed  with  STABL  were  compared  to  results 
obtained  with  other  methods  of  slope  stability  analysis.  Results  by 
STABL  are  in  general  conservative,  but  the  difference  is  usually  less 
than  1055.  However,  for  deep  circular  failure  surfaces  in  soil  with 
strength  defined  mainly  in  terms  of  a  strength  intercept  c,  the  simpli- 
fied Janbu  method  of  slices  normally  applied  in  the  STABL  solution  may 
give  erroneous  and  nonconservative  results.   In  such  cases,  the  modified 
Bishop  method  should  be  used. 
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The  modified  STABL  program,  STABL2,  is  written  in  FORTRAN  IV  source 
language,  developed  for  a  CDC  6500  computer,  and  i3  also  available  in  a 
form  compatible  vith  the  IBM  360  (or  370 )  computer  system.   Both  ver- 
sions are  verified. 


I   INTRODUCTION 

This  report  presents  investigations  conducted  with  the  computer 
program  STABL,  which  performs  analyses  of  general  slope  stability- 
problems  . 

STABL  was  developed  by  R.  A.  Siegel  (1975a,  1975b)  under  a  previous 
research  program  conducted  for  and  in  cooperation  with  the  Indiana  Gtate 
Highway  Commission  (ISHC).   This  current  research  program  is  a  continu- 
ation and  extension  of  the  work  done  by  Siegel  to  develop  computerized 
slope  stability  solutions  to  meet  the  needs  of  the  ISHC  for  efficient, 
versatile  and  accurate  computational  methods. 

1.1  PURPOSE  OF  THE  PRESENT  STUDY 

The  purpose  of  the  present  study  was  originally  stated  as  follows. 

(1)  To  perform  a  parametric  study  (sensitivity  analysis)  of  the 
existing  general  program  developed  by  Siegel,  or  variations 
thereof.   This  is  needed  to  assess  the  precision  with  which 
input  data,  particularly  soil  data,  need  to  be  defined. 
High  sensitivity  to  an  input  parameter  requires  more 
attention  to  be  given  to  that  item. 

(2)  To  develop  a  hierarchy  of  special  and  simplified  programs 
to  treat  specific  and  common  slope  stability  problems.   It 
was  felt  that  the  general  program,  being  extremely  versatile, 
but  accordingly  rather  bulky  and  time  consuming,  would  be 
uneconomic  to  use  for  simple  problems. 


(3)  To  compare  the  results  from  all  programs  developed  with 
existing  ISHC  solutions.   This  will  increase  the 
confidence  in  the  new  programs ,  and  indicate  whether  they 
are  in  general  more  or  less  conservative  than  present 
solutions . 
During  the  course  of  this  research,  objective  (2)  was  somewhat 
altered.   Instead  of  developing  special  and  simplified  programs  for 
common  slope  stability  problems,  it  was  found  more  desirable  to  retain 
and  use  the  general  STABL  program. 

In  addition  to  above  objectives,  the  author  was  assigned  as 
consultant  on  computer  conversions  of  the  STABL  program  from  the  CDC  6500 
computer  at  Purdue  University,  where  the  program  is  developed,  to  the 
IBM  360  computer  at  the  ISHC  computer  center,  Indianapolis. 

1.2  ORGANIZATION  OF  REPORT 

The  report  is  arranged  as  follows. 

Chapter  II  gives  a  short  introduction  to  slope  stability  analysis 
in  general,  and  presents  some  of  the  methods  that  are  being  used.   It 
also  contains  a  brief  description  of  the  STABL  program. 

The  sensitivity  study  is  reported  in  Chapter  III  with  additional 
illustrations  in  Appendix  A. 

Modifications  of  the  STABL  program  are  described  in  Chapter  IV, 
including  an  example  of  data  setup  and  output  for  the  modified  program, 
along  with  some  general  advise  for  users. 

Chapter  V  presents  a  case  study  for  the  design  of  an  embankment 
on  varved  clay,  comparing  STABL  solutions  with  others.  Related  computer 
plots  are  contained  in  Appendix  B. 


Comparisons  of  STABL  with  existing  ISHC  procedures  for  slope 
stability  analysis  are  reported  in  Chapter  VI,  with  related  additional 
analysis  and  computer  plots  in  Appendix  C. 

Chapter  VII  summarizes  goals,  results  and  conclusions  of  the 
research,  and  gives  some  suggestions  for  future  studies  in  this  area. 

Appendix  D  describes  smaller  changes  and  errors  that  have  been 
corrected  in  the  original  STABL  program.   A  total  listing  of  the  revised 
(corrected  and  modified)  program  is  presented  in  Appendix  E. 

Finally,  Appendix  F  reports  some  of  the  experiences  in  converting 
the  STABL  program  from  the  CDC  computer  to  the  IBM  computer.   Included 
in  Appendix  F  is  a  listing  of  the  changes  necessary  to  transfer  the 
STABL  program  from  the  CDC  6500  computer  to  an  IBM  360  (or  370) 
computer  system. 

Because  of  the  large  number  of  tables  and  figures  included  in  the 
appendices,  each  of  these  are  arranged  with  the  text  (if  any)  first, 
then  the  tables  in  consecutive  order,  followed  by  the  figures  in  con- 
secutive order.   This  was  done  to  avoid  interrupting  the  text,  and 
because  most  of  the  text  readily  can  be  interpreted  apart  from  the 
tables  and  figures. 


II   SLOPE  STABILITY  ANALYSIS 

2.1  GENERAL  DESCRIPTION 

Earthen  masses  located  beneath  a  sloping  ground  surface,  whether 
natural  or  man-made ,  have  a  tendency  to  move  downward  and  outward 
under  the  influence  of  gravity.   Unless  this  tendency  is  suitably 
counteracted  by  the  shearing  resistances  within  the  mass  ,  a  landslide 
occurs.   Other  forces  such  as  seepage  forces  and  surcharge  loads  may 
add  to  the  instability.   Avoiding  such  instabilities  is  a  major  concern 
of  the  geotechnical  engineer. 

The  investigation  of  the  stability  of  a  slope  has  three  main 
steps : 

(1)  Identification  of  the  forces  which  cause  and  resist  failure. 

(2)  Determination  of  the  shear  strengths  of  the  soils  beneath  the 
slope . 

(3)  Application  of  a  model  to  investigate  the  stability  based  on 
(1)  and  (2). 

The  forces  acting  on  the  slope  depend  on  the  profile  and  density  of 
the  soils;  the  position  of  the  phreatic  surface,  and  associated  seepage 
forces ;  and  external  applied  loads . 

The  shear  strength  of  a  particular  soil  depends  upon  its  type  and 
origin,  water  content,  prior  stress  history,  method  and  type  of 
compaction,  and  the  like. 


The  shear  strength  of  a  soil  was   initially  expressed  by  Coulomb 
(1776)   as 

Tf  =   c  +  of  tan  <J>  (2.1) 

where   c  is   the  strength  intercept,   <j>  is   the  strength  angle,    and  a     is 
the  stress   normal  to  the  shear  plane. 

The  equation  commonly  used  today  is  modified  to  express  T  in 
terms  of  effective  stresses  rather  than  total  stresses  and  has  the 
form 

Tf  =   c'    +   (a     -  u)    tan  <J>' 

or  t     =  c"   +  a*      tan  <f>'  (2.2) 

where   c'    and  <J)'    are  the  strength   intercept  and  strength   angle  with 
respect  to  effective  stresses,   also  named  effective  strength  parameters, 
and  a*      =   (o     -  u)   is   the  effective  stress  normal  to  the  shear  plane. 

In  effective  stress   analysis   of  sands   and  gravels,    c'    ~  0,    and 
equation    (2.2)    is   simplified  accordingly.      In  total  stress   analysis, 
and  with  particular  reference  to  un drained  shear  in   clayey  soils,   <J>  =  0, 
and  equation   (2.1)   is   reduced  to  Tf  =   c. 

Instability   of  a  slope  in  sands   and  gravels   occurs   if  the  slope 
inclination  exceeds   a  certain  value,  known  as   the   angle  of  repose.      This 
angle  is    constant   for  a  particular  material  at   a  particular  void  ratio 
and  is   independent   of  the  height   of  the  slope,   provided  the  slope  is  not 
subjected  to  seepage   forces.      In   contrast,   a  slope   in  a  clayey  soil  has   a 
critical  height,  which   depends   on  the  slope  inclination  and  the  soil 
properties . 
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The  shearing  behaviour  of  a  particular  material  ir,   determined 
empirically,    i.e.,  by   field  tests   or  laboratory  tests.     Most  common  are 
the  laboratory  triaxial   and  direct  shear  tests   and  the   field  vane  test. 

The   results   are  then  applied  to  the  slope  stability  analysis  by 
assuming  certain  mechanisms   and  modes   of  failure. 

A   common  approach  to  slope  stability  analysis   is   to  investigate 
failure   along  an  assumed  slip  surface,   and  express   the  stability   in  terms 
of  a  safety   factor  F,    defined  as   the   ratio  of  the  available  shearing 
forces   along  the  slip  surface,  S    ,   to  the  shearing   forces   required  for 
static  equilibrium,   S    : 

S 

F  =  <T 

r  (2.3) 

An  F  value  equal  to  or  below  unity   indicates   that  the  slope   is 
unstable   and  failure   can  be  expected. 

To  identify  the  more   critical  slip  surfaces    (lower  F  values)    it 
is   generally  necessary  to   conduct  the  analysis   for  a  number  of  surfaces 
of  different   shape  and  location. 

2.2     DIFFERENT  METHODS  OF  SLOPE  STABILITY  ANALYSIS 

2.2.1     Early  Development 

Until  Coulomb    (1776)    [Skempton   (19^6) ]   soils  were  considered  as 
strictly   frictional,    i.e.,   a  slope  was   stable  if  the  slope  angle  was 
less  than  the   angle  of  repose  of  the  material.      Coulomb   analysed  the 
stability  of  a  vertical  bank  of  cohesive  soil  with  respect  to  a  plane 
slip  surface,    and  obtained  the  equation 

H     =  k     [-    cos.^  J  (2.10 

c         y       1  -  sin  <f> 


where  H     is  the   critical  height,   c   is   the  strength   intercept,   <J>  is   the 
strength   angle,    and  y  is   the  unit  weight   of  the  soil. 

In  1820  Francais   gave  a  solution  to  the  more  general  problem  of  a 
sloping  bank.      He   also  assumed  a  plane  slip  surface   and  found 


ti     _  if£     r      cos   4>  sin  B      ■<  (o  k\ 

Hc  "     Y      [1  -  cos    (8  -  <fr)]  (2'5) 


where  0  is   the  slope  inclination. 

Collin   (181+6)   introduced,  based  on  his   observations   of  deep  curved 
slip  surfaces   in   clay  slopes,   an   analysis   assuming  a  cycloidal  slip 
surface  and  a  <J>  value  of  zero. 

However,   not  until  the  twentieth   century   did  greatly  improved 
measurements   of  shear  strength   and  refined  procedures   of  analyses 
appear. 

Wright    (I969)   and  Siegel   (1975a)   have  described  in  detail  the 
different  methods   of  slope  stability  analysis   available  today.      As   this 
study  is   a  continuation  of  that  begun  by  Siegel,   only  a  resume  of  the 
most   common  methods  will  be  given  here . 

2.2.2     The  <ft  =  0  Method 

Based  on   field  observations   in  Sweden,   a  circular  slip  surface  was 
introduced  by  Petterson   (1916)    [Skempton   (l9*+6)],  which   combined  with 
the   assumption  that  <j>  equals   zero,   gave  what  we  today  know  as  the  <J>  =  0 
method.      Figure  2.1  shows  the  principle.     When  $  equals    zero,   the  shear 
strength  T_  equals   c   (Equation  2.1)    and  is   independent  of  the  actual 
value  of  the  normal  stress   at   any  point   of  the   failure  surface.      Equilib- 
rium of  tte  potential  sliding  mass   is   satisfied  by  moment   equilibrium 


L. 

u 

h 


3 

S 

i    5 


o 
o 

I 
I- 


o 

H 


Id 

I 
h 

o 

id 

_j 
cl 

o 

z 

a: 


CM 
Id 

u. 


around  the  center  of  the  circle.   The  forces  acting  normal  to  the 
failure  surface  pass  through  the  center  of  the  circle  and  do  not  con- 
tribute to  the  moment,  hence  the  shear  stress  required  for  equilibrium 
(c  )  is  determined  directly  by  the  moment  equation,  and  the  factor  of 
safety  F  can  be  found  from  Equation  2.3.  Different  centers  and  radii  are 

investigated  to  determine  F  .  .  Taylor  (1937)  made  a  chart  that  directly 

min 

c 
gives   the  stability  number  — —     as    a  function  of  the  slope   angle   3   and 

r  yli 

depth  factor  n.   Knowing  c,  y   and  H,  F  can  be  found.   Gee  Figure  ^.7- 

The  <}>  =  0  method  is  commonly  used  for  short  term  stability  analysis 
of  saturated  clay  slopes,  where  a  failure  would  occur  undrained. 

2.2.3  The  Logarithmic  Spiral  Procedure 

When  4> '(<!>)  is  not  taken  as  zero,  t  is  not  a  constant,  but  depends  on 
the  normal  stress  O'(a)  along  the  potential  failure  surface,  and  a  dis- 
tribution of  0'{a)   has  to  be  assumed  in  order  to  determine  F.  To  over- 
come this  problem  the  logarithmic  spiral  was  introduced  as  a  proposed 
failure  surface.   Such  a  surface  has  the  characteristic  that  the  resultants 
of  the  normal  stresses  0'  and  frictional  components  of  shear  strength 
a'  tan  4>'  pass  through  the  center  of  the  spiral.   Consequently  moment 
equilibrium  around  the  center  will  only  involve  the  body  forces  and 
strength  intercept  force  (C)  of  the  soil  as  in  the  <j)  =  0  method.   See 
Figure  2.2. 

The  equation  of  a  logarithmic  spiral  can  be  expressed  in  polar 
coordinates  as 


etanO'/FJ  (2.6) 


r  =  r     e  <j> 

o 
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where  r  is  the  radius  from  the  center,  r  is  the  initial  radius,  0  is 

o 

the  angle  between  initial  radius  r  and  radius  r,  (J)1  is  the  effective 

strength  angle  and  F  is  the  factor  of  safety  with  respect  to  the  strength 

angle.   If  F  is  taken  as  unity  [Taylor  (1937)]  the  procedure  directly 

gives  F  ,  the  factor  of  safety  with  respect  to  the  strength  intercept,  c'. 

In  most  cases,  however,  one  will  require  an  F  larger  than  unity  for  both 

<J>'  and  c'.   Assuming  F ,  =  F  ,  several  spirals  with  the  same  center  and 
©    c 

initial  radius  r  are  investigated  to  find  the  value  of  F.  To  find  F  .  , 
o  min 

the  center  0,  and  initial  radius  r  have  to  be  varied. 

o 

2.2.U  The  Friction  Circle  Method 

For  a  circular  shear  surface  the  resultants  of  the  normal  stresses 
and  the  frictional  components  of  shear  strength  will  lie  on  lines  tangent 
to  a  circle  of  radius  r  sin  <}>',  called  the  friction  circle.   Hence  the 
moment  about  the  center  of  the  circle  will  involve  the  normal  stress  dis- 
tribution which  is  unknown.   The  assumption  commonly  made  is  that  the 
resultant  of  all  normal  stresses  and  frictional  shear  stresses  also  lies 
tangent  to  the  friction  circle.   This  assumption  is  equivalent  to  assuming 
that  the  normal  stresses  are  concentrated  at  a  single  point  [Wright  (1969)] 
and  the  solution  obtained  is  a  lower  bound  solution  (conservative).  The 
principle  is  shown  in  Figure  2.3.   Taylor  (1937)  also  presented  charts 

for  the  stability  number  =-rr  based  on  the  friction  circle  method.  See 

J  FyH 

Figure  U.6. 

2.2.5     Method  of  Slices 

The  methods   described  thus   far  are  particularly  suited  to  simple 
slope  stability  problems  with  homogeneous   soil   conditions   and  no   forces 
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other  than  those  due  to  weight  and  shearing  resistance.  For  problems 
involving  changing  soil  conditions  and  additional  forces  such  as  seepage 
and  earthquake  forces,  several  limiting  equilibrium  methods  have  been 
developed,  all  based  on  the  method  of  slices. 

The  soil  mass  above  an  assumed  slip  surface  is  divided  into  slices 
in  such  a  way  that  each  slice  has  a  minimum  change  in  variables  (Figure 
2 ,U) .  The  system  of  equations  and  unknowns  associated  with  complete 
equilibrium  of  the  entire  soil  mass  in  terms  of  n  slices  is  summarized 
in  Table  2.1  (see  also  Figure  2.5).   For  any  value  of  n,  other  than  unity, 
there  are  (2n  -  2)  more  unknowns  than  equations,  and  the  problem  is 
indeterminate.  The  differences  in  the  various  methods  of  stability 
analyses  are  contained  in  the  different  assumptions  made  with  respect 
to  these  unknowns  in  order  to  render  the  problem  determinate. 
2.2.5.1  Ordinary  Method  of  Slices 

This  method  assumes  a  circular  slip  surface.  The  assumptions  made 
by  Fellenius  (1926)  to  eliminate  the  unknowns  are 

(1)  The  normal  force  N  on  the  base  of  each  slice  acts  at  the 
center  of  the  base. 

(2)  The  interslice  forces  are  zero. 

These  two  assumptions  eliminate  (^n  -  3)  unknowns,  which  is  more  than 
is  necessary  to  make  the  problem  statically  determinate.  Remaining  are 
(n  +  l)  unknowns,  namely  F  and  the  n  normal  forces  at  the  base  of  each 
slice.   The  equations  used  are:   (a)  force  equilibrium  normal  to  the 
base  of  each  slice,  which  directly  gives  the  normal  forces  N,  and  (b) 
overall  moment  equilibrium  around  the  center  of  the  circle  forming 
the  slip  surface,  which  directly  gives  the  shear  strength  required  for 
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Table  2.1 
Equations  and  Unknowns  Associated  with  the  Method  of  Slices 

Equations 

n      Moment  equilibrium  for  each  slice. 
2n       Force  equilibrium  in  2  directions  for  each  slice. 
n      Mohr-Coulomb  relationship  between  shear  strength  and 
normal  effective  stress. 


l»n      Total  number  of  equations. 


Unknowns 


1  Factor  of  Safety  F. 

n  Normal  force  at  base  of  each  slice  N. 

n  Position  of  N. 

n  Shear  force  at  base  of  each  slice  S. 

n  -  1  Horizontal  interslice  forces  E. 

n  -  1  Vertical  interslice  forces  X. 

n  -  1  Location  of  interslice  forces  (line  of  thrust) 


6n  -  2     Total  number  of  unknowns. 


*   n  is  the  number  of  slices. 
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FIGURE  2.5    FORCES   ACTING  ON  A  SLICE 
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equilibrium.     The   factor  of  safety  F  will  be  expressed  in  terms   of  an 
average  value   (summing  over  the   n  slices)    of  the  shear  strength   avail- 
able divided  by  the   shear  strength   required. 

Because  the  problem  is   overdetermined  the  values   obtained  for  the 
normal   forces  will  usually  not  satisfy  equilibrium  in   directions   other 
than  those  normal  and  parallel  to  the  base  of  each  slice.      Hence   the 
neglect   of  interslice   forces  may  lead  to  unrealistic  results. 
2.2.5.2     Bishop's  Method 

Bishop   (1955)    introduced  a  method  satisfying  complete  equilibrium. 
His   primary  equations   are  no.    of  eq. 

(1)  Overall  moment   equilibrium  1 

(2)  Vertical   force  equilibrium  for  each  slice  n 

(3)  Force  equilibrium  parallel  to  base  of  each  slice  n 
(h)      Shear   force   as   function   of  normal   force  and  F                    n 

Total  of  3n  +   1 

In  addition  he  has   two  equations  with   respect   to  overall  equilibrium 
of  interslice   forces   X  and  E. 

In   the  primary   analysis   Bishop   does   not  make  any  assumptions 
about  the  line  of  thrust   (see  Figure  2.5)-      This   line  may  be   determined 
by   adding  the   remaining  (n  -  l)   moment   equations    for  the   individual 
slices,  but   the  primary  equations   can  be  solved  independent   of  the 
location  of  the  line  of  thrust.      The  normal   force,   as  well  as   the 
resultant  weight   of  each  slice,   is   assumed  to  act  through   the   center 
of  the  base  of  the  slice.     The   remaining  unknowns   are 
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no.  of  eq. 

(1)  Factor  of  safety  F  1 

(2)  Normal  force  at  base  of  each  slice  V.  n 

(3)  Shear  force  at  base  of  each  slice  n 
(1+)  Vertical  interslice  forces  X                        n  -  1 
(5)   Horizontal  interslice  forces  E                      n  -  1 

Total  of       Un  -  1 
This  total  of  unknowns  is  compared  to  ( 3n  +  l)  equations,  which  gives 
an  indeterminacy  of  (n  -  2).   The  procedure  is  to  assume  values  of  X 
that  satisfy  the  two  overall  interslice  force  equations  (E  can  be  expressed 
in  terms  of  X  from  the  force  equation  parallel  to  the  base  of  a  slice)  . 
To  judge  whether  a  particular  solution  is  reasonable,  one  may  calculate 
the  corresponding  line  of  thrust  using  the  remaining  (n  -  l)  moment 
equations  for  the  individual  slices. 
2.2.5.3  Modified  Bishop  Method 

The  simplest  solution  to  Bishop's  method  is  to  assume  X  equal  to 
zero  which  directly  satisfies  equilibrium  with  respect  to  the  vertical 
interslice  forces,  but  generally  not  equilibrium  with  respect  to  the 
horizontal  interslice  forces.   The  solution  is  called  the  modified  or 
simplified  Bishop  method.   It  can  be  directly  arrived  at  by  neglecting 
the  interslice  forces  as  for  the  ordinary  method  of  slices,  but  using 
vertical  force  equilibrium  for  each  slice  instead  of  force  equilibrium 
normal  to  the  base  of  each  slice. 
2.2.5.^  Other  Procedures 

Many  other  methods  have  been  developed  by  different  investigators 
making  various  assumptions  with  respect  to  the  interslice  forces. 
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Spencer  (1967)  assumed  parallel  interslice  forces  and  developed  a 
statically  determinate  solution  with  3n  equations  and  3n  unknowns. 
Morgenstern  and  Price  (1965)  assumed  a  relationship  between  the 
horizontal  and  vertical  interslice  forcer,  dependent  on  a  function  f(x) 
along  the  shear  surface.   They  extended  their  solution  to  failure 
surfaces  of  general  shape. 

Janbu  (195^+a)  formulated  an  expression  for  F  in  terms  of  vertical 
force  equilibrium  for  each  slice  combined  with  overall  horizontal  force 
equilibrium.   As  with  the  Bishop  method,  this  procedure  requires 
assumptions  of  the  vertical  interslice  forces  X.   Janbu  made  his 
assumptions  so  that  they  satisfy  interslice  force  equilibrium  as  well 
as  moment  equilibrium  for  each  slice,  assuming  a  position  of  the  line 
of  thrust.  A  simplified  procedure  results  if  X  is  assumed  equal  to 
zero.   This  procedure  is  similar  to  the  modified  Bishop  method,  the 
difference  being  that  overall  horizontal  force  equilibrium  is  used 
versus  the  use  of  overall  moment  equilibrium  in  the  modified  Bishop 
procedure. 

Carter  (1971)  reached  the  same  expression  as  Janbu's  simplified 
procedure,  but  arrived  at  it  in  a  different  manner  by  considering  the 
modified  Bishop  procedure  applied  to  slip  surfaces  of  general  shape. 
In  the  case  where  the  slip  surface  is  not  a  circular  arc,  no  specific 
point  exists  that  is  obvious  to  use  as  moment  center.   Carter  found 
that  the  value  of  the  factor  of  safety  was  independent  of  the  x-coordinate 
of  the  moment  center,  but  varied  with  the  y-coordinate,  and  reached  a 
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# 
minimum  value  for  y  =  °°.    lie  concluded  that  y  =  °°  was  the  only  consis- 
tent moment  center  to  use  for  an  irregular  as  well  as  a  circular  slip 
surface. 

2.2.6  Methods  other  than  Limiting  Equilibrium 

An  alternative  method  of  slope  stability  analysis  is  to  investigate 
the  shear  stresses  imposed  under  the  slope  by  means  of  the  theory  of 
elasticity  [Perloff  and  Baron  (1976),  Romani  (1970),  Romani,  Lovell  and 
Harr  (1972)].  The  criterion  is  that  the  shear  stress  at  any  point  does 
not  exceed  the  shear  strength.  The  factor  of  safety  is  defined  as  the 
shear  strength  divided  by  the  shear  stress  at  the  point  where  this  ratio 
is  the  least,  hence  it  gives  the  safety  at  the  most  critical  point,  whereas 
F  defined  by  a  limiting  equilibrium  method  gives  the  average  value  of 
safety  over  the  entire  slip  surface. 

The  method  may  be  useful  when  dealing  with  soils  where  progressive 
failure  is  likely  to  occur,  because  it  will  discover  the  areas  of  stress 
concentration  from  which  a  failure  is  likely  to  initiate.   On  the  other 
hand,  it  does  not  take  into  account  the  redistribution  of  stress  which 
occurs,  when  the  stress  level  at  a  point  approaches  the  strength. 

In  a  case  with  concentrated  surface  loads ,  methods  based  on  an 
elastic  model  can  take  into  account  the  distributional  effect  of  the 
slope  material.  This  cannot  be  done  with  the  method  of  slices  in  the 
limiting  equilibrium  models,  and  an  error  may  be  introduced  by  applying 


*Further  investigations  by  the  author  show  that  this  is  not  always  true, 
see  Chapter  VI  and  Appendix  C.II. 
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such  a  model  in  this  case.   For  important  and  complex  problems  it  may- 
be advisable  to  perform  both  types  of  analysis. 

2.2.7  Probabilistic  Approach 

All  the  methods  presented  thus  far  have  been  discussed  from  a 
deterministic  viewpoint,  as  if  the  parameters  involved  in  the  slope 
stability  analysis  were  single  fixed  values.  This  is  a  gross  simplifi- 
cation of  the  real  situation,  where  uncertainties  exist  with  respect  to 
soil  properties  as  well  as  the  soil  profile. 

Several  people,  e.g.,  Wu  et  al  (1970),  Athanasiou-Grivas  (1976), 
and  Catalan  and  Cornell  (1976),  have  worked  on  probabilistic  models  for 
slope  stability  analysis,  which  make  use  of  slip  surface  and  limiting 
equilibrium  models,  but  which  operate  in  terms  of  probability  of  failure 
rather  than  the  traditional  factor  of  safety.   In  this  approach  the  soil 
properties  are  defined,  not  as  exact  values,  but  as  distributions  with 
a  mean  value  and  standard  deviation. 

The  relationship  between  the  factor  of  safety  as  defined  by  the 
traditional  methods,  and  the  probability  of  failure  was  investigated  in 
a  case  study  by  Athanasiou-Grivas  (1976)  ,  where  a  slope  having  a  con- 
ventional factor  of  safety  of  1.2  failed.  The  probability  of  failure 
was  found  to  be  approximately  0.3%  with  zero  pore  pressure  and  1%   for  a 
pore  pressure  parameter  (r  =  u/yh)  equal  to  0.6. 


*The  significance  of  varying  soil  properties  may  be  investigated 
within  the  deterministic  model  by  a  sensitivity  study  (see 
Chapter  III) 
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2.2.8  Comparison  of  Methods 

Wright  (1969)  has  done  a  comprehensive  study  of  the  many  methods 
available  for  slope  stability  analysis.   He  found,  in  agreement  with 
others  [Bishop  (1955),  Lowe  and  Karafiath  (1967)],  that  for  homogeneous 
soil  the  modified  Bishop  procedure  is  as  accurate  as  any  of  the  more 
rigorous  methods  (Spencer,  Janbu,  Morgenstern  and  Price),  whereas  the 
ordinary  method  of  slices,  as  well  as  Janbu' s  simplified  procedure,  gave 
factors  of  safety  as  much  as  10$  and  5%   lower  than  the  lower  bound  value 
given  by  the  friction  circle  method.   Some  of  the  results  are  shown  in 
Table  2.2.  The  difference  is  larger  if  pore  pressures  are  present,  and 
for  cases  with  non-homogeneous  soil  it  may  be  even  larger. 

In  any  case  where  the  slip  surface  can  be  approximated  by  a  circular 
arc,  Bishop's  modified  procedure  is  recommended,  as  it  combines  reason- 
able accuracy  with  simplicity  and  small  computer  time  requirements. 

If  the  slip  surface  cannot  be  approximated  by  a  circular  arc,  Bishop's 
method  cannot  be  used  without  modifications.   The  method  that  is  most 
convenient  for  irregular  slip  surfaces  is  undoubtedly  Janbu' s  simplified 
procedure.   Unfortunately  it  is  not  so  accurate  and  may  lead  to  overly 
conservative  (and  hence  uneconomic)  design.   If  this  is  a  problem, 
Spencer's  procedure  seems  to  be  the  more  rigorous  method  that  is  the 
least  time  consuming.   However,  it  requires  analysis  of  the  position  of 
the  line  of  thrust,  together  with  soil  properties,  to  judge  whether  the 
result  is  conservative  or  not  (see  example  in  Appendix  C.l). 

A  study  comparing  the  modified  Bishop  method  with  a  linear  elastic 
finite  element  method  is  reported  by  Wright  et  al  (1973)  with  the 
following  results: 
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Table  2.2 

Factor  of  Safety  for  Different  Methods  of  Slope  Stability  Analysis 

cot  6 

Method      X                1.5/1  2.5/1  3.5/1 

1  2               1.171  1.1*85  1.768 

2  1.125  1.1»02  1.659 

3  l.lVf  1.M3  1.711* 
h                                                       1.103  1.351  1.595 

1  5              1.831  2.U80  3.095 

2  1.738  2.331*  2.911* 

3  1.785  2.1*09  3.008 
h                                                    1.703  2.277  2.850 

1  8                                    2.U30  3.399  1».333 

2  2.301  3.210  U.105 

3  2.368  3.309  k.226 
k                                                       2.262  3.151*  't.0U5 

1  20                                    2.309  3MU  U.503 

2  2.197  3.268  U .332 

3  2.257  3.3'+7  1*.1425 
U                                                  2.176  3.2U2  U.307 

1  50                                    1.93U  2.988  1+.035 

2  1.863  2.901  3.936 

3  1.902  2.9^8  3.991 
h                                                       1.855  2.893  3.928 

X  ~   c/yh  '  P  =  slope  ansle' 

F 
The  stability  number  N  =  ■■>■        is  a  function  of  X  and  B  only, 

and  is  constant  for  constant  X  and  3  although  c,  $,  y   and  h  may 

vary. 

The  methods  compared  are  (l)  Simplified  (or  Modified)  Bishop 

Method,  (2)  Ordinary  Method  of  Slices,  (3)  Friction  Circle  and 

(h)   Janbu's  simplified  procedure. 
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(1)  The  distribution  of  normal  stress  along  the  failure  surface 
was  approximately  identical  for  the  two  methods. 

(2)  Although  the  factor  of  safety  by  the  finite  element  method 
varies  along  the  failure  surface,  the  average  value  is  nearly 
the  same  as  for  the  modified  Bishop  method. 

(3)  Because  the  finite  element  method  involves  a  completely 
different  set  of  assumptions  than  those  employed  in  the 
modified  Bishop  method,  the  close  agreement  between  the  two 
methods  indicates  that  they  both  are  reliable.   Because  the 
limiting  equilibrium  model  is  much  simpler  than  a  finite 
element  analysis,  the  former  will  usually  be  preferred. 

2.3  STABL  -  A  COMPUTER  PROGRAM  FOR  GENERAL  SLOPE  STABILITY  ANALYSIS 

Many  programs  exist,  based  on  one  of  the  previous  described  methods 
of  slope  stability  analysis.   When  STABL  was  developed,  the  major  purpose 
was  to  develop  a  program  with  the  least  number  of  restrictions  imposed  on 
the  input  data,  and  the  maximum  utilization  of  the  computer  in  the  search 
for  the  most  critical  failure  surfaces,  combined  with  ease  of  use. 

STABL  can  handle  almost  any  problem  with  respect  to  boundary  geometry 
and  loading,  soil  conditions  including  anisotropic  soils,  equilibrium 
or  excess  of  equilibrium  pore  pressures  and  pseudo-static  earthquake 
forces.   It  has  routines  that  search  for  the  more  critical  failure  sur- 
faces by  generating  surfaces  of  circular,  sliding  block  or  irregular 
shape  and  compare  the  corresponding  safety  factors.   It  can  as  well 
handle  specified  failure  surfaces.   It  is  limited  to  2-dimensional  analy- 
sis,  which  is  the  case  with  almost  all  methods  available  today. 


*For  an  evaluation  of  3-dimensional  effects  on  slope  stability  analysis, 
see  Baligh  and  Azzouz  (1975)- 
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The  method  chosen   for  the  calculation  of  the   factor  of  safety   is 
Janbu's   simplified  procedure  satisfying  vertical   force   equilibrium  for 
each  slice   and  overall  horizontal   force  equilibrium,   assuming  horizontal 
interslice   forces.      For  the   simple   case  of  dry  soil  and  no  external 
loading,   the  equations   are  as    follows: 

Vertical   force  equilibrium  for  a  slice 

N  cos  a  +  S  sin  a  -  W  =  0  (2.7) 

Overall  horizontal  force  equilibrium 

£   (U  sin  a  -  S  cos  a)   =  0  (2.8) 

Mohr-Coulomb   relationship 


C    +  N'    tan  d>' 


(2.9) 


F 
From   (2.7)    and  (2.9)      (N  =  N1) 

N  =  W/cos   a  -  S   tan  a  (2.10) 

_  C    +  W  tan  (ji'/cos   a  ,  x 

F  +  tan  a  tan  (J)'  K£-u.i 

Substitute   (2.10)   into   (2.8) 

£   (W  tan  a  -  S/cos  a)   =0  (2.12) 

Substitute    (2.11)    into   (2.12) 

_   rif   .  C*   +  W  tan  <j)'/cos   a -i   _  n 

Z    [W  tan  a  -  - — ■ r=,  -   7 - tttJ        ° 

cos  a  (F  +  tan  a  tan  <p') 

Rearranging  gives 

j,  W  cos   g   (F  tan  a  -  tan   4>')    -  C    _  Q 
cos  a  (F  +  tan  a  tan  <t>') 


„  C'/cos  a  +  W  (tan  <j>'    -  F  tan  a)   _  Q  (2  13\ 

F  +  tan  a  tan  (b1 
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This  expression  is  the  same  as  that  derived  by  Carter  [Siegel  (1975a)] 
A  similar  expression  including  surface  loads,  water  table  and  earth- 
quake forces  is  presented  by  Siegel  (1975a);  see  also  Chapter  IV, 
Equation  (*+.l). 

This  method  was  chosen  for  its  simplicity  and  because  it  can  be 
applied  to  failure  surfaces  of  any  shape.   It  is  usually  conservative 
and  may  in  some  cases  give  factors  of  safety  which  are  too  low.   A  com- 
parison was  made  with  the  modified  Bishop  method  from  Wright  (1969). 
STABL  gave  values  of  F  from  2.3%  to  8.8%  lower  than  those  obtained  by 
the  modified  Bishop  method,  the  largest  difference  occurring  for  low 
values  of  X   (see  Table  2.3).   Siegel  compared  STABL  with  Spencer's  pro- 
cedure and  found  a  maximum  difference  of  15%  for  steep  slopes  and  high 
values  of  A,  with  a  pore  pressure  parameter  r  =  0.5,  STABL  giving  the 
lower  value. 

In  most  slope  stability  problems,  however,  there  is  an  uncertainty 
in  the  determination  of  the  input  data  such  as  the  soil  profile  and 
particularly  the  soil  properties  in  terms  of  density  and  strength.  The 
factor  of  safety  may  be  more  sensitive  to  variations  in  these  properties 
than  to  different  methods  of  analysis  [Singh  (1970)].   This  point  is 
further  investigated  in  Chapter  III. 

A  detailed  discussion  of  how  the  program  performs  its  functions  is 
given  by  Siegel  (1975a).  Later  modifications  of  the  program  are  reported 
in  Chapter  IV  of  this  report. 

^Important  exceptions,  see  Chapter  VI,  Problems  1  and  k,    and  Appendix  C.II. 
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Table  2.3 

Comparison  of  STABL  with  Modified  Bishop  Method  after  Wright  (1969), 
Using  Failure  Surfaces  Obtained  by  the  Modified  Bishop  Analysis 

cot  6 

Factor  of  Safety    A          1.5/1  2.5/1  3.5/1 

Mod.  Bishop        2          1.171  1.1*85  1.768 

STABL                        1.106  1.367  1.625 

%   difference                 5-9  8.6  8.8 

Mod.  Bishop        5          1.831  2.1*80  3.095 

STABL                        1.720  2.318  2.890 

%   difference                 6.5  7.0  7.1 

Mod.  Bishop        8          2.1*30  3-399  1+.333 

STABL                        2.290  3.209  1+.099 

%   difference                 6.1  5-9  5.7 

Mod.  Bishop       20          2.309  3 -UlU  1+.503 

STABL                         2. 201*  3.275  U.3U6 

%   difference                 1+.8  U.2  3-6 

Mod.  Bishop       50          1.931+  2.988  1*.035 

STABL                        1.873  2.912  3.91+5 

%   difference                 3.3  2.6  2.3 

X  =  *?*  4*  ,   6  =  slope  angle  (see  also  Table  2.2). 
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III   SENSITIVITY  ANALYSIS 

3.1  INTRODUCTION 

A  sensitivity  study  was  performed  on  the  STABL  program  to  assess 
the  precision  with  which  input  data  need  to  be  defined,  particularly 
soil  data.   Although  the  particular  results  apply  to  this  program,  the 
trend  is  believed  to  be  applicable  to  slope  stability  analysis  in 
general. 

Little  research  has  been  published  on  sensitivity  analysis.   An 
article  by  Singh  (1970)  concludes  that  the  factor  of  safety  for  slope 
stability  analysis  is  much  more  sensitive  to  the  choice  of  strength 
parameters  as  interpreted  from  soil  tests  than  to  the  choice  of  the 
computational  method  of  analysis.   This  conclusion  is  of  particular 
interest  with  respect  to  the  STABL  program,  because  the  procedure  of 
analysis  used  in  this  program  usually  is  somewhat  conservative  compared 
to  other  methods  of  analysis  (see  Chapters  II,  V  and  VI). 

3.2  PROCEDURE  FOR  ANALYSIS 

The  procedure  initially  chosen  for  the  sensitivity  analysis  is  as 
follows  : 

(l)  For  a  given  slope  (geometry,  soil  profile,  water  table,  seep- 
age, load  conditions),  a  STABL  searching  routine  is  used  to 
generate  a  failure  surface,  which  will  be  close  to  the  criti- 
cal one,  and  the  safety  factor  is  calculated. 
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(2)  Using  this  specified  surface,  the  factor  of  safety  is  calcu- 
lated for  a  small  change  in  an  input  parameter  equal  to  the 
expected  variation  in  the  field  or  scatter  in  test  results. 

(3)  The  influence  of  the  change  in  the  factor  of  safety  due  to  a 
change  in  one  of  the  input  parameters,  relative  to  the  total 
change  in  the  factor  of  safety  due  to  changes  in  all  parame- 
ters, is  investigated. 

(k)     A  conclusion  may  be  reached  with  respect  to  the  required  pre- 
cision in  definition  of  the  different  variables  to  limit 
uncertainties  in  the  factor  of  safety  to  tolerable  levels. 

The  steps  are  illustrated  on  the  flow  chart  of  Figure  3.1. 

3.3  CASES  STUDIED 

The  sensitivity  analysis  was  performed  for  the  idealized  case  of 
a  slope  in  homogeneous  soil  with  no  pore  water  pressure,  and  for  an 
embankment  founded  on  a  soil,  which  has  a  weak  layer. 

3.3.1  Sensitivity  Analysis  of  homogeneous  Slopes 

It  was  decided  to  perform  the  sensitivity  study  for  a  relatively 
wide  range  in  parameters  for  slopes  in  homogeneous  soils,  to  get  a  good 
feel  for  how  the  factor  of  safety  is  likely  to  change,  when  different 
parameters  are  changed. 

For  homogeneous  soils  the  factor  of  safety  depends  on  5  to  6  inde- 
pendent (to  some  extent)  parameters,  which  are: 

Slope  geometry  defined  by 

Height  h 

Slope  angle  6 


STEP  I 
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Slope  geometry 
Soil  characteristics 

Water  conditions 
Loads 


STABL 


Locate 
critical 
failure 
surface 


STEP  2 


Variation  in 
problem  variables 


'STABL     w 
critical    surface 


Variation   in 
Factor  of  Safety 


STEP  3 


AF(Axi)  vs.  E      AF(Axj) 


STEP  4 


Restrictions  on  Ax;     to  limit     AF(Ax:) 


FIGURE    3.1        FLOW  CHART  SHOWING  STEPS  IN 
SENSITIVITY    ANALYSIS 
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Soil  properties  defined  by- 
Unit  weight  Y 
Strength  intercept         c' 
Strength  angle  (J> ' 
Pore  water  pressure,  if  present,  is  important,  as  well.   This  parameter 
was  neglected  in  this  study. 

Figure  3.2  shows  the  range  in  these  parameters,  covered  by  the 
sensitivity  analysis.  Five  different  slopes  were  examined.  Height  (h) 
and  unit  weight  (y)  were  chosen  arbitrary  as  100  feet  and  100  pcf,  re- 
spectively.  However,  the  results  can  be  applied  to  any  values  of  h  and 

c ' 
Y  by  introducing  the  dimensionless  number  —  [ Taylor ( 1937) ] .  The  slopes 

were  analysed  for  all  combinations  of  (c',  <f>')  shown  in  Figure  3.2,  with 
the  exception  of  those  where  a  combination  of  low  strength  intercept  c  * , 
low  strength  angle  <f>' ,  and  a  steep  slope,  gave  values  of  the  factor  of 
safety  much  below  unity. 
3.3.1.1  Analysis  by  STABL 

For  each  problem  case  (slope  geometry,  soil  properties)  limitations 
for  initiation  and  termination  of  failure  surfaces  were  set  according 
to  the  particular  problem,  e.g.,  all  failure  surfaces  for  slopes  with  a 
slope  angle  8  greater  than  53°  passed  through  the  toe;  for  6  less  than 
53°  and  <j>  =  0,  always  a  base  failure.  A  hundred  surfaces  were  generated 
for  each  problem.   If  the  critical  surfaces  were  located  on  the  boundary 
of  the  search  limitations,  these  were  changed  accordingly,  and  a  new 
search  was  carried  out.  An  exception  to  this  procedure  was  the  <J>  =  0 
case,  where  the  critical  failure  surface  theoretically  extends  to  an 
infinite  depth.  Only  the  CIRCLE  search  was  used  because  of  the  homoge- 
neity of  the  soil  (see  also  p.  69). 
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Vertical  Slope 


Slope  1/2 


Slope  l/l 


Slope   2/1 


Slope   3/1 


Height  of  slopes 
Soil  characteristics; 
unit  weight 
strength  intercept 
strength  angle 


h  *■  100  feet 

y*  100  pcf 

c'- 1000-3000  psf 

<£«  0°-50° 


FIGURE  3.2        RANGE  IN  VARIABLES  FOR  SENSITIVITY 
ANALYSIS  OF  HOMOGENEOUS  SLOPES 
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The  most  critical  failure  surfaces  generated  by  STABL  for  each  par- 
ticular problem  were  plotted  for  comparison  with  the  others,  as  shown 
on  Figures  3.3  to  3.7.  Where  the  location  of  a  failure  surface  deviated 
from  the  logical  position  with  respect  to  the  other  failure  surfaces ,  a 
more  reasonable  trial  surface  was  introduced  by  the  author  and  analysed 
as  a  specified  surface.   If  this  surface  was  significantly  more  critical 
(lower  factor  of  safety)  than  the  most  critical  surface  generated  by 
STABL,  the  former  surface  was  used.  These  surfaces  are  drawn  with 
dashed  lines  on  Figures  3.3  through  3.7.  The  reason  why  this  can  happen 
is  explained  in  the  word  "random".   STABL  generates  failure  surfaces 
randomly  within  the  defined  limits,  and  selects  the  10  most  critical  of 
the  surfaces  generated.   However,  it  is  possible  that  there  are  some 
more  critical  failure  surfaces  than  those  found  by  STABL. 
3.3.1.2  Results  of  Stability  Analysis  of  Homogeneous  Slopes 

The  results  of  the  stability  analysis  in  terms  of  factors  of  safety, 
F,  are  given  in  Table  3.1.  Following  the  procedure  outlined  in  Section 
3.2,  critical  failure  surfaces  (Figures  3.3  to  3.7)  were  used  as  speci- 
fied surfaces ,  and  the  factor  of  safety  calculated  for  small  changes  in 
c  '  and  <f '  • 

Another  approach  to  determine  AF(Ax  )*  is  to  plot  F  vs.  x^  keeping 
other  parameters  constant.  This  method  has  the  advantage  that  it  also 
reflects  changes  in  F  due  to  a  slight  change  in  the  critical  failure 
surface  when  the  soil  parameters  change.   It  was  possible  to  apply  the 
method  in  the  case  of  homogeneous  slopes,  because  a  relatively  wide 


The  change  in  the  factor  of  safety  due  to  a  change  in  x  ,  keeping 
other  parameters  constant. 


3^4 


J*-°- 


# 


0.2 


yfc'-°-3 


FIGURE   3.3  CRITICAL  FAILURE    SURFACES  FOR 

VERTICAL    SLOPES   IN   HOMOGENEOUS 
SOIL 
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7%  -0.1 


7%-02 


7^  =  0.3 


FIGURE  3.4         CRITICAL   FAILURE    SURFACES  FOR  SLOPES 
1/2    IN  HOMOGENEOUS   SOIL 
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j£-o. 


3^-0.2 


^-0.3 


FIGURE  3.5  CRITICAL   FAILURE  SURFACES  FOR  SLOPES 

l/l    IN  HOMOGENEOUS   SOIL 
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^K-o.3 
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FIGURE  3.6  CRITICAL   FAILURE  SURFACES  FOR    SLOPES 

2/1     IN    HOMOGENEOUS  SOIL 
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0  =20'- 50° 


rh  -o.i 


7R-0.2 


FIGURE  3.7  CRITICAL   FAILURE    SURFACES  FOR   SLOPES 

3/1    IN    HOMOGENEOUS   SOIL 
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rfr-o.3 


4>'-o° 


XhE0' 


NolS  Dashed  circle  was  forced  deep,  but  result  is  the 

same  as  the  shallower  circle. 

For   yfi  =0.2  and  0.3   the  failure  surface  is  at  limit 
of  boundaries.  (Theoretical  surfaces  extend  to 
infinite  depth.) 


FIGURE   3.7 


(CONT.) 
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Table   3.1 
Factors   of  Safety  for  Homogeneous   Slopes 


1 

c 
Yh 

Nscot  6 

0 

0.5 

1.0 

2.0 

3.0 

0.1 

0 
10 
20 
30 
40 
50 

.677 

.792 

.915 

1.057 

.787 

.987 

1.196 

1.1+31 

1.719 

.936 
1.210 
1.531+ 
1.883 
2.331 

1.152 
1.618 
2.181+ 
2.808 
3.61+2 

0.757 
1.388 
2.087 
2.871 
3.806 
5.011 

0.2 

0 
10 
20 
30 
40 
50 

.801 
.967 
1.120 
1.261+ 
1.1+08 
1.598 

1.088 
1.31+7 
1.577 
1.818 
2.087 
2.1+26 

1.106 
1.575 
1.891 
2.230 
2.622 
3.115 

1.129 
1.817 
2.366 
2.928 
3.587 
1+.519 

1.11+3 
2.039 
2.819 
3.620 
1+.582 
5.879 

0.3 

0 
10 
20 
30 
40 
50 

1.202 
1.368 
1.521 
1.683 
1.860 
2.075 

1.633 
1.886 
2.11+1+ 
2.1+06 
2.701 
3.066 

I.658 
2.232 
2.51+1 
2.888 
3.303 
3.901 

1.707 
2.1+37 
3.021+ 
3.631 
1+.327 
5-232 

1.711+ 
2.67I+ 
3.1+90 
1+.351+ 
5.331 
6.627 

Note:      The  underlined  results   are  obtained  from  specified  failure 

surfaces,   shown     with  dashed  lines  on  Figures    3.3  through   3.6, 
see   also  top  of  page   33. 
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range  in  the  independent  variables  was  studied.  The  dimensionless 

c' 

number  —  is  used  as  an  independent  variable,  reducing  the  number  of 

c' 
independent  variables  to  three:   -r*  ,  <p'  and  8.  The  dependent  variable 

is  the  factor  of  safety,  F.   The  location  of  the  more  critical  failure 

surface  may  be  considered  another  dependent  variable,  see  Figure  3.8. 

From  the  graphical  plots,  the  sensitivity  of  F  to  changes  in  the 
independent  variables  can  be  found.  The  graphs  may  be  plotted  as  func- 
tions in  3-dimensional  space,  giving  F  as  a  function  of  two  independent 
variables,  while  other  independent  variables  (one  in  this  case)  are  held 
constant,  see  Figure  3-9.  For  any  practical  purposes,  the  two-dimension- 
al projections  of  the  three  dimensional  plots  are  the  most  useful 
(Figures  3.10  to  3.12). 

The  variables  are  combined  in  different  ways  and  graphical  results 
obtained,  which  are  convenient  to  use  for  a  quick  estimate  of  the  in- 
fluence of  change  in  an  independent  variable  x  on  the  safety  factor  F. 
These  graphs  are  presented  in  Appendix  A,  Figures  A.l  through  A. 29. 


3.3.1.3  Sensitivity  of  the  Factor  of  Safety  to  Input  Variables 

AF 
To  determine  AF(Ax. )  for  a  particular  problem,  ^    was  determined 

over  the  total  range  of  parameters  investigated.  This  was  done  by  deter- 


mining the  tangents  or  secants  to  the  F  vs.  x.  curves  at  particular 

points.   Since  these  curves  are  based  on  empirical  values,  the  values 

AF  AF 

of  v—  are  subjected  to  some  uncertainty.   Examples  of  variations  of  -^ 

i  i 

with  x  are  illustrated  in  Figures  3.13  through  3.17.  Tables  and  more 

J 

figures  are  included  in  Appendix  A.  The  results  are  summarized  in  Table 
3.2,  page  52. 
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A  general  observation  is  that  t—  (x  =  —  ,  <j> '  ,  or  cot  8)  is 


c  *  c  ' 

relatively  insensitive  to  variations  in   —  for  values  of  — —     larger 

Yh  Yh 

AF 
than  0.1.  An  exception  is  the  t — 7—5-  ratio  for  slopes  steeper  than 

a '  AF       At*1 

1*5°  •   For  —  values  between  0  and  0.1,  -rry     and  -r — — ^  increases  with 

c'  AF  c* 

— —  ,  whereas  ; —  is  constant  or  decreases  with  -r-  . 

Yh  A(cl}  Yh 

Yh 

Figures  3.13  to  3.16  show  that  the  factor  of  safety  F,  is  most 

c' 

sensitive  to  variations  in  the  —  parameter  as  well  as  $'  for  relatively 

Yh      AF 
flat  slopes.  For  8  angles  up  to  1+5° »  w^  decreases  slightlv  with 

Wh;  AF 

increasing  8;  for  8  angles  larger  than  U50  ,  the  decrease  in  a/c\  is 

AF  ^ 

larger.  Also  ttt  decreases  with  increasing  value  of  8;  however,  con- 

AF  _ 
trarily  to  a/C\  ,  the  decrease  is  large  for  8  angles  up  to  U50  to  50°, 

^Yh 
and  less  for  8  angles  larger  than  50° . 

Figure  3.17  shows  the  sensitivity  of  F  to  variations  in  cot  8,  8 

being  the  slope  angle.   Cot  8  was  chosen  as  independent  variable  instead 

of  8,  because  the  relationship  between  F  and  cot  8  was  more  simple  than 

c' 
that  between  F  and  8.   For  values  of  —  higher  than  0.1  and  slope 

Yh 

AF 
angles  8  less  than  U50  ,  7 — —5-  is  almost  independent  of  the  value  of 

Acot  p 

c ' 

—  as  well  as  R,  but  increases  with  increasing  ip.  For  steeper  slopes, 

AF  c '       , ,     . .,  „ 

increases  with  —    as  well  as  with  p. 


Acot  8  yh 

To  investigate  actual  changes  in  factors  of  safety  due  to  variation 

in  the  variables,  specific  variances  were  assigned  to  the  variables  as 

follows : 

A<(>'  =  5°   (constant) 

=T~  =  0.20  =   20£ 
c 


5U 


^  =  0.05  =   5% 
n 

^  =  0.05  =   5% 
Assuming  5#  variance  in  height  h  and  base  of  slope  b,  we  find 

Acot  6  -  rrf  -  I  *  -  °-1 1  -  -  °-1  cot  6 

or     ^2t6=     0.10     =     10* 
cot  8 

~-     =     20$     gives      A(^-)c  =   205$ 

%   .    At   ,    „    gives    M^)h  -    <>,  .    5* 

c ' 
The  total  change     A(— — )   due  to   c',   y,   and  h  will  be 
Yh 


<>" 


1  -   0.2  _ 

Li.052    "    . 


°-     =     -   0.271.   ^     .     27?, 


which  is  less  than  the  sum  of  the  three  (20  +  5  +  5  =  30$).  Hence  the 
procedure  will  overestimate  E  AF  a  little.  For  positive  values  of  M— ) , 
it  is  opposite;  however,  all  variations  are  investigated  with  respect 
to  a  decrease  in  F,  because  this  is  the  critical  situation.   In  prac- 
tice, the  coincidence  of  all  variances  being  in  the  direction  of  less 
stability  is  unlikely,  but  in  the  deterministic  approach,  we  design  for 
the  combination  that  gives  the  most  critical  situation. 

The  specific  variances  are  applied  to  the  five  slopes  investigated, 
in  the  range  of  factors  of  safety  from  0.8  to  3.  The  calculations  are 
tabulated  in  Appendix  A,  Tables  A. 5  to  A. 13.   Figures  3.18  through  3.22 
show  the  relative  sensitivity  of  F  to  the  variables. 
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FIGURE  3.18  RELATIVE  SENSITIVITY   SAF     FOR 
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RGURE  3.20         RELATIVE   SENSITIVITY    ZAF    FOR 
SLOPE  l/l 


100% 


7U- 


< 

M 

S  50% 
o 


0% 


Slope  2/1 


ZAF 


=  28  to  41%  ,  av.  33% 


F  =0.8  to  3 


0.4 


*FXi 
FIGURE  3.21  RELATIVE    SENSITIVITY  EaF  FOR 

SLOPE    2/1 


59 


100% 


to 


o 


0% 


Slope  3/1 


-£^E  =30  to  43%  ,  av.  33% 


F  =  0.8  to  3 


FIGURE  3.22  RELATIVE    SENSITIVITY     EAF    FOR 

SLOPE  3/1 


60 


It  is  seen  that  for  slopes  with  6  >  1*5°  and  —  >  0.1,  c*  is  the 
most  important  parameter  vith  respect  to  variation  in  the  factor  of 

safety  F,  contributing  30%  to  60%   of  the  total  change  in  F. 

c  '    /C  ' 
For  flatter  slones  and  not  too  high  values  of  —  ,  (—  <  0.1  cot  8), 

Yh    yh 

<$>'  is  the  more  significant,  contributing  35%>   to  60%   of  the  total  change 
in  F. 

The  total  change  in  F,  if  variances  in  all  parameters  are  assumed 

to  decrease  F,  amounts  to  2h%   to  1*3$  of  F,  which  is  quite  large.   There 

1   AF 
is  no  relationship  between  F  and  „   . 

If  different  values  of  the  variances  are  assumed,  different  results 
may  be  obtained.  The  values  chosen  here  are  thought  to  be  reasonable. 
Very  little  information  was  found  in  the  literature  about  actual  vari- 
ances in  soil  parameters  as  obtained  from  field  and/or  laboratory  tests. 
Singh  (1970)  uses  maximum  variances  of  50$  for  c'  and  10°  for  <f>',  which 
shifts  the  results  in  the  direction  of  c'  being  most  significant  with 
respect  to  variations  in  F  for  a  larger  range  of  slopes . 
3.3.1.1*  Conclusion  and  Recommendations 

The  results  of  this  sensitivity  study  show  that  for  homogeneous 
slopes  the  factor  of  safety  is  most  sensitive  to  variations  in  the 
strength  parameters  c'  and  <J>*.  These  two  parameters  may  contribute  65$ 
to  80$  of  the  total  change  in  F  due  to  variances  in  the  independent 

variables. 

c  ' 

For  steep  slopes  and  —  >  0.1,  variation  in  c'  has  a  larger  influ- 
Yh 

c' 
ence  on  F  than  has  variation  in  <j>*.  For  flatter  slopes  with  —  < 

0.1  cot  8,  variation  in  <f>*  has  the  larger  influence  on  F. 
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These  two  parameters  are  of  course  the  most  difficult  to  determine 

with  high  accuracy.   This  is  already  reflected  in  the  assumption  of 

—-  =  20£,  whereas  ^  and  ~  ■  5*. 
c  Y      h 

The  fact  that  it  is  impractical  to  obtain  precise  values  of  c'  and 
4>'  is  already  accounted  for,  to  some  extent,  in  the  requirement  that 
the  value  of  the  safety  factor  F  should  lie  in  the  range  of  1.2  to  1.5 
for  slope  stability.  These  values,  although  they  are  based  on  experi- 
ence, may  not  always  assure  a  safe  design.   This  depends  on  how  reliable 
the  soil  parameters  are. 

In  some  countries,  difficulties  in  the  determination  of  soil 
strength  parameters  have  led  to  the  introduction  of  "partial  coeffi- 
cients". These  coefficients  are  factors  that  are  applied  to  the  soil 
parameters,  as  well  as  structural  and  load  parameters,  if  applicable. 
The  factored  variables  are  used  in  the  stability  calculations  and  the 
safe  design  requires  P  i  1.   This  method  has  the  advantage  over  the 
traditional  method  of  stability  analysis,  in  that  the  factors  can  vary 
from  one  variable  to  another,  and  hence  account  for  the  fact  that  some 
variables  are  more  difficult  to  determine  than  others ,  and  that  varia- 
tion in  some  parameters  has  more  effect  on  the  result  than  variation  in 
others.  Table  3.3  shows  some  of  the  partial  coefficients  used  in 
Denmark  [from  Lundgren  and  Brinch  Hansen  (1965)].  For  stability  and 
earth  pressure  the  partial  coefficients  for  c'  (f  )  and  <j>'  (f.  )  are  1.5 
and  1.2,  respectively,  indicating  that  c'  is  the  one  of  the  two  parame- 
ters believed  to  have  the  most  uncertainty  and  the  highest  influence  on 

the  result.   The  value  of  c'  used  in  stability  analysis  would  be  c '  = 

c' 

lab 

f 
c 
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Table  3.3 

Typical  Values  of  Partial  Coefficients  used  in  Denmark  for 
Geotechnical  Designs* 


Symbol 

Partial  Coefficient  for 

formal 

Extra- 
ordinary 

f 
g 

Dead  load,  water  pressure 

1.0 

1.0 

f 
P 

Live  loads 

1.5 

1.5 

f 

V 

Wind  loads 

1.5 

1.0 

r* 

Friction  (stability  & 

earth  pressures) 

1.2 

1.1 

Friction  (bearing  capacity, 
foundations  &  piles) 

1.25 

1.15 

f 
c 

Cohesion  (stability  & 

earth  pressures) 

1.5 

l.i» 

Cohesion  (bearing  capacity, 
foundations) 

1.75 

1.6 

Cohesion  (bearing  capacity, 
piles) 

2.0 

1.8 

f 

a 

Skin  friction,  piles 

2.0 

< 

1.8 

*  From  Lundgren  and  Brinch  Hansen  (1965) 
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For  practical  purposes  it  is  recommended  that  Figures  A. 11  through 
A. 15  he  used  for  determination  of  AF  due  to  variances  in  soil  data. 
For  a  particular  problem,  the  variances  Ac'  and  A<K  should  be  estimated 
from  the  different  possibilities  of  interpreting  the  test  data  (differ- 
ent ways  of  fitting  the  Mohr-Coulomb  envelope  to  the  test  data).   The 
same  curves  may  also  be  used  to  determine  the  maximum  variation  in  the 
c'  and  <J> '  parameters  that  can  be  tolerated  if  a  certain  minimum  factor 
of  safety  is  required  (Step  h   in  Figure  3.1). 

Figures  A. 25  through  A. 29  may  be  used  in  a  similar  way  to  determine 
the  change  in  F  due  to  changes  in  the  geometry  parameters  h  and  3,  when 
the  soil  parameters  are  considered  fixed. 

3.3.2  Sensitivity  Analysis  of  an  Embankment  Founded  on  a  Soil  with  a 
Weak  Seam 
Moving  from  the  simple  problem  of  a  slope  in  homogeneous  soil  to 
problems  encountering  inhomogeneous  soils ,  the  number  of  independent 
variables  increases  significantly.  The  following  items  must  be  deter- 
mined. 

(a)  Independent  variables  x. . 

(b)  Range  of  each  of  these  variables.   The  range  of  some 
variables  may  influence  the  range  of  others. 

(c)  Range  of  variation  in  each  of  the  variables,  Ax  . 
(Constant,  percentage  of  mean,  etc.) 

The  16  variables ,  which  apply  to  an  embankment  built  on  a  founda- 
tion with  a  weak  soil  layer,  are  shown  in  Figure  3.23.  For  some  of  the 
variables,  e.g.,  the  soil  parameters,  certain  limits  in  the  values  can 
easily  be  established.  The  range  in  geometric  variables  may  impose 


6k 


<>;<:,<*»* 


Variables 

Height  and  width  of  embankment 

Depth  to  weak  layer  from  base  of  embankment , 
measured  at  center  line  of  embankment 

Inclination  of  ground  surface 

Inclination  of  weak  soil  layer 

Inclination  of  left  and  right  sloDe 
of  embankment 

Soil  parameters  of  embankment  fill 

Soil  parameters  of  foundation  soil 

Soil  parameters  of  weak  soil  layer 


"ymbols 
h,  B 

D 

i 

a 


6L>eR 


(Y,c,(f>)] 
(Y,c,*)] 
(Y,c,(f>) 


FIGURE  3.23         VARIABLES   FOR  AN  EMBANKMENT  BUILT 
ON  A  FOUNDATION  SOIL  WITH  A  WEAK 
SOIL  LAYER 
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further  limits  on  soil  strength  parameters,  when  we  are  only  interested 
in  problems  having  factors  of  safety  less  than  3.  Many  of  the  limits 
in  the  parameters  must  be  determined  by  trial.  For  example,  at  what 
depth  will  a  weak  seam  in  a  foundation  soil  have  no  effect  on  the  sta- 
bility of  an  embankment? 

Because  of  the  excessive  number  of  possible  combinations  of  the 
variables  (if  each  of  the  16  variables  were  given  3  values,  the  number 
of  combinations  would  be  3   =  l+3,0U6,72l)  ,  it  was  necessary  to  cover 
only  changes  in  those  variables  which  are  considered  to  be  most  impor- 
tant. 

The  variables  that  are  held  constant  are  mainly  those  with  respect 
to  geometry,  particularly  those  which  are  easy  to  get  information  about 
in  a  real  situation,  such  as  the  surface  profile  (i)  and  the  geometry 
of  the  embankment  proper.   The  ground  surface  was  taken  as  horizontal; 
the  embankment  geometry  was  assumed  as:  20  feet  in  height,  ho   feet  of 
crown  width,  and  slopes  1.5/1  (33.7°).   The  number  of  soil  parameters 
is  reduced  by  using  an  average  unit  weight  y  for  &H  soils,  and  taking 
<(>  as  equal  to  zero.   After  this  reduction,  the  following  parameters 
remain  as  variables: 

D        Depth  to  weak  seam 

a        Inclination  of  weak  seam 

y        Average  unit  weight  of  soil 

(c,  <{>)„    Strength  intercept  and  strength  angle 
of  embankment  soil 

(c,  <)>)_,    Strength  intercept  and  strength  angle 
of  foundation  soil 

c        Strength  of  weak  soil  seam 
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The  relationship  between  the  factor  of  safety  F  and  the  unit  weight 

Y  is  relatively  simple,  which  means  that  y  can  he  held  constant,  leaving 

seven  variables . 

In  this  study,  c  has  a  maximum  of  five  values.   Ten  combinations 
w 

of  (c,  <J>)  and  (c,  <J>)p  are  investigated,  for  up  to  nine  different  depths 
D,  and  three  different  angles  a.   See  Table  3.1*. 

3.3.2.1  Factor  of  Safety  for  Failure  Determined  by  a  Weak  Soil  Seam 
The  factor  of  safety  for  an  embankment  founded  on  a  soil  with  a 
weak  layer  should  be  less  than 

(a)  the  factor  of  safety  for  the  embankment  itself 
(toe  failure  with  9  >  0), 

(b)  the  factor  of  safety  for  an  embankment  founded  on 
a  similar  soil  with  no  weak  layer. 

If  the  factor  of  safety  becomes  larger  than  either  (a)  or  (b),  the  weak 
seam  does  not  influence  the  stability  of  the  embankment. 

Where  the  geologic  conditions,  such  as  Joints  in  a  rock  slope  or  a 
weak  soil  seam  beneath  an  earth  slope,  determine  the  shape  and  location 
of  the  critical  failure  surfaces,  the  factor  of  safety  F  can  be 
expressed  as 


F  =  F  +  P.  (3.1) 

c    <p 

where        Fc  =  ^[^tkj  ^  ...  ♦  *n  ej  (3.2) 

and  F.    =     a,    tan  <J>,    +  a0  tan  <f>0  +    •••    +  a     tan  $  (3.3) 

<p  1  1  2  c  n  n 
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Table  3.k 

Range  in  Variables  Studied  for  an  Embankment  on  a  Foundation  with 
a  Weak  Layer 


Variable 


Height  of  embankment 

Width  of  embankment 

Depth  to  weak  layer  from  base  of  embankment 
(measured  at  centerline  of  embankment) 

Inclination  of  ground  surface 

Inclination  of  weak  soil  layer 

Inclination  of  left  and  right  slope  of 
embankment 

Average  unit  weight  of  soil 
Strength  of  embankment  soil 

Strength  of  foundation  soil 


Symbol 

Range 

h 

20  feet 

B 

HO  feet 

D   0  to  60  feet 


a   0  to  11.3 


V*R 


33.7 


Y       130pcf 
(100  to  lUOpcf) 

c  ,  <J>    lOOOpsf,  0° 
2000psf,  0 
500psf,  10C 

lOOOpsf,  0 
2000psf,  0° 
500psf,  10^ 

0   ,  30: 

200psf,  30 
,  32C 


CF'  ^F 


0 


Strength  of  weak  soil 


c    200  to  lOOOpsf 
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The  variables  are 

n  Number  of  soil  types  that  the  failure 

surface  passes  through 

Y  Mean  value  of  unit  weights  of  soils 

above  the  failure  surface 

h  Height  of  slope  or  some  other 

characteristic  length 

c  to  c       Strength  intercepts  of  soil  No.  1  to  soil  No.  n 

(J>1  to  <|>      Strength  angles  of  soil  No.  1  to  soil  No.  n 

a  ,  k.  to  a  ,  k   Geometry  dependent  variables ,  constant  for  a 
certain  slope  geometry  and  soil  profile. 

3.3.2.2  Sensitivity  Analysis 

Using  the  above  expression  for  the  safety  factor  F,  it  may  be  seen 

that  F,  is  independent  of  the  unit  weight  y   as  well  as  h,  the  height  of 

the  embankment.  Any  change  in  y   or  h  will  only  be  reflected  in  F  . 

It  is  easy  to  show  that  the  following  relationships  hold: 

F 
AF        c  /  \  \ 

-r—    =     -  —  ,  for  other  parameters  held  constant       (3.h) 

AF       F 

77-  =  -  : —  ,  for  other  parameters  held  constant       (3.5) 

Ah       h 

From  equations  (3.2)  and  (3.3)  the  sensitivity  of  F  with  respect 
to  the  soil  strength  parameters  c.  and  <f>.  can  be  defined  as 

#    =  ST*.  (3.6) 

Ac .       yh  i 


i 
AF 


=  a. 


A(tan  <f>.  )      i 


(3.7) 
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Sensitivity  of  F  with  respect  to  the  relative  geometry  can  be 
expressed  as 

3F  3F, 

AF  =  E  ~  Ak.  +  I  s-*  Aa, 

3k.    i     3a.    i 
i  1 

=  —■  E  c,  Ak.  +  I  tan  <J>,  Aa.  (3.8) 

Yh    l   l         i   l 

3.3.2.3  Results  of  Stability  Analysis  for  Embankments 

Different  search  routines  in  the  STABL  program  were  applied  to  the 
problem  to  locate  the  shape  and  position  of  the  most  critical  failure 
surfaces.   It  was  found,  as  expected,  that  the  sliding  block  mode  of 
failure  was  the  most  critical,  with  the  base  of  the  block  located  along 
the  weak  soil  seam.   See  Figure  3. 2U  and  Figures  A. 55  through  A. 53. 
However,  this  was  also  the  case  for  homogeneous  soil,  see  Table  A. 1^4 
and  Figures  A. 59  and  A.60.   When  comparing  these  results  with  results 
obtained  in  the  comparative  analyses,  Chapters  V  and  VI,  it  may  be  con- 
cluded that  the  method  of  slices  used  in  the  STABL  program  (Janbu's 
simplified  procedure)  is  likely  to  give  more  conservative  results  for 
irregular  failure  surfaces  than  for  circular  failure  surfaces. 

It  was  also  found  that  a  sliding  block  type  of  failure  surface  with 
active  and  passive  portions  of  the  surface  inclined  at  angles  corre- 
sponding to  the  Rankine  theory,  1+5°  +  <j>'/2  for  the  active  portion  and 
**5°  -  <J>'/2  for  the  passive  portion,  was  more  critical  than  the  sliding 
block  with  randomly  generated  active  and  passive  portions  of  the  failure 
surface  (Figures  A. 57  and  A. 58).   For  this  reason,  the  BL0CK2  sliding 
block  search  was  introduced  in  addition  to  the  BLOCK  search,  see 
Chapter  IV. 


TO 


weak  seam 


FIGURE  3.24      A    SLIDING  BLOCK  TYPE  OF  FAILURE 
SURFACE  IS   MOST  CRITICAL 
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When  the  weak  soil  layer  is  inclined  with  the  horizontal,  the  sta- 
bility should  be  investigated  for  both  sides  of  the  embankment.   This 
was  done  in  a  few  cases,  but  in  most  cases  the  stability  was  investi- 
gated only  to  the  side  of  the  embankment  where  the  weak  layer  is  dipping 
down,  since  this  is  most  difficult  to  change.   On  the  side  where  the 
weak  soil  layer  is  inclined  upwards,  it  may  more  easily  be  replaced 
with  stronger  soil,  if  the  stability  is  threatened. 

The  results  for  the  few  cases  where  stability  was  investigated  to 
both  sides  show,  that  for  foundation  soil  of  the  same  strength  as  the 
embankment  soil,  failure  is  more  critical  towards  the  down-dipping  weak 
seam.   In  the  case  of  stronger  foundation  soil,  it  may  be  more  critical 
towards  the  side  where  the  weak  seam  is  inclined  upwards  ,  see  Figures 
3.25  and  3.26.   In  most  cases  the  foundation  soil  will  be  weaker  than 
the  embankment  soil,  so  that  failure  towards  the  side  where  the  weak 
layer  dips  down  is  the  most  likely. 

The  failure  surfaces  remain  the  same,  or  very  close  to  the  same, 
for  different  foundation  soils  within  a  relatively  large  range  in  soil 
parameters  (y,  c,  <J>),  so  long  as  the  soil  geometry  and  relative  strength 
of  the  weak  soil  remain  unchanged,  see  Figures  A.6l  through  A . 6U .   But 
the  safety  factor  varies. 

As  the  strength  of  the  weak  seam  increases,  the  active  and  passive 
portions  of  the  failure  surface  move  closer  to  each  other.  The  volume 
involved  in  the  failure  decreases,  as  shown  in  Figure  3.27-   Despite 
this  change  in  the  failure  surface,  the  initial  equations  [(3.1),  (3.2), 
(3.3)]  can  still  be  applied  to  some  extent,  because  the  change  in  shape 
and  location  of  the  failure  surface  is  relatively  small.   Exceptions  to 
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FIGURE  3.26  COMPARISON  BETWEEN  FAILURE  TO  LEFT 
AND  FAILURE  TO  RIGHT  SIDE  OF  EMBANK- 
MENT, D=IOFT,   a-H.3?  <*>  =  0 


7^ 


weak  seam 


Type  I  surface 
Type  2  surface 


FIGURE   3.27 


CHANGE  IN  CRITICAL  FAILURE    SURFACE 
WHEN   STRENGTH  OF  WEAK   SOIL  SEAM 
INCREASES 
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this  are  the  situations  where  the  active  part  of  the  failure  surfaces 
intersects  the  ground  surface  beyond  the  crown  of  the  embankment,  see 
Figure  3.27  (surface  type  2).   Tn  these  cases  the  geometry  of  the 
problem  changes  so  much  that  new  values  of  "a"  and  "k"  parameters  apply 
in  equations  (3.2)  and  (3.3).   Type  1  surfaces  (Figure  3.27)  occur  for 
small  relative  depths  of  the  weak  soil  layer  and  small  inclinations  of 
a.   Type  2  surfaces  occur  for  larger  relative  depths  of  the  weak  soil 
layer  and  for  larger  inclinations  of  a. 

From  the  results  (Table  A.ll+),  curves  are  plotted  showing  varia- 
tions of  the  normalized  factor  of  safetv  F  vs.  the  normalized  strength 

n 

Fyh  Cw 

of  the  weak  seam  c  ,  where  F  =  — ' —  and  c  =  —  .   These  curves 
n        n     c£       n    cE 

show  a  very  obvious  change  when  the  failure  surface  changes  from  type  1 

to  type  2,  as  a  break  appears  in  the  otherwise  linear  relationship 

between  F  and  c  ,  see  Figures  3.28  and  3.29.  From  these  curves  the 
n      n 

"k"  and  "a"  values  can  be  found.   Table  3.5  summarizes  the  results. 
Curves  of  "k"  and  "a"  values  are  plotted  vs.  a  for  constant  depth  D 
(Figures  3.30  and  3.31)  and  vs.  depth  D  for  constant  a  (Figures  3.32 
through  3.36),  where  a  is  the  inclination  of  the  weak  soil  seam. 

The  a-graphs  show  that  the  strength  angle  <(>_,  for  the  foundation 
soil  has  a  larger  influence  on  the  stability  than  has  the  strength 
angle  <$>„   of  the  embankment  soil.   The  difference  increases  with 
increasing  depth  D,  as  a_  increases  and  &-   decreases  with  increasing 
depth.   This  seems  reasonable,  as  the  frictional  strength  increases 
with  the  overburden,  and  the  overburden  is  larger  for  the  part  of  the 
failure  surface  that  passes  through  the  foundation  soil. 
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FIGURE  3.28      NORMALIZED  FACTOR  OF  SAFETY  Fn    VS. 
NORMALIZED  STRENGTH  c„   FOR 
EMBANKMENT  20  FEET  HIGH ,  SLOPES 
1.5/1  ,  DEPTH  TO  WEAK  SEAM   5  FEET, 
INCLINED  AT    a=  11.3°   TO  HORIZONTAL 
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FIGURE  3.29      EXAMPLE  OF  Fn  VS.  cn    FOR  EMBANKMENT 
20  FEET    HIGH  ,    SLOPES  1.5/1 ,   DEPTH 
TO  WEAK  SEAM   5  FEET,  INCLINED  AT 
a  -11.3°  WITH  HORIZONTAL 


Table   3.5 
Geometry  Parameters  k  and  a   for  Embankment   Stability  Analysis 
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Surface 
type* 

D 
ft 

tan  a 

kE 

kF 

kw 

aE 

P 

(2) 

0 

0.20 

2.1+1 

1.01+ 

2.80 

1 

0 

0.20 

2.1+1+ 

1.1+7 

1.57 

0.71 

0.86 

1 

5 

0 

2.71+ 

1.06 

1.87 

0.65 

0.72 

1 

5 

0.10 

2.31+ 

I.69 

1.33 

0.83 

1.01 

1 

5 

0.20 

1.90 

2.11 

1.11 

0.1+7 

1.51+ 

(2) 

5 

0.20 

1.50 

1.66 

2.72 

0.51+ 

0.92 

1 

10 

0 

1.99 

1.72 

1.50 

0.82 

1.20 

1 

10 

0.10 

1.69 

2.19 

1.20 

0.59 

1.60 

2 

10 

0.20 

1.15 

2.25 

2.00 

0.1+5 

1.71 

2 

15 

0.20 

1.05 

2.50 

1.1+6 

1 

20 

0 

1.51+ 

2.1+6 

0.88 

0.1+7 

2.25 

2 

20 

0.20 

0.83 

2.82 

1.07 

0.33 

3.01 

2 

30 

0.20 

0.58 

3.28 

0.61+ 

2 

i+o 

0 

O.56 

3. 31+ 

0.5!+ 

0.12 

1+.73 

2 

i+o 

0.20 

0.1+1 

3.60 

0.1+1+ 

0.08 

(5.1+3) 

2 

50 

0.20 

0.27 

3.9*+ 



0.3!+ 

Surface  types  1  and  2  are  illustrated  in  Figure  3.27. 
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FIGURE  3.30      VARIATIONS  OF  k   AND  "a"  COEFFICIENTS 
WITH    a  FOR  CONSTANT  DEPTH  D«0.25h 
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FIGURE  3.31        VARIATIONS  OF  k  AND  "a"  COEFFICIENTS 
WITH   a  FOR  CONSTANT  DEPTH  D-0.5h 
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FIGURE  3.32        VARIATION  OF  k  COEFFICIENTS  WITH 
DEPTH   D   FOR  CONSTANT   a  =0° 
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FIGURE  3.33       VARIATION  OF  k  COEFFICIENTS  WITH 
DEPTH    D   FOR  CONSTANT   0=5.7° 
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FIGURE   3.34       VARIATION  OF  k  COEFFICIENTS  WITH 
DEPTH   D   FOR  CONSTANT   a -11.3° 
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FIGURE    3.35       VARIATION  OF  "a"  COEFFICIENTS  WITH 
DEPTH  D  FOR  CONSTANT  a  *  0° 
(AND  a- 5.7°) 
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DEPTH  D   FOR  CONSTANT   a  =  11.3° 
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With  respect  to  the  strength  intercepts  c,  the  k-graphs  show  that 
for  small  angles  of  a  and  relatively  small  depths  (D  =  0.3h  to  0.5h), 
the  strength  intercept  c„  of  the  embankment  soil  has  a  larger  influence 

on  the  stability  than  the  strength  intercepts  c„  and  c  of  the  founda- 

r       w 

tion  and  the  weak  soil,  respectively.   For  greather  depths,  c^  has  the 

r 

larger  influence.  Comparing  k^  and  k  ,  it  is  seen  that  for  very  small 

angles  of  a   (=  0° )  and  depths  D  =  0.5h,  the  strength  intercept  c  of 

w 

the  weak  seam  has  a  larger  influence  on  the  stability  than  has  the 
strength  intercept  c„  of  the  foundation  soil.  For  larger  depths  and/or 

r 

larger  angles  a,  c,,  has  the  larger  influence. 

r 

Over  the  range  of  variables  investigated,  the  factor  of  safety  F 
decreases  with  increasing  angle  a   or  depth  D  of  the  weak  seam,  to  a 
certain  depth,  then  F  increases  with  increasing  a  or  D,  other  parameters 
being  held  constant.   The  depth  at  which  the  change  occurs  depends  upon 
the  strength  of  the  foundation  soil  relative  to  the  embankment  soil, 
and  to  some  extent  on  the  relative  strength  of  the  weak  soil  seam.  For 
foundation  soil  of  equal  or  lesser  strength  than  the  embankment  soil, 
the  shift  for  F  vs.  a  occurs  at  a  depth  D  approximately  equal  to  twice 

the  embankment  height,  whereas  for  stronger  foundation  soil,  the  shift 

CF 
occurs  at  much  smaller  depths.  Figure  3.37  shows  F  vs.  c  for  —  =  1 

E 
and  2,  and  a  =  0°,  5.7°  and  11.3°.   The  shift  for  F  vs .  D  occurs  at 

depths  D  between  0.25h  and  2h  for  a  foundation  soil  of  the  same  strength 

as  the  embankment  soil.   For  a  stronger  foundation  soil,  F  increases 

with  depth  D  over  the  total  range  of  D.  Figure  3.38  shows  the  variation 

in  F  with  respect  to  a  and  D  for  foundation  soil  of  the  same  strength 

as  the  embankment  soil . 
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FIGURE     3.38        EXAMPLE  OF  VARIATION  OF  THE   FACTOR 
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3 .  3 .  2 .  '*  Results  Applied  to  Sensitivity  Analysis  of  Embankments 
Sensitivity  of  F  with  respect  to  y 

From  Equations  (3.l),  (3.2)  and  (3.3)  the  following  expression 
was  derived 

£  -  -r  <*•*> 

Example:   An  embankment  has  a  factor  of  safety  with  respect  to  strength 
intercept  F  =1.2,  and  a  unit  weight  Y  =  120  pcf.  This  gives 

&    =  .  L£  =  _  o  oi 
Ay     120     U,U1 

If  120  ucf  is  the  average  value  of  several  measurements  of  y, 
and  say  the  standard  deviation  of  these  measurements  is  +  2.5, 
this  gives  a  range  in  F  of 

F  =  1.2  +  AF  =  1.2  +  (-0.0l)-2.5  =  1.2  +  0.025 

or  f-    =   ~  f-    =  lf|o  =  2'1%>   or  F  =  1-115   to  1>225, 


Sensitivity  of  F  with  respect  to  height  h 


£■  -  -r1  (3-5) 

Ah       h 

Example :   Again  we  assume  an  embankment  with  F  =  1.2,  and  a  height  h 
of  20  feet. 

AZ  =  .  1^2  =  _  0>o6 
Ah       20 

Considering  whether  we  can  add  say  5  feet  to  the  height ,  we 
find  AF  =  -  0.06*5  =  -  0.3,  which  gives  F  =  1.2  -  0.3  =  0.9, 
so  the  answer  is  no. 
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Considering  an  uncertainty  in  the  measurement  of  the  height 
during  construction,  which  leads  to  a  height  that  is  21  feet 
instead  of  20  feet.   This  decreases  the  safety  with 
AF  =  -  0.06-1  =  -  0.06,  hence  F  =  l.lU. 
Sensitivity  of  F  with  respect  to  strength  intercept  c 

f~    -     ik-  (3.6) 

Ac .     Yh  l 

l 

Example:   An  embankment  20  feet  in  height  with  slopes  1.5/1  and  a  ^0 

feet  wide  crest  is  built  on  a  foundation  soil  which  has  a 

weak  seam  located  10  feet  below  the  original  horizontal 

ground  surface  at  the  center  of  the  embankment ,  and  inclined 

5.7°  with  the  horizontal.   This  problem  has  a  maximum  of  3 

types  of  soil,  and  the  k  parameters  found  were  k„  =  1.69, 

k»  »  2.19,  k  =1.20.   If  the  variation  in  strength  intercepts 
■jr         v 

is  Ac  =  +  200  psf ,  and  Y  =  125  pcf 

AF  .  AF_  +  AF_  +  AF.  =  J.         +    , 

Ac     Ac_   Ac„   Ac      yh   E    T    w 

E     F     w 


5.08 
125'20 


=   0.002 


and  AF  =  O.Hl,  which  is  quite  high.   If  the  embankment  is 
designed  with  a  safety  of  F  =  1.5,  this  variation  in  c  may 
result  in  an  F  value  of  only  1.1. 

Notice  that  k^  is  almost  twice  as  large  as  kw>  which  means 
that  greater  effort  should  be  used  to  determine  cp  than  cw 
in  this  particular  case.   However,  complete  neglect  of  a 
weak  seam  could  be  catastrophic. 
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Sensitivity  of  F  with  respect   to  strength  angle  <j> 

AF 
AltalT^T     =     ai  (3-T) 

Example :   For  the  same  embankment  as  described  above,  it  was  found  that 
aE  =  0.59  and  a^  =  1.60. 
If  A*.  =  5°,  we  find  for  <J>.  ~   30°  to  35°, 

A(tan  <f>  )  =  0.123  +  AF  =  0.27 

and  for  <j».  -   10°  to  15°, 

A(tan  (J).)  =  0.092  ■*■     AF  =  0.20, 

which  is  also  a  significant  variation  in  the  factor  of  safety, 
although  the  strength  intercept  c  in  this  particular  case  has 
a  larger  influence  on  F  than  has  the  strength  angle  <J>. 

From  Tables  3.5  and  3.6  it  can  be  seen  that  although  k„,  k„  and  k 

h        r  w 

vary  with  location  of  the  weak  soil  layer  (as  represented  by  a  and  D), 
the  sums  of  the  k  parameters  vary  only  a  little.   Hence  the  result  for 
AF(Ac)  is  almost  independent  of  the  problem  geometry,  and  depends  main- 
ly on  the  values  of  Ac.   This  is  not  the  case  with  respect  to  AF(A<}>)  as 
a-j,  increases  almost  linearly  with  depth  D  to  the  weak  seam,  and  a  de- 
creases only  slightly  with  D.   So  for  greater  depths  D  of  the  weak  seam, 
AF(A(J>)  may  be  larger  than  AF(Ac). 
Sensitivity  of  F  with  respect  to  relative  geometry 


AF  =  -r  E  c.  Ak.  +  I   tan  <J>.  Aa.  (3.8) 

Yh    1   1         1   1 
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Table  3.6 
statistical  Analysis  of  Geometry  Parameters  k  and  a 


Dominating 
Failure 
Surface 

D 
ft 

tan  a 

kE  ♦   kp 

*?  +  ku 

kE   +   kF  +  \l 

^  +  *F 

(2) 

1 

0 
0 

0.20 
0.20 

3.H5 
3.91 

3.81+ 
3.01+ 

6.25 
5.1t8 

1.57 

1 

1 

1 

(2) 

5 
5 
5 
5 

0 
0.10 
0.20 
0.20 

3.80 
1+.03 
It.  01 
3.16 

2.93 
3.02 
3.22 
U.38 

5.67 
5.36 
5.12 
5.88 

1.37 
1.81+ 
2.01 

1.1+6 

1 
1 
2 

10 
10 
10 

0 
0.10 
0.20 

3.71 
3.88 
3. '40 

3.22 
3.39 
U. 25 

5.21 
5.08 
5.1+0 

2.02 
2.19 
2.16 

2 

15 

0.20 

3.51* 

3.96 

5.00 

- 

1 
2 

20 
20 

0 
0.20 

i+.oo 
3.65 

3. 31+ 
3.89 

1+.88 
U  .72 

2.72 
3.31+ 

2 

30 

0.20 

3.86 

3.92 

1+.50 

- 

2 
2 

1*0 
1+0 

0 
0.20 

3.90 
1+.01 

3.88 
It.Olt 

1+.1+1+ 
It.  1+5 

1+.85 
(5.51) 

2 

50 

0.20 

U.21 

it. 28 

U.55 

- 

Mean 

3.78 

3.66 

5.12 

2.32 

Variance 

0.08 

0.21+ 

0.29 

1.03 

Standard 

Dev. 

0.28 

0.1+9 

0.53 

1.01 

Variance 
Coef. 

iM 

13.3$ 

10.lt* 

1+3.75? 
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Example:   Consider  the  same  embankment  as  in  the  previous  example, 

with  a  safety  factor  F  =  F  +  F,  corresponding  to  c^  =  1000 

c    (J)  E 

psf,  c„  =  500  psf,  $_  =  10°,  and  c  =  200  psf,  which  gives 

r  r  W 

F  =   J:    (1.69-1000  +  2.19-500  +  1.2-200)  +1.6  tan  10° 

=  1.U9. 
Let  us  assume  that  there  is  an  uncertainty  in  the  determina- 
tion of  the  position  of  the  weak  soil  seam,  reflected  as 
+  2  feet  in  the  depth  D  and  +0.05  in  the  inclination  (tan  a). 
We  now  have  to  determine  which  sign,  plus  (+)  or  minus  (-), 
will  decrease  the  factor  of  safety.   Since  cv   is  larger  than 

c„  and  c  ,  the  value  of  k„  is  compared  to  k_  and  k  .   We  are 

r        w  Jii  r        W 

first  investigating  the  influence  of  a  change  in  the  depth  D. 

Figure  3.33  shows  that  for  a  =  5-7°,  k„  decreases  with  in- 

r, 

creasing  depth  D,  whereas  k„  increases  approximately  the  same 
amount.   Since  c„  =  2c_,  the  result  is  a  decrease  in  F. 

L        r 

Further,  k  also  decreases  with  increasing  depth.   On  the 
w 

other  hand,  a„  increases  very  much  with  increasing  depth  (see 

Figure  3.35),  so  it  is  not  obvious  whether  F  =  F  +  F  will 

decrease  of  increase  with  increasing  depth.  Therefore  both 

situations  have  to  be  investigated. 

AD  =  +  2  feet 

AkE  =  -  0.25,  Akp  =  0.20,  Akw  =  -  0.05,  Aa^,  =  0.2U 

AF  =  lgc. go  (-  0.25-1000  +  0.20-500  -  0.05-200) 

+  0.2^  tan  10°  =  -  0.06  +  0.0H  =  -  0.02, 
which  shows  that  increased  depth  decreases  the  stability, 
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although  the  change  is  negligible.   For  AD  =  -  2  feet,  the 
changes  in  the  "a"  and  "k"  parameters  are  of  opposite  sign, 
but  of  the  same  magnitude,  so  AF  =  +  0.02. 
Applying  the  same  procedure  to  a  change  in  the  angle  a, 
Figure  3.31  shows  that  k„  decreases  for  increasing  a;  k^,  in- 
creases, but  is  almost  constant  for  a   between  5.7°  and  11.3° 

at  this  -particular  depth  (D  =  10  feet  =  0.5h),  whereas  k 

w 

increases  for  increasing  as  well  as  decreasing  a,  because  the 

failure  surface  for  increasing  a  changes  from  type  1  to  type  2. 

However,  a„  increases  with  increasing  a,  so  like  the  previous 

case,  it  is  not  obvious  whether  a  decrease  of  an  increase  in 

a  will  result  in  a  decrease  in  F. 

Atan  a  =  +  0.05 

Ak  =  -  0.25,  Akj,  =  0.08,  Akw  =  0.16,  Aa^  =  0.09 

AF  =  jpc.gQ  (-0.25-1000  +  0.08-500  +  0.16-200) 

+0.09  tan  10°  =  -  0.07  +  0.02  =  -  0.05 
Atan  q  =  -  0.05 
Ak£  =  0.19,  Akp  =  -  0.19,  Akw  =  0.08,  Aa^,  =  -  0.16 

AF  =  ,2c.g0  (0.19-1000  -  0.19-500  +  0.03'200) 

-  0.16  tan  10°  =  0.01+  _  0.03  =  0.01, 
so  F  decreases  for  increasing  q. 
Comparing  the  changes  in  the  factor  of  safety  F  due  to  variations 
in  the  independent  variables,  we  find  that  the  only  significant  change 
in  F  is  due  to  variation  in  the  soil  strength  parameters  c  and  (J). 
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The  above  results  are  based  upon  a  rather  limited  study,  particu- 
larly with  respect  to  the  range  in  the  a  angles  investigated.   For  very 
steep  angles  of  a,  the  influence  of  the  weak  seam  may  be  negligible. 
However,  this  was  not  investigated. 
3.3.2.5  Summary  and  Conclusions 

The  task  of  performing  a  sensitivity  analysis  on  a  slope  stability 
problem  becomes  almost  impossible  when  moving  from  the  simple  case  of 
slopes  in  homogeneous  soils  to  more  practical,  but  also  more  complica- 
ted problems,  involving  different  types  of  soils  and  changing  soil 
profiles.   The  number  of  variables  increases  drastically,  and  the  pos- 
sible combinations  of  these  variables  increases  even  more.   In  this 
study,  the  number  of  independent  variables  was  reduced  from  sixteen  to 
seven,  compared  to  the  three  independent  variables  that  apply  to  the 
homogeneous  slone  stability  problem.   The  results,  in  terms  of  safety 
factors,  occupy  ten  tables,  compared  to  one  for  homogeneous  slopes,  and 
even  so  they  cover  only  a  small  range  in  the  possible  combinations  of 
the  independent  variables. 

Despite  the  limited  investigation,  the  analysis  did  show  that 
almost  the  same  conclusion  may  be  reached  with  respect  to  problems  en- 
countering inhomogeneous  soils,  as  was  reached  for  homogeneous  slope 
stability  problems ,  namely  that  the  dominating  influence  on  the  stabil- 
ity and  changes  therein  are  determined  by  the  soil  strength  parameters 
c  and  (j).   However,  the  neglect  of  a  weak  soil  seam  may  have  a  large  in- 
fluence on  the  result  of  a  stability  analysis,  and  it  may  lead  to 
failure  of  a  slope  that  otherwise  would  have  been  stable  with  a  safety 
factor  of  F  =  1.2  to  1.5.   But  this  is  in  fact  caused  by  using  an  in- 
correct c  parameter  (or  cf>  parameter). 
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IV  MODIFICATIONS  OF  THE  STABL  PROGRAM 

This  Chapter  describes  the  major  modifications  of  the  GTABL 
program  that  were  performed  during  the  course  of  this  research.   Other 
adjustments  and  corrections  are  reported  in  Appendix  D. 

Also  included  in  this  Chapter  is  some  general  advice  to  users, 
specifically  regarding  problems  that  frequently  have  been  raised  by 
users  of  the  STABL  program. 

k.l     THE  MODIFIED  SLIDING  BLOCK  ROUTINE  (BL0CK2) 

During  the  sensitivity  analysis  of  an  embankment  founded  on  a  soil 
with  a  weak  soil  layer,  as  described  in  Section  3-3.2,  it  was  observed 
that  a  sliding  block  type  of  failure  surface  with  active  and  passive 
portions  of  the  sliding  block  surface  inclined  at  (1+5°  +  <j>'/2)  and 
U5°  -  <j>'/2)  with  the  horizontal,  respectively,  gave  the  lowest  factor 
of  safety.   This  type  of  failure  surface  corresponds  to  Rankine' s  earth 
pressure  theory  [Rankine  (l857)]>   See  Figure  k.l. 

As  the  original  STABL  sliding  block  search  (BLOCK)  generates  the 
active  and  passive  portion  of  the  sliding  block  more  randomly,  it  was 
felt  that  a  sliding  block  using  the  Rankine  theory  for  active  and 
passive  portions  of  the  failure  surfaces  would  be  useful  in  the  analy- 
sis of  relatively  simple  slope  stability  problems. 

The  original  sliding  block,  on  the  other  hand,  was  still  felt  to 
be  a  better  solution  for  cases  of  soil  profiles  that  are  complex. 


97 
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FIGURE  4.1  SLIDING   BLOCK  FAILURE  SURFACE 

WITH   RANKINE   ACTIVE  AND  PASSIVE 
ZONES 
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The  new  sliding  block  procedure  (BL0CK2)  is  modified  from  the 
original,  with  much  of  the  original  procedure  retained.  This  applies 
to  the  random  technique  used  to  pick  control  points  on  the  sliding 
block  failure  surface.   Only  the  technique  for  generation  of  active 
and  passive  portions  of  the  sliding  surface  is  changed.   The  data  cards 
are  unchanged  except  for  the  command  code. 

k . 1 . 1  Generation  of  Active  and  Passive  Portions  of  the  Sliding  Block 
According  to  Rankine's  Theory 

The  procedure  for  generation  of  the  active  and  passive  portions  of 
the  new  sliding  block  is  as  follows. 

The  passive  portion  of  the  failure  surface  starts  from  the  point 
of  the  failure  surface  chosen  randomly  within  the  first  box  defined, 
see  Siegel  (1975a,  p.  96).   The  number  of  the  nearest  soil  boundary  lo- 
cated above  this  point  is  determined  by  checking  for  intersection  with 
a  vertical  line  through  the  point.  When  this  soil  boundary  is  deter- 
mined, the  soil  type  below  is  known,  and  the  direction  of  the  first 
line  segment  of  the  passive  wedge  is  determined  as  (U50  -  <f>'  /2 )  with 
respect  to  the  horizontal,  where  4>'  is  the  strength  angle  of  the  soil. 

Each  time  a  line  segment  is  generated,  intersection  with  soil 
boundaries  is  investigated.   If  an  intersection  occurs,  the  previously 
described  procedure  is  repeated  to  determine  the  new  <{>'  angle  to  use. 

The  active  portion  of  the  failure  surface  is  generated  similarly 
and  from  the  point  chosen  within  the  last  box  specified.   The  angle 
from  the  horizontal  is  (1*5°  +  <J>'/2),  instead  of  (U50  -  4>72). 

The  <!>'  angle  used  is  that  defined  as  the  isotropic  strength  angle. 
If  the  soil  is  specified  with  anisotropic  strength  parameters,  the 
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active  and  passive  zones  of  the  sliding  block  BL0CK2  will  still  be 
generated  according  to  the  <J>'  angle  defined  as  the  isotropic  <f>'  angle. 
This  means  that  one  should  specify  the  isotropic  <J>'  parameter  with  the 
value  of  the  anisotropic  <j>'  parameter.   In  the  active  wedge  range 
between  approximately  +  ^5°  and  +  65°,  and  in  the  passive  wedge  range 
between  approximately  -  15°  and  -  1*5°  •   See  Figure  k.2.      If  (J)'  differs 
considerably  from  the  active  to  the  passive  regions,  one  may  have  to 
divide  the  anisotropic  soil  into  two  soil  types,  one  in  the  active 
region,  and  another  in  the  passive  region.   Although  the  failure  sur- 
face will  be  generated  according  to  the  isotropic  <j>'  value,  the  aniso- 
tropic <(>'  (if  different)  will  be  used  in  calculation  of  the  available 
shear  strength  along  the  failure  surface. 

Further,  limitation  boundaries  are  ignored  in  this  failure  surface 
searching  routine.  The  program  prints  out  a  warning  if  they  are  speci- 
fied, but  generates  failure  surfaces  as  if  they  did  not  exist. 

U.1.2  Application  of  the  New  Sliding  Block  Routine 

The  new  sliding  block  routine,  BL0CK2,  is  particularly  suited  for 
relatively  simple  soil  profiles,  where  a  sliding  block  type  of  failure 
is  obviously  the  most  likely.   This  is,  for  example,  the  case  with  the 
problem  analysed  in  Chapter  III,  an  embankment  built  on  a  foundation 
with  a  weak  soil  layer.  If  the  soil  profile  is  complex,  it  is  recom- 
mended that  the  original  routine,  BLOCK,  be  used. 

k.2     OPTION  OF  SEVERAL  PIEZOMETRIC  SURFACES 

The  original  STABL  program  had  three  possibilities  for  defining  the 
pore  water  pressure  within  a  slope.  These  were:  the  groundwater  table; 
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Anisotropic   direction  ranges 
for  soil  types  I  and  2 

♦90° 


Soil  type  2 
/  Define  isotropic 


soil  type  3 
■  weak  seam 


FIGURE   4.2         BLOCK2    SLIDING   BLOCK  ANALYSIS 
APPLIED  TO  ANISOTROPIC  SOIL 
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the  pore  water  pressure  parameter  r  =  — ,  which  may  be  used  to  present 

excess  pore  water  pressures  due  to  shear;  and  finally  the  pore  pressure 

constant  u  ,  that  defines  a  constant  pore  pressure  for  a  defined  zone. 

This  would  seem  to  be  sufficient  to  model  any  type  of  pore  water 
pressures,  equilibrium  and  excess  of  equilibrium,  that  would  be  encoun- 
tered in  a  slope  stability  analysis.   However,  a  modification  which 
allows  different  piezometric  levels  to  be  assigned  to  different  layers 
of  soil,  was  proposed  by  the  Division  of  Materials  and  Tests,  ISHC. 

This  approach  is  primarily  applied  to  end-of-construction  stability 
analysis  of  embankments  built  upon  clayey  foundations,  where  the  con- 
struction of  the  embankment  creates  excess  pore  pressures  in  soil  layers 
with  poor  drainage  characteristics. 

Using  the  r  parameter  in  such  cases  would  give  pore  pressures  that 
were  proportional  to  the  total  overburden  from  the  soil  layer  considered 
to  the  gound  surf ace ,  whereas  the  excess  pore  pressures  imposed  by  the 
embankment  is  believed  to  be  proportional  only  to  the  increase  in  over- 
burden which  the  embankment  represents. 

The  pore  pressure  constant  u  might  be  applied,  but  this  would  re- 
quire an  increase  in  the  number  of  soil  types  to  be  defined  in  the  area 
where  the  embankment  slopes,  and  a  corresponding  increase  in  soil  bound- 
aries.  It  was  felt  that  it  would  be  much  simpler  to  introduce  a  second 
piezometric  surface  to  represent  the  excess  pore  pressures  built  up  in 
the  less  permeable  soil  layers  due  to  construction  of  an  embankment. 

1+.2.1  Program  Modifications 

The  option  of  defining  more  than  one  phreatic  surface  required 
significant  changes  in  the  subroutines  WATER  and  WEIGHT  in  STABL.  WATER 
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is  the  subroutine  that  reads,  checks  and  stores  the  data  for  the 

piezometric  surfaces,  if  any.  WEIGHT  calculates  the  resultant  forces 

acting  on  each  slice.  These  include  the  total  weight  of  the  slice  and 

water  forces  acting  on  the  top  and  bottom  of  the  slice.   To  be  sure 

that  the  correct  unit  weight,  moist  or  submerged,  is  used  for  each  soil 

type  within  a  slice,  it  is  necessary  to  calculate  and  store  the  level 

of  each  piezometric  surface,  where  it  intersects  the  center  line  of  the 

slice.   Only  the  piezometric  surfaces  that  apply  to  the  top  soil  and 

bottom  soil  of  a  slice  are  used  in  the  calculation  of  the  water  forces 

U  and  U  ,  where  U  is  the  water  force  acting  on  the  bottom  of  a  slice, 
a     8        a 

and  U0  is  the  water  force  acting  on  the  top  of  a  slice. 

p 

H.2.2  Change  in  Data  Cards 

To  define  more  than  one  piezometric  level,  an  integer  is  added  to 
the  soil  parameters,  indicating  which  level  applies  to  the  particular 
soil.   If  no  piezometric  surfaces  are  defined,  any  number  may  be  used. 
Also  the  data  cards  for  WATER  are  changed  to  allow  up  to  10  piezometric 
surfaces  to  be  defined.   These  changes  are  all  included  in  the  revised 
STABL  User  Manual. 

U.2.3  Application  of  more  than  one  Piezometric  Level 

The  option  of  defining  more  than  one  piezometric  level  can  be 
applied  to  problems  where  either  artesian  water  pressures  or  a  perched 
water  table  exist,  see  Figure  k.3.      The  program  will  treat  the  perched 
water  table  problem  as  follows.   For  a  slip  surface  located  within  the 
clayey  soil,  the  water  pressure  on  the  bottom  as  well  as  on  the  top  of  a 
slice  will  be  calculated  as  the  hydrostatic  pressure  corresponding  to  the 
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perched  water  table.   If  a  slip  surface  extends  down  into  the  sand  layer, 
the  water  pressure  at  the  bottom  of  a  slice,  if  located  in  the  sand, 
will  be  calculated  as  the  hydrostatic  pressure  corresponding  to  the 
ground  water  table.  The  top  of  the  slice  will  still  have  a  water  pressure 
corresponding  to  the  perched  water  table. 

Another  stability  problem,  where  the  option  of  defining  two  water 
tables  was  applied,  is  that  of  a  reservoir  embankment  analysed  for  rapid 
drawdown,  see  Figure  U.k.      A  slight  modification  of  the  method  proposed 
by  Taylor  (1937)  is  applied.   Taylor  assumes  that  during  the  rapid  draw- 
down no  drainage  takes  place  in  the  embankment.  While  the  water  pressure 
on  the  surface  of  the  slope  is  removed,  the  water  pressure  within  the  em- 
bankment remains  the  same  as  it  was  under  the  submerged  condition.  This 
is  considered  a  very  conservative  assumption  [Dvinoff  and  Harr  (1971)]. 
The  modified  assumption  used  in  the  above  is  that  within  the  clayey  soil 
core  no  change  in  water  pressure  occurs  during  the  rapid  drawdown,  whereas 
in  the  outer  shell  consisting  of  granular  soil  of  very  high  permeability, 
the  water  pressure  changes  according  to  the  drawdown  water  table.   The 
stability  analysis  is  performed  defining  the  drawdown  water  table  for  the 
outer  shell  of  granular  material,  and  the  original  high  reservoir  water 
table  for  the  impermeable  core  material. 

The  strength  of  the  core  material  has  to  be  defined  in  terms  of 
effective  strength  parameters  <t>'  and  c'.   If  the  strength  is  defined 
solely  in  terms  of  a  total  strength  c,  with  <J>  =  0,  the  presence  or  ab- 
sence of  one  or  several  water  tables  does  not  influence  the  result , 
unless  the  water  table  that  applies  to  the  top  soil  extends  above  the 
ground  surface ,  acting  as  a  surcharge . 
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(a)         Embankment  With    Full  Reservoir 


(b)        Water    Table  After  Drawdown 
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FIGURE    4.4  RAPID    DRAWDOWN    STABILITY 

ANALYSIS    FOR  EMBANKMENT 
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This  may  also  be  seen  directly  from  the  expression  for  the  factor 
of  safety  (Section  U.3,  Equation  k.l),   where  the  water  force  at  the 
bottom  of  a  slice,  U  ,  will  disappear  from  the  equation  if  <f>'  =  0, 
whereas  U„,  the  water  force  at  the  top  of  a  slice,  if  any,  still  is 
contained  in  the  expression  for  A  . 

it.  3  MODIFIED  BISHOP  FACTOR  OF  SAFETY  FOR  CIRCULAR  FAILURE  SURFACES 

During  the  comparative  analysis  described  in  Chapter  VI,  it  was 
discovered  that  the  simplified  Janbu  method  of  slices  used  in  the  STABL 
program  does  not  always  provide  conservative  results  compared  to  other 
methods.   In  fact  it  may  give  erroneous  and  unconservative  results 
where  a  failure  surface  intersects  the  ground  surface  at  the  top  of  the 
slope  at  very  steep  angles ,  and  the  strength  of  the  soil  is  defined 
mainly  in  terms  of  a  strength  intercept  c  (or  c').   Surfaces  with  such 
steep  angles  are  seldom  the  most  likely  failure  surfaces,  so  in  most 
situations  we  will  not  meet  this  problem.   But  it  may  become  important 
where  deep  circular  failure  surfaces  are  analysed,  e.g.,  for  an  embank- 
ment on  a  soft  foundation.   This  problem  is  investigated  in  Section  6.k. 

As  a  result  of  the  findings,  it  was  decided  to  add  the  option  of 
using  the  modified  Bishop  factor  of  safety  in  addition  to  the  simplified 
Janbu  factor  of  safety  for  circular  failure  surfaces. 

1+.3.1  Program  Modifications 

The  modification  of  the  program  was  relatively  simple.  The  Janbu 
factor  of  safety  by  the  method  of  slices  is  calculated  by  iteration  of 
the  following  equation 

A   -FA 
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where 


C ' 

A.,  = +  tan  <J>*  (W  (1  -  k  -  k,  tan  a] 

1    cos  a  v    h 


U 
—  +  LL  (cos  B  +  tan  a  sin  B) 


+  Q  (cos  6  +  tan  a  sin  6)} 


A  =  W  {(1  -  k  )  tan  a  +  k,  } 
d.  v  n 


+  U„  (tan  a  cos  f>  -  sin  B) 
+  Q  (tan  a  cos  6  -  sin  6) 

A  =  tan  a   tan  (J)' 

[see  Siegel  (1975a) ] 

This  equation  is  derived  by  taking  the  moment  of  the  forces  acting 
on  the  slices  around  infinity  (overall  horizontal  force  equilibrium). 
If  the  moment  is  taken  around  an  arbitrary  point,  the  equation  becomes 

A,  -PA 

where  y  is  the  vertical  distance  from  the  center  of  the  bottom  of  a 
slice  to  the  moment  center. 

If  we,  for  a  circular  shaped  failure  surface,  take  the  moment 
around  the  center  of  the  circle,  y  =  r  cos  a  (see,  e.g.,  Figure  0.1*0, 
and  Equation  k.2   becomes 


A   -  F  A2 
E  r  cos  a     -  -  -    =  0 
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or  since  r  =  constant, 

A   -  F  A 
Z   cos  a       x?  -   -  £  =  0  (U.3) 

which  is  the  equation  used  for  the  modified  Bishop  factor  of  safety- 
calculation  . 

A  more  detailed  description  and  discussion  of  the  equations  and 
assumptions  applied  in  the  Bishop  method  as  well  as  the  Janbu  method 
may  be  found  in  Chapter  II  of  this  report,  by  Siegel  (1975a),  and 
Wright  (1969). 

U.3.2  CIRCLE2  and  SURBIS 

To  use  the  modified  Bishop  factor  of  safety  option  for  a  circular 
failure  surface  search,  one  should  use  the  command  word  CIRCL2  in  place 
of  CIRCLE.   It  may  also  be  applied  to  a  specified  surface  by  using  the 
command  word  SURBIS  in  place  of  SURFAC.   To  avoid  erroneous  use,  a 
warning  is  printed  out  that  this  method  should  be  applied  only  to 
failure  surfaces  that  approximate  a  circular  shape. 

k.k     EXAMPLE  PROBLEM  FOR  REVISED  STABL  PROGRAM 

The  problem  analysed  in  Chapter  VI,  Section  6.2,  is  used  to  illus- 
trate the  data  setup  and  output  for  the  revised  STABL  program.   Figure 
6.k   shows  the  problem  with  soil  parameters  and  positions  of  piezometric 
surfaces. 

Table  U.l  shows  the  data  setup  for  a  sliding  BL0CK2  analysis.   Note 
that  the  unit  weight  of  water  (the  second  variable  defined  under  the  com- 
mand word  WATER)  is  specified  as  0.,  which  means  that  the  default  value 
of  62. k   pcf  is  used. 
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Table  h.l 
Data  Setup  for  Example  Problem  for  Revised  GTABL  Program 


PRDFIL 

EXAMPLE  PROBLEM 
23  11 

0.  46.  10.  46.  4 
10.  46.  17.  42.  4 
17.  42.  26.  42.  4 
26.  42.  34.  46.  4 
34.  46.  56.  45.  4 
56.  45.  90.  56.  1 
90.  56.  145.  75.  2 
145.  75.  156.  76.  2 
156.  76.  240.  105.  1 
240.  105.  259.  108,  1 
259.  10S.  320.  103.  1 
156.  76.  320.  78.  2 
90.  56.  320.  58.  1 
56.  45.  216.  38.  4 
216.  38.  320.  43.  4 
0.  36.  20.  36.  5 
20.  36.  136.  29.  5 
136.  29.  216.  32.  5 
216.  32.  320.  28.  5 
0.  22.  136.  22.  4 
0.  17.  31.  17.  1 
31.  17.  136.  22.  1 
136.  22.  320.  22.  1 
SDIL 


5 

125.    125.    0.    35.    0.    0. 

1 

125.    125.    0.    33.    0.    0. 

1 

125.    125.    0.    35.    0.    0. 

1 

122.    122.    320.    29.    0. 

0. 

2 

122.    122.    740.    16.    0. 

Q. 

2 

UATER 

2   0. 

5 

0.    45.    56.    45. 

126.    51.    216.    47. 

320.    47. 

-? 

0.    45.    56.    45. 

126.    51.    136.    60. 

150.    64.    250.    112. 

320.    110. 

BLQCK2 

50  2  50. 

50.    29.    90.    26.    8. 

140.    26.    280.    26.    6. 
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Table  k.2   is  a  listing  of  the  corresponding  output.   The  computer 
plot  is  shown  in  Figure  U.5. 

This  problem  illustrates  the  use  of  the  BL0CK2  sliding  block  search, 
and  of  the  WATER  command.  An  example  of  the  CIRCL2  search  can  be  found 
in  Chapter  VI,  Section  6 .h ,  with  data  setup  in  Appendix  C,  Table  C.7. 

*t. 5   SOME  ADVICE  IN  THE  USE  OF  STABL 

This  section  is  intended  to  give  general  advice  in  the  use  of  the 
STABL  program,  and  to  discuss  some  of  the  questions  and  problems  that 
frequently  have  been  raised  by  users  of  the  program. 

First  it  should  be  emphasized  that  a  computer  program  gives  results 
that  are  no  better  than  the  input  data.   But  let  us  assume  that  we  have 
determined  our  problem  parameters  with  reasonable  accuracy,  or  at  least 
with  a  known  variability.   How  do  we  get  the  best  results  from  the  com- 
puter with  the  least  amount  of  effort? 

U . 5 . 1  Preliminary  Analysis  -  Stability  Charts 

Whether  the  problem  is  relatively  simple  or  more  complicated,  it  is 
usually  a  good  time  investment  to  try  to  determine  an  approximate  solu- 
tion by  means  of  stability  charts  for  homogeneous  slopes  as,  for  example, 
presented  by  Taylor  (1937,  19^8)  ;  see  Figures  k.6   and  1+.7. 

A  chart  similar  to  Figure  h.6   was  constructed  based  on  STABL  results 
from  the  analysis  of  homogeneous  slopes  presented  in  Chapter  III.   It  was 
found,  however,  that  charts  constructed  from  equal  contour  lines  of  factors 
of  safety  F  =  2  and  F  =  3  differed  somewhat  from  the  chart  based  on  F  =  1 
contours.   See  Figure  U.8. 

The  STABL  stability  chart  was  compared  to  Taylor's  chart  (Figure 
k.6)  .      Both  charts  are  superposed  on  Figure  h.9.      It  is  seen  that  for 
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Table  k.2 
Output  of  Example  Problem 


—SLOPE  STABILITY  ANALYSIS- 
SIMPLIFIED  JANBU  METHOD  OF  SLICES 
IRREGULAR  FAILURE  SURFACES 


PROBLEM  DESCRIPTION   EXAMPLE  PROBLEM 


BOUNDARY  COORDINATES 


11  TDP   BOUNDARIES 
23  TOTAL  BOUNDARIES 


BOUNDARY 

X-LEFT 

Y-LEFT 

X-RIGHT 

Y-RIGHT 

SOIL  TYPE 

MO. 

<FT> 

<FT) 

<FTJ 

<FT> 

BELOW  BND 

1 

0 

46.00 

10.00 

46.00 

4 

2 

10.00 

46.00 

17.00 

42.00 

4 

3 

17.00 

42.00 

£6.00 

42.00 

4 

4 

£6.00 

42.00 

34.00 

46.00 

4 

5 

34.00 

46.00 

56.00 

45.00 

4 

6 

56.00 

45.00 

90.00 

56.00 

1 

7 

90.00 

56.00 

145.00 

75.00 

2 

8 

145.00 

75.00 

156.00 

76.00 

2 

9 

156.00 

76.00 

£40.00 

105.00 

1 

10 

240.00 

105.00 

259. 00 

108.00 

1 

11 

259.00 

108.00 

320. 00 

108.00 

1 

12 

156.00 

76.  00 

320.00 

78.00 

2 

13 

90.00 

56.00 

320.00 

53.00 

1 

14 

56.00 

45.00 

216.00 

38.00 

4 

15 

216.00 

38.  00 

3£0. 00 

43.00 

4 

16 

0 

36.00 

£0.00 

36.00 

5 

17 

20.00 

36.  00 

136.00 

29.00 

5 

13 

136.00 

£9.00 

216.00 

32.00 

5 

19 

216.00 

32.00 

320.00 

£8.00 

5 

20 

0 

22.00 

136.00 

£2.00 

4 

21 

0 

17.00 

31.00 

17.00 

1 

22 

31.00 

17.00 

136.00 

22.00 

1 

£3 

136.00 

££.00 

320.00 

22.00 

1 

ISOTROPIC  SOIL  PARAMETERS 
5  TYPE<S>  DF  SOIL 


SOIL   TOTAL 
TYPE  UNIT  UT. 
NO.    <PCF> 


125.0 
125.0 
125.0 
1££.0 
1££.0 


SATURATED 

UNIT  WT. 

<PCF) 

125.0 
125.0 
125.0 

182.0 
122.0 


COHESION 

INTERCEPT 

<PSF> 


320 
740 


FRICTION    PORE 
ANGLE    PRESSURE 
<DEG)   PARAMETER 


35.0 
33.0 
35.0 
£9.0 
16.0 


PRESSURE 

CONSTANT 

<PSF> 

0 
0 
0 
0 
0 


PIEZQMETRIC 
SURFACE 
NO. 

1 
1 

1 
£ 
2 
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Table  k.2   (cont.  ) 

£  PIEZOMETRIC  SURFACE<S>  HAVE  BEEN  SPECIFIED 
UNITUEIGHT  OF  WATER  =  62.40 

PIFZOMETRIC  SURFACE  NO.   1  SPECIFIED  BY  5  COORDINATE  POINTS 


POINT 
NO. 

1 
2 
3 
4 
5 


X-UATER 

Y-UATER 

<FT> 

<FT> 

0 

45.00 

56.00 

45.  00 

126.00 

51.  00 

216.00 

47.00 

320.00 

47.00 

PIEZOMETRIC  SURFACE  NO.   2  SPECIFIED  BY  7  COORDINATE  POINTS 


POINT 

X-WATER 

Y-UATER 

NO. 

<FT) 

<FT> 

1 

0 

45.00 

2 

56.00 

45.  00 

3 

126.00 

51.00 

4 

136.00 

60.00 

5 

150.00 

64.00 

6 

250.00 

112.00 

7 

320.00 

110.00 

A  CRITICAL  FAILURE  SURFACE  SEARCHING  METHOD*  USING  A  RANDOM 
TECHNIQUE  FOR  GENERATING  SLIDING  BLOCK  SURFACES?  HAS  BEEN 
SPECIFIED. 


THE  ACTIVE  AND  PASSIVE  PORTIONS  DF  THE  SLIDING  SURFACES 
ARE  GENERATED  ACCORDING  TO  THE  RANKINE  THEORY. 


50  TRIAL  SURFACES  HAVE  BEEN  GENERATED. 
2  BOXES  SPECIFIED  FDR  GENERATION  DF  CENTRAL  BLOCK  BASE 


LENGTH  OF  LINE  SEGMENTS  FOR  ACTIVE  AND  PASSIVE  PORTIONS  OF 
SLIDING  BLOCK  IS  50. 0 


BOX        X-LEFT     Y-LEFT    X-RIGHT    Y-RIGHT      WIDTH 
NO.         <FT>       <FT>       <FT>       <FT>        <FT) 


50.00      29.00      90.00 
140.00      26.00     280.00 


26.00 
26.00 


8.00 
6.00 
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Table  k.2    (cont. ) 


FOLLOWING  ARE  DISPLAYED  THE  TEN  MOST  CRITICAL  OF  THE  TRIAL 
FAILURE  SURFACES  EXAMINED.   THEY  ARE  ORDERED  -  MOST  CRITICAL 
FIRST. 


FAILURE  SURFACE  SPECIFIED  BY  9  COORDINATE  POINTS 


POINT 
NO. 

1 
2 
3 
4 
5 
6 
7 
8 
9 


X-SURF 

Y-SURF 

<FT) 

<FT) 

52.43 

45.16 

73.46 

32.77 

86.96 

£2.61 

£11.84 

26.69 

215.84 

31.99 

219.47 

38.17 

229.39 

57.21 

240.15 

77.03 

256.03 

107.53 

1.658   *** 

FAILURE  SURFACE  SPECIFIED  BY  9  COORDINATE  POINTS 


>OINT 

X-SURF 

Y-SURF 

NO. 

<FT> 

<FT> 

1 

39.19 

45.76 

2 

59.85 

33.60 

3 

69.02 

26.68 

4 

212.99 

25.05 

5 

218.17 

31.92 

6 

221.92 

38.28 

7 

231.78 

57.23 

8 

242.55 

77.06 

9 
( 

258.62 
1.673  *** 

107.94 

FAILURE  SURFACE  SPECIFIED  BY   9  COORDINATE  POINTS 


PDINT 
NO. 

1 
2 

3 

4 
5 
6 
7 
8 
9 


X-SURF 

Y-SURF 

<FT> 

<FT) 

£9.76 

43.83 

45.78 

34.44 

55.09 

27.43 

218.26 

23.08 

224.74 

31.66 

228.83 

38.62 

238.55 

57.29 

249.33 

77.14 

265.39 

103.00 

1.691    *** 
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Table  k . 2  ( cont . ) 


FAILURE  SURFACE  SPECIFIED  BY  9  COORDINATE  POINTS 


PDINT 
NO. 

1 
2 
3 
4 
5 
6 
7 
8 
9 


X-SURF 

Y-SURF 

(FT) 

(FT) 

32.78 

45.39 

51.99 

34.  07 

59.31 

£3.  56 

220. 04 

£3.65 

226.04 

31.61 

230.20 

33.63 

239.89 

57.30 

250.67 

77.15 

266.73 

103.00 

1.696   *** 

FAILURE  SURFACE  SPECIFIED  BY  10  COORDINATE  POINTS 


POINT 
NO. 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


X-SURF 

Y-SURF 

(FT) 

(FT) 

58.03 

45.66 

59.59 

44.34 

80.83 

32.  33 

86.95 

£7.72 

225.52 

£5.34 

£30.13 

31.46 

£34.51 

33.39 

£44.11 

57.34 

£54.90 

77. £1 

£70.93 

103.00 

1.697   *'** 

FAILURE  SURFACE  SFECIFIED  BY  10  COORDINATE  POINTS 


POINT 
NO. 

1 
£ 
3 

4 
5 
6 
7 
8 
9 
10 


X-SUPF 

Y-SURF 

<FT) 

(FT) 

53.60 

45.84 

60.60 

44.80 

81.88 

32.27 

89.76 

£6.3£ 

£30.83 

£6.86 

£34.18 

31.30 

£38.77 

39.  09 

£48.29 

57.38 

£59. 08 

77.26 

£75.09 

103. CO 

1.701    *** 
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Table  k.2   (cont.) 


FAILURE  SURFACE  SPECIFIED  BY  9  COORDINATE  POINTS 


POINT 

X-SURF 

Y-SURF 

NO. 

<FT) 

<FT> 

1 

39.15 

45.77 

£ 

59.81 

33.60 

3 

67.11 

£8.09 

4 

££8.59 

£8.53 

5 

£30.78 

31.43 

6 

£35.19 

38.  9£ 

7 

£44.73 

57.  35 

8 

£55.57 

77.31 

9 

* 

£71.59 
1.710  ksk 

103. 00 

FAILURE  SURFACE  SPECIFIED  BY  9  COORDINATE  POINTS 


POINT 

X-SURF 

Y-SURF 

m. 

<FT> 

<FT> 

J 

£8. SO 

43.40 

£ 

43.80 

34.56 

3 

5£.7<* 

£7.83 

4 

£££.£8 

£5.37 

5 

££6.96 

31.58 

6 

£31.17 

38.  73 

7 

£40.84 

57.31 

8 

£51. 6£ 

77.17 

9 

« 

£67.67 
1.714  *** 

108.00 
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Table  k.2   (cont. ) 


FAILURE  SURFACE  SPECIFIED  BY  10  COORDINATE  POINTS 


POINT 
NO. 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


*xx 


FAILURE  SURFACE  SPECIFIED  BY  9  COORDINATE  PDINTS 


PDiriT 
NO. 


X-SURF 

Y-SURF 

<FT) 

<FT> 

56.32 

45.10 

56.56 

44.98 

77.71 

32.52 

84.40 

27.48 

238.81 

27.56 

241.41 

31.02 

246.38 

39.46 

255. ?4 

57.44 

266.55 

77.35 

282.51 

108.00 

1.725  *** 

X-SURF 

Y-SURF 

<ft;> 

<FT> 

35.03 

45.95 

55.57 

33.  85 

59.55 

30.85 

212.46 

24.48 

218.07 

31.92 

221.81 

38.28 

£31.68 

57.23 

242.44 

77.05 

258.51 

107.92 

1.728  *** 

Table  h .2   ( cont .  ) 
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FIGURE  4.6       STABILITY  CHART  FOR  SLOPES  IN 

HOMOGENEOUS    SOIL  (AFTER  TAYLOR,  1948) 
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FIGURE  4.7       STABILITY  CHART  FOR  SLOPES  IN 
HOMOGENEOUS  SOIL,   *  =  0 
(AFTER  TAYLOR,  1937) 
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FIGURE    4.8    STABILITY  CHART  OF  THE  TAYLOR  TYPE 
FOR  HOMOGENEOUS  SLOPES,  BASED  ON 
RESULTS   FROM  STABL  PROGRAM 
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FIGURE   4.9        COMPARISON  BETWEEN  STABL  AND 
TAYLOR  STABILITY   CHART  FOR 
HOMOGENEOUS    SLOPES 
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slopes  with  slope  angles  3  <  50°,  the  GTABL  solution  is  conservative 

compared  to  Taylor's  solution,  particularly  for  low  values  of  * ' , 

r 

whereas  Taylor  is  more  conservative  for  steeper  slopes.   Taylor's  chart 
is  based  on  a  friction  circle  method  with  assumed  sinusoidal  stress 
distribution  along  the  shear  surface.  This  stress  distribution  should 
give  a  more  accurate  result  than  the  conventional  friction  circle  method 
described  in  Chapter  II,  Section  2.2.H. 

These  charts  may  also  be  used  to  estimate  the  variability  in  F  due 
to  variation  in  problem  parameters  (see  also  Chapter  III). 

Where  water  pressures  are  present,  the  stability  coefficients  pre- 
sented by  Bishop  and  Morgenstern  (i960)  may  be  used. 

The  preliminary  estimation  of  the  factor  of  safety  for  a  problem 
by  simple  means,  provides  a  good  base  from  which  to  start  the  ''TART, 
analysis.   It  gives  a  check  on  the  computer  output,  and  helps  avoid 
some  troubles  that  frequently  occur  when  problems  having  very  low  safety 
factors  (<  1)  are  analysed  by  GTABL. 

it. 5. 2  Use  of  STABL 

The  next  step  is  to  set  up  the  data  for  the  GTABL  analysis.   The 
basic  steps  are  described  in  the  GTABL  User  Manual  [Giegel  (1975b)], 
but  some  additional  information  may  be  useful. 
H.5.2.1  How  to  Define  Limitations  for  Gearch  Routines 

The  simple  stability  charts  can  provide  us  with  preliminary  estimates. 
Again  we  look  at  Figure  h.6.      The  full  lines  on  the  chart  represent  cases 
where  the  most  critical  failure  surfaces  pass  through  the  toe  of  the  slope. 
If  our  problem  is  located  here  ,  we  need  only  investigate  such  types  of 
failures,  hence  the  left  as  well  as  the  right  limitation  point  for 
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initiation  of  the  failure  surfaces  should  he  specified  as  the  x- 
coordinate  to  the  toe  of  the  slope. 

The  dashed  line  on  the  chart  divides  the  problem  into  zone  A  and 

zone  R.   Zone  A  represents  those  problems  where  the  critical  failure 

> 

surface  does  not  extend  below  the  toe  of  the  slope:  a  -  0  for  the  total 

failure  surface,  where  a  is  the  inclination  of  the  failure  surface  with 
the  horizontal. 

If  our  problem  lies  in  zone  A,  we  can  initially  define  our  direction 
limitations  accordingly.   Since  the  default  values  used  in  the  program 
are  -  ^5°  for  the  clockwise  direction  limit,  and  B  -  5°  for  the  counter- 
clockwise direction  limit,  the  use  of  these  values  for  problems  located 
in  zone  A  of  Taylor's  stability  chart  does  not  provide  the  most  efficient 
search  limitations.   Instead  we  should  specify,  for  example,  the  clock- 
wise direction  limit  as  0°,  and  the  counterclockwise  as  B  -  5°. 
k . 5 . 2 . 2  Clockwise  and  Counterclockwise  Direction  Limits 

At  this  point  it  may  be  appropriate  to  discuss  some  of  the  confusion 
that  has  arisen  with  respect  to  the  use  of  the  term  "clockwise"  and 
"counterclockwise  direction  limits". 

The  terms  are  used  relative  to  each  other,  not  to-  the  X-axis,  whereas 
the  angles  are  defined  with  respect  to  the  X-axis. 

An  angle  is  defined  positive  counterclockwise  from  the  positive 
X-axis,  and  negative  clockwise  from  the  positive  X-axis,  see  Figure  k. 10(a) 

A  counterclockwise  direction  limit  does  not  have  to  be  a  counterclock- 
wise angle  and  vise  versa.   See  Figure  it. 10  (b)  .  Both  angles  are  clock- 
wise with  respect  to  the  positive  X-axis,  so  they  are  both  negative. 


1P5 


positive  angle 

*   X-axis 

negative  angle 


(a)       Definition  of  positive  ana  negative  angles 


Clockwise 
direction 

limitation 

«-60° 


-  ->  X-axis 


Counter  clockwise  airection  limitation  «-40< 


soft  clay 


(b) 


FIGURE    4.10      ILLUSTRATION  OF  CLOCKWISE  AND 
COUNTER  CLOCKWISE  DIRECTION 
LIMITATIONS   FOR  INITIATION  OF 
FAILURE  SURFACES 
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1+.5.3  Strength  along  Goil  Boundaries 

If  we  do  not  know  whether  one  soil  is  stronger  than  another,  it 
may  be  tempting  to  specify  a  surface  or  a  sliding  block  search  on  the 
boundary  between  two  soil  layers,  assuming  that  the  program  will  ana- 
lyse the  stability  for  both  types  of  soils  and  give  the  most  critical 
factor  of  safety. 

It  must  be  clearly  understood  that  this  is  not  the  case.   If  a 
failure  surface  is  specified  or  generated  exactly  on  the  boundary  be- 
tween two  types  of  soil,  the  strength  parameters  of  the  soil  above  the 
failure  surface  will  be  used  in  the  stability  analyses,  independent  of 
which  soil  is  the  weaker.   One  should  always  specify  a  failure  surface 
or  a  sliding  block  search  totally  above  or  totally  below  a  soil  boundary. 
If  one  does  not  know  which  case  will  be  the  most  critical,  both  situa- 
tions should  be  analysed. 

k . 5 . U  Analysis  Involving  Earthquake  Loadings 

A  slope  may  be  analysed  for  the  influence  of  earthquake  forces  by 
pseudo-static  earthquake  coefficients.   Some  uncertainty  has  existed 
with  respect  to  this  option.   For  example,  why  is  the  vertical  earth- 
quake coefficient  defined  positive  upwards?  What  value  for  the  cavi- 
tation pressure  should  be  used,  and  how  does  it  affect  the  stability? 
k . 5 . k . 1  Effect  of  Vertical  Earthquake  Loading  on  the  Stability  of  a 
Slope 

Let  us  first  consider  the  case  of  dry  soil.  Specifying  a  vertical 
earthquake  coefficient  k  has,  for  dry  soil,  exactly  the  same  effect  as 
changing  the  unit  weight  of  the  soil  mass  from  Y  to  yd  -  ky). 
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For  granular  soil,  where  the  strength  in  solely  dependent  on  the 
stress  (a*  tan  <f>',  c'  =  0),  the  factor  of  safety  is  independent  of  y, 
so  a  vertical  earthquake  loading  alone  will  neither  decrease,  nor  in- 
crease the  safety  of  the  slope. 

c ' 
Where  c'  ?   0,  the  factor  of  safety  F  is  proportional  to  — ,  see 

Figure  h.Q;   so  a  change  in  Y  causes  a  change  in  F.   F  decreases  with 

increasing  Yi  so  the  critical  direction  of  the  vertical  earthquake 

loading  is  downward,  hence  k  should  be  specified  negative. 

If  part  of  the  sliding  surface  extends  below  a  water  table,  the 

result  changes  for  the  cases  where  <ji  4   0.   The  earthquake  loading  is  a 

dynamic  force,  momentarily  applied.   If  therefore  the  sliding  mass  is 

submerged,  the  component  of  earthquake  forces  normal  to  the  base  of  a 

slice  will  be  carried  by  an  increase  or  decrease  in  the  pore  water 

pressure.  This  means  that  a  soil  with  strength  defined  in  terms  of  c' 

and  4>'  ,  will  act  as  a  soil  with  <)>  =  0  with  respect  to  the  earthquake 

forces,  for  the  part  of  the  soil  that  is  under  the  water  table.   This 

again  means  that  a  granular  soil  (c'  =  0)  that  is  submerged  will  react 

exactly  as  a  cohesive  soil  (<t>  =  0)  with  respect  to  the  earthquake  forces, 

so  k  should  be  defined  negative  also  in  this  case, 
v 

At  this  point  it  should  be  mentioned  that  the  way  the  earthquake 
forces  are  transmitted  to  the  pore  water  in  the  case  of  a  partially  or 
totally  submerged  slice  is  approximate.  Figure  l+.ll  illustrates  the 
principle.   Since  the  Janbu  procedure  as  well  as  the  modified  Bishop 
procedure  satisfies  vertical  force  equilibrium  for  individual  slices, 
equation  (b)  in  Figure  U.ll  is  the  most  correct.   However,  it  has  AS^, 
the  change  in  required  shear  strength,  as  an  unknown,  and  cannot  be 
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(a)  Force  equilibrium  normal  to  base  of  slice: 


AU  +  KHsina  +  Kvcosa  s  0 


(b)  Vertical   force  equilibrium 


AUcosa  +  ASrsina  +  Kv-  0 


FIGURE  4. 


TRANSMISSION  OF  EARTHQUAKE 
FORCES  TO  THE  PORE  WATER 
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solved  directly  with  respect  to  the  change  in  the  pore  pressure  force 
All.   Hence  equation  (a),  the  projection  of  forces  normal  to  the  slice 
base,  is  used  to  give  an  approximate  value  for  AU. 
1+. 5.^.2  Cavitation  Pressure 

The  pressure  in  the  pore  water  between  soil  particles  that  are 
submerged  can  be  negative;  but  it  cannot  go  below  the  cavitation 
pressure,  where  water  starts  to  evaporate.   This  pressure  is  found  to 
be  about  -  1  atmosphere  for  granular  soil,  which  is  equivalent  to 
-  2115  psf  [Ceed  and  Lee  (1967)]. 

Cavitation  is  unlikely  to  occur.   Only  a  large  upward  directed 
vertical  earthquake  force  would  be  able  to  reduce  the  pore  pressure 
to  a  cavitation  pressure  below  -  1  atmosphere,  and  the  critical  direc- 
tion of  the  vertical  earthquake  force  was  found  to  be  downward. 

Horizontal  earthquake  forces,  on  the  other  hand,  would  result  in 
a  decrease  in  pore  water  pressure  for  slices  with  positive  a's,  when 
defined  positive  (outward  directed),  which  is  the  critical  direction, 
see  Figure  '+.12.   However,  this  decrease  in  the  pore  water  pressure  is 
small  (except  for  slices  with  a  >  U50)  compared  to  the  component  of  the 
force  parallel  to  the  base  of  the  slice,  which  acts  as  a  driving  force. 
It  is  most  likely  that  this  driving  force  will  decrease  the  stability 
of  the  slope  until  the  point  of  failure,  long  before  the  component 
normal  to  the  base  of  the  slice  will  create  a  pore  water  pressure  that 
reaches  the  cavitation  pressure. 

If,  however,  the  pore  pressure  does  reach  the  cavitation  pressure, 
the  soil  particles  will  have  to  "carry"  the  remaining  unloading.  This 
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Ct(-) 


AU  «-KMsina>0 


a(+) 


Outward  directed  earthquake  force  (  +  ) 


KH  »  Wkh  ,    W-  weight  of  slice 


FIGURE     4.12  CHANGE  IN  PORE  WATER    PRESSURE 

DUE  TO    A  HORIZONTAL  EARTHQUAKE 
FORCE 
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will  result  in  a  decrease  in  the  effective  normal  stresses,  and  a  de- 
crease in  the  normal  stress  dependent  strength  (<b  strength). 
h . 5 . k . 3  Illustrative  Example 

A  slope  kO   feet  high,  inclined  26.6°  with  the  horizontal  (slope  2/1), 
was  analysed  for  the  effect  of  vertical  as  well  as  horizontal  earthquake 
loadings.   Three  different  situations  were  considered. 

In  case  A  the  strength  of  the  soil  is  defined  in  terms  of  strength 
intercept  c  only  ((f)  =  0).   In  case  B  the  strength  of  the  soil  is  defined 
in  terms  of  strength  angle  <J)'  alone  (c'  =  0);  the  soil  is  assumed  dry. 
Case  C  is  similar  to  case  B,  hut  a  water  tahle  is  introduced. 

The  slopes  are  first  analysed  for  normal  conditions  with  no  earth- 
quake forces,  and  then  for  five  different  combinations  of  earthquake 
loadings,  see  Table  U.3.   Figure  H.13  shows  the  three  slopes  with  the 
critical  failure  surfaces.   The  critical  failure  surface  for  each  case 
remained  the  same  or  very  close  to  the  same  for  all  analyses. 

The  results  agree  with  the  previous  considerations.   For  cases  A  and 
C  the  factor  of  safety  increases  with  a  positive  (upward  directed)  verti- 
cal earthquake  loading  and  decreases  for  a  negative  (downward  directed) 
vertical  earthquake  loading,  whereas  no  change  results  for  case  B  (dry- 
granular  soil). 

A  horizontal  outward  directed  (positive  k  )  earthquake  loading  re- 
duces the  stability  drastically  for  all  three  cases ,  but  most  for  case  A 
(the  (J)  =  0  case)  . 

A  combination  of  horizontal  and  vertical  earthquake  forces  does  not 
give  the  same  result  as  a  direct  superposition  of  the  two  forces  individ- 
ually acting,  but  the  tendency  is  the  same.  An  upward  directed  (positive) 
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Table   It.  3 
Stability  Analysis   of  Slopes   Subjected  to  Earthquake  Forces* 


Factors  of  safety 

Case  A 

Case  B 

Case  C 

Earthquake  loadings 

c  =  lOOOpsf 

0=0 

c*  =  0, 
dry  soil 

*'  =   35° 
water  table 

(1)    ky  =  kh  =  0 

1.2U 

1.1+2 

1.21 

(2)   ky=  0.1,  ^  =  0 

1.38 

1.U2 

1.33 

(3)  ky  =  -  0.1,  kh  =  0 

1.13 

1.1*2 

1.11 

(h)    ky  =  o,  kh  =  o.i 

0.86 

1.12 

0.98 

(5)  kv=kh=0.l 

0.92 

1.09 

1.05 

(6)  ky  =  -o.l,  kh  =  0.1 

0.80 

l.llt 

0.92 

A(3)  in  %   of  (1) 

-  9 

0 

-  8 

A(l»)  in  %   of  (1) 

-  31 

-  21 

-  19 

A(5)  in  %   of  (1) 

-  26 

-  23 

-  13 

A(6)  in  %   of  (1) 

1 

-  35 

-  20 

~2k 

The  slopes   are   Uo   feet  high,   inclined  2/1    (B  =  26.6° 
Y  =   125pcf. 


Case  A 


Case  B 


Case  C 
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c'«0 
</>'«  35° 


FIGURE    4.13        CRITICAL   FAILURE  SURFACES  FOR 
SLOPES  ANALYSED  FOR  INFLUENCE 
OF  EARTHQUAKE   LOADINGS 
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vertical  force  increases  the  safety  compared  to  the  case  where  only  a 
horizontal  force  is  acting,  and  a  downward  directed  (negative)  vertical 
force  decreases  the  safety.   There  is  one  exception,  namely  case  B,  the 
slope  in  granular  dry  soil.   In  this  case  a  combination  of  upward  and 
outward  directed  earthquake  loadings  was  a  little  more  critical  (lower 
factor  of  safety)  than  the  outward  directed  (horizontal)  alone,  and  the 
outward  and  downward  directed  earthquake  forces  were  slightly  less 
critical  than  the  outward  directed  alone.   This  does  not  appear  logical, 
and  may  be  due  to  some  of  the  assumptions  made  in  the  analytical  model, 
e.g. ,  interslice  shear  forces  are  ignored. 
k.5.h.k     Dynamic  Earthquake  Analysis 

The  method  of  earthquake  analysis  provided  in  the  STABL  program  in 
terms  of  a  pseudo-static  loading,  is  the  conventional  approach. 
However,  in  recent  years  dynamic  earthquake  analyses  have  been  intro- 
duced, which  more  accurately  similates  the  forces  that  are  involved 
[Newmark  (I965);  Seed  (1966)].   The  stability  in  a  dynamic  earthquake 
analysis  is  assessed  in  terms  of  embankment  (or  slope)  deformations, 
rather  than  a  single  value  of  the  factor  of  safety.   For  more  important 
structures,  where  the  consequence  of  failure  is  great,  this  type  of 
analysis  should  be  used  [Seed  et  al  (1975)].   Dynamic  analysis  cannot 
be  applied  in  the  STABL  program. 

^•5-5  Misleading  Factors  of  Safety 

Most  users  of  the  STABL  program  will  sooner  or  later  experience  an 
output  with  the  comment  that  the  factor  of  safety  for  a  failure  surface 
is  misleading. 
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There  are  two  reasons  why  a  factor  of  safety  calculated  by  STABL 
can  be  misleading. 

A  problem  occurs  when,  for  some  reason,  the  iteration  process  for 
the  factor  of  safety  does  not  converge.   This  is  usually  associated 
with  a  factor  of  safety  which  is  either  very  high  or  very  low  (below 
unity).   However,  in  most  cases  STABL  will  handle  factors  of  safety  as 
high  as  10  and  as  low  as  0.6. 

In  one  case,  where  a  user  had  troubles  with  no  convergence  for  a 
low  factor  of  safety,  the  program  was  modified  to  define  the  initial 
value  of  the  factor  of  safety  for  the  iteration  process  to  1.2  instead 
of  1.5  (arbitrarily  defined  in  the  program),  and  this  change  made  it 
work  [Maunsell  8c  Partners  Pty.  Ltd.,  Australia  (1977)]. 

A  misleading  factor  of  safety  may  also  occur  if  the  denominator  of 
the  expression  for  the  factor  of  safety  (F  +  tan  a  tan  4>'  )  is  zero  or 
negative.  This  may  happen  for  deep  failure  surfaces,  where  the  angle  a 
(see  Figure  U.12)  has  a  relatively  large  negative  value  for  some  of  the 
slices.  The  normal  force  on  the  base  of  a  slice  will  thus  become  very 
large  or  negative.  This  implies  an  unreasonable  state  of  stress,  hence 
the  result  is  not  reliable. 

The  STABL  program  therefore  always  checks  for  the  sign  of 
(F  +  tan  a  tan  <j>'  )  for  each  slice,  and  if  the  expression  turns  negative, 
the  message  of  a  misleading  factor  of  safety  is  printed,  together  with 
the  coordinates  of  the  failure  surface  in  question,  and  the  factor  of 
safety.  This  factor  of  safety  will  often  be  either  very  high  or  very 
low.   Frequently  an  asterisk  will  be  printed,  indicating  that  the  value 
is  either  larger  or  smaller  than  the  computer  is  programmed  to  handle. 
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The  negative  value  of  (F  +  tan  a  tan  <f>')  will  often  occur  for  prob- 
lems having  an  actual  factor  of  safety  below  unity.   Hence  an  output 
with  misleading  factors  of  safety  should  be  a  warning  that  the  slope 
may  be  unsafe.   Whitman  and  Bailey  (1967)  suggests  the  requirement 
[cos  a  +  (tan  <j>'  sin  a/F)  <  0.2]  be  used  as  a  basis  for  disregarding  a 
result. 

In  no  case  should  a  factor  of  safety  that  is  labeled  as  "misleading1 
be  used  as  a  result.   It  could  be  close  to  an  actual  factor  of  safety, 
but  it  can  also  be  very  far  from  the  true  value. 

Even  if  some  of  the  factors  of  safety  are  misleading,  STABL  may 
give  a  correct  result  for  the  rest  of  the  surfaces  examined.   These 
will  be  listed  under  the  ten  most  critical  surfaces.   If  some  of  these 
surfaces  give  a  factor  of  safety  of  '500",  this  indicates  that  they  are 
surfaces  giving  a  "misleading"  factor  of  safety.   See  also  Appendix  D, 
Section  D.3.3.2. 

If  all  results  by  STABL  are  misleading,  another  method  of  analysis 
has  to  be  used. 


137 


V  COMPARATIVE  DESIGN  STUDY 

A  comparison  was  made  among  results  obtained  by  STABL  and  those  ob- 
tained by  other  methods  of  analysis  for  an  example  design  problem  posed 
by  Ladd  (1975).   All  results  in  terms  of  factors  of  safety  are  sum- 
marized in  Tables  5-**  and  5.5  (pp.  167-168). 

5.1  DESCRIPTION  OF  PROBLEM 

The  problem  is  stated  in  Figure  5.1.  A  30  feet  high  embankment 
for  a  four  lane  highway  is  to  be  constructed  over  80  feet  of  medium  to 
soft  varved  clay.   The  embankment  will  extend  between  bridge  abutments, 
and  a  two-dimensional  analysis  is  appropriate. 

The  design  process  includes  a  consolidation-settlement-preload 
analysis,  influenced  by  drains  to  satisfy  requirements  of  maximum 
allowable  settlements  and  available  construction  time.   However  only 
the  slope  stability  analyses  are  considered  here. 

5.2  TYPES  OF  ANALYSIS 

Analyses  are  performed  for  the  following  situations 

(1)  Undrained  shear  for  end  of  construction  stability. 

This  case  assumes  that  no  consolidation  takes  place  during 
construction. 

(2)  Drained  shear  for  long  term  stability. 

This  case  assumes  that  all  excess  pore  pressure  has  disappeared, 
consolidation  is  completed. 
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(3)  Partially  drained  conditions  for  end  of  construction  sta- 
bility.  This  case  assumes  that  some  consolidation  takes 
place  during  construction,  increasing  the  strength  of  the 
soil. 
Undrained  shear  is  often  examined  by  total  stress  analysis  while 
drained  shear  is  examined  by  effective  stress  analysis.   However,  both 
total  and  effective  stresses  may  in  general  be  applied,  regardless  of 
drainage ,  although  the  results  expressed  in  terms  of  factors  of  safety 
will  be  different  for  the  two  types  of  analysis  unless  the  factor  of 
safety  equals  one.   This  point  is  illustrated  in  Figure  5.2: 
Total  undrained  shear  strength  is  defined  by 


s  =  c  +  a.  tan  d>  (5.1) 

u    u    f      u 

For  saturated  clay  d>     =0,sos     =c;s      is  normally   defined  as 
Tu  u         u       u 

su  =  qf  =  l/2(a1  -  a3)f  (5.2) 

where  l/2(a  -  a  )„  is  the  maximum  shear  stress  at  failure.   In  this 
design  problem,  however,  s  is  defined  as 

s  =  q _  cos  *'  (5.3) 

u    f 

where  q.  cos  <J>'  is  the  maximum  shear  stress  on  the  failure  plane  at 

failure. 

Effective  shear  strength  is  defined  by 

s  _,  =  c*  +  a'  tan  <J>*  (5.M 

el  i         f 

where  a'  is  the  effective  normal  stress  on  the  potential  failure  plane, 
and  c'  and  <f>'  are  the  strength  intercept  and  strength  angle,  respect- 
ively, in  terms  of  effective  stress.   For  normally  consolidated  clay 
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c'  =  0,  hence  s  __  =  a'  tan  <f>'  . 
ell    I 

In  the  case  where  F  =  1,  a'  =  a'  ,  and 

sgff  =  C  +  a^f  tan  <t>'  =  Tff  =  qf  cos  (J)1  =  su   (5-5) 

where  a'   is  the  effective  normal  stress  on  the  failure  plane  at 

failure,  and  T-_  is  the  shear  stress  on  the  failure  plane  at  failure. 

In  design  situations,  however,  we  require  F  >  1.   This  does  not 

change  the  value  of  s  ,  but  the  effective  strength  will  change  since 

a'  is  not  the  same  as  a'   (corresponding  to  undrained  failure)  unless 

the  pore  pressure  parameter  A  =  l/2(l  -  sin  (J)').   For  A  >  1/2  (l  - 

sin  <j>'),  s  will  be  less  than  s    .   Circle  1  (Figure  5-2)  represents 

the  effective  stress  conditions  for  a  point  on  the  potential  failure 

plane.  Assuming  that  a  sample  of  the  soil  can  be  brought  to  the  same 

state  of  effective  stress  and  then  sheared  undrained,  the  sample  will 

fail  when  the  effective  stress  conditions  given  by  circle  2  is  reached. 

This  gives  a  value  of  s  =  q„  cos  (J)'.   However,  an  effective  stress 
u    f 

analysis  would  use  the  strength  s  f„  (see  Figure  5.2)  which  is  usually 
larger  than  s  .   The  actual  shear  stress  on  the  potential  failure  plane, 
T  ,  is  independent  of  the  type  of  analysis,  hence  total  stress  analysis 
will  in  most  cases  yield  a  lower  factor  of  safety  than  effective  stress 
analysis.  The  above  considerations  apply  to  the  ideal  case,  where  un- 
drained and  drained  shear  tests  result  in  the  same  value  of  <j>'.   Fre- 
quently this  is  not  the  case.   In  this  analysis,  the  drained  and  undrained 
direct  simple  shear  tests  gave  values  of  <£'  equal  to  21°  and  9  .^°  , 
respectively.   This  explains  the  very  large  difference  in  results  by 
effective  and  total  stress  analysis  (see  Section  5.5.2.2). 
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5.3  METHODS  OF  ANALYSIS 

In  the  report  by  Ladd  (1975)  two  methods  of  analysis  were  applied: 
the  simplified  Bishop  method  with  corrected  field  vane  strengths  (see 
Section  5-^.1) »  and  the  Morgenstern-Price  analysis  using  SHANSEP*  an- 
isotropic strengths  (see  Section  5.^.2). 

5-3.1  The  Simplified  Bishop  Method 

The  simplified  Bishop  method  assumes  a  circular  failure  surface 
and  applies  the  method  of  slices  with  the  simplifying  assumption  that 
side  forces  between  slices  are  zero.  The  method  satisfies  overall 
moment  equilibrium,  but  not  moment  equilibrium  for  individual  slices. 
It  satisfies  vertical  force  equilibrium,  but  not  horizontal. 

In  Ladd's  report  the  analyses  are  performed  by  a  computer  program 
ICES  LEASE**.  A  grid  system  combined  with  a  search  routine  is  used  to 
locate  the  center  of  the  most  critical  circle. 

5.3.2  Morgenstern-Price  Analysis 

The  Morgenstern-Price  (MP)  method  of  slices  analyses  failure  sur- 
faces of  any  shape  composed  of  straight  line  segments .  Various  assump- 
tions can  be  made  with  respect  to  interslice  forces,  and  the  analysis 
satisfies  complete  equilibrium.   Different  failure  surfaces  have  to  be 
specified  until  the  more  critical  surfaces  are  found.   The  line  of 
thrust  and  the  magnitude  of  interslice  forces  are  examined  to  judge  the 
reasonableness  of  the  result  (see  also  Appendix  C.l).  If  either  one  is 
deemed  unreasonable,  a  new  set  of  side  forces  must  be  assumed  until  a 
satisfactory  solution  is  found. 


*  Stress  History  And  Normalized  Soil  Engineering  Properties. 

**  Integrated  Civil  Engineering  System  -  Lease  Electronic  Accounting 
System. 
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Such  a  procedure  is  rather  time  consuming,  hence  is  only  recom- 
mended for  problems  where  the  simplified  Bishop  circular  arc  does  not 
provide  reasonable  results.   This  may  be  the  case  when  dealing  with 
non-homogeneous  or  anisotropic  soil  conditions. 

5.3.3  Method  Used  in  the  STABL  Program 

The  STABL  program  uses  a  modification  of  the  simplified  Bishop 
method  of  slices  in  which  overall  horizontal  force  equilibrium  is 
satisfied  instead  of  overall  moment  equilibrium.*  This  modification 
is  introduced  in  order  to  be  able  to  analyse  irregular  as  well  as  cir- 
cular failure  surfaces. 

In  case  of  a  circular  failure  surface,  STABL  will  in  most  cases 
give  a  factor  of  safety  somewhat  lower  than  the  simplified  Bishop 
method,  provided  the  same  (or  approximately  the  same)  most  critical 
failure  surface  is  found.** 

The  advantage  of  the  STABL  program  is  most  significant  when  deal- 
ing with  cases  that  require  failure  surfaces  of  irregular  shape.   In 
this  example  design  problem,  the  RANDOM  and  BLOCK  search  routines  were 
used  in  connection  with  the  SHANSEP  anisotropic  soil  strength  para- 
meters where  the  Morgenstern-Price  analyses  were  used  in  Ladd's  report. 
BLOCK  generates  sliding  block  type  failure  surfaces.  RANDOM  generates 
surfaces  of  more  irregular  shape.   The  advantage  of  the  STABL  program  is 
that  it  searches  for  the  more  critical  failure  surfaces,  whereas  the 
engineer  has  to  conduct  this  search  himself  when  using  the  MP-analysis. 


*  This  method  is  also  known  as  Janbu's  simplified  procedure  (Wright,  1969) 
**  Exceptions,  see  Chapter  VI,  Problem  3  and  h,   and  Appendix  C.II. 
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From  the  results  obtained  in  this  comparative  study,  it  is  seen  that  the 
STABL  program  gives  results  very  close  to  those  obtained  by  the  MP- 
procedure . 

5.U  SOIL  STRENGTH  PARAMETERS  USED  IN  STABILITY  ANALYSES 

5.^.1  Corrected  Field  Vane  Shear  Strength 

Bjerrum  (1972)  introduced  the  corrected  field  vane  shear  strength 
to  evaluate  the  undrained  stability  of  embankments  constructed  on  soft 
clay  foundations  by  circular  arc  analysis .  The  correction  factor  u  is 
based  upon  case  studies  of  embankment  failures  and  y  was  plotted  versus 
Plasticity  Index  (Pi).   However,  many  other  case  studies  show  a  wide  de- 
viation from  Bjerrum' s  recommended  curve.   Special  problems  arise  when 
the  clay  is  layered  and  anisotropic:   Which  PI  should  be  used?  How  far 
is  the  circular  arc  from  actual  failure  surfaces  that  would  develop  in 
these  types  of  soils? 

Instead  of  using  Bjerrum' s  curve,  the  correction  factor  for  the 
varved  clay  was  evaluated  based  on  case  studies  of  embankment  failures 
on  varved  clays . 

A  further  check  of  the  stability  analysis  performed  by  the  simpli- 
fied Bishop  method  and  corrected  field  vane  shear  strengths  was  made  by 
analysing  the  undrained  situation  with  the  Morgenstern-Price  method  and 
SHANSEP  anisotropic  shear  strengths. 

5.U.2  SHANSEP  Anisotropic  Shear  Strength 

The  SHANSEP  method  is  described  by  Ladd  and  Foott  (197U).  The 
method  provides  an  evaluation  of  the  undrained  (and  partially  drained) 


H*5 


stress-strain-strength  properties  of  clays  and  consists  of  three  basic 
steps : 

(1)  Determine  the  stress  history  of  the  clay  foundation,  i.e.,  varia- 
tion with  depth  of  the  effective  overburden  stress  a'  ,  the  maxi- 

vo 

mum  past  pressure  a'  ,  and  the  overconsolidation  ratio  OCR  =  a1  /a1  . 
vm  vm  vo 

(2)  Perform  K  consolidated-undrained  tests  on  undisturbed  samples  to 

o 

find  the  relationship  between  the  undrained  strength  ratio  s  /a* 

u  vc 

and  OCR  for  representative  modes  of  failures  (e.g.,  as  shown  on 

Figures  5.3  and  5.*+). 

(3)  Compute  the  variation  in  s  with  depth  for  the  different  modes  of 
failure  from  the  stress  history  and  normalized  strength  data  ob- 
tained in  steps  (l)  and  (2). 

The  method  is  applied  to  initial  in  situ  undrained  shear  strengths 
and  is  also  used  to  predict  increase  in  strength  due  to  consolidation. 
The  anisotropy  of  the  varved  clay  is  taken  into  account  by  the  modes  of 
failure  investigated  in  the  testing  program  (Figure  5.3)  and  applied  in 
the  stability  analysis. 

5.5  RESULTS  AND  COMPARISONS 

5.5.1  Preliminary  Stability  Analyses 

Preliminary  analyses  were  performed  on  the  initial  design  cross- 
section  (Figure  5.1)  for  undrained  end-of-construction  stability  and 
drained  long  term  stability  using  the  simplified  Bishop  circular  arc 
analysis  with  corrected  field  vane  strengths.   For  the  undrained  analyses, 
various  embankment  heights  were  investigated  with  2:1  and  k:l   slopes, 
in  order  to  determine  the  maximum  height  that  could  be  constructed 
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FIGURE    5.3  STRESS  VS.  STRAIN  FROM  CKoU 
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OF  FAILURE  (After  Ladd,  1975) 
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rapidly  without  consolidation  (i.e.,  assuming  no  increase  in  strength 
due  to  consolidation). 

For  comparison,  the  CIRCLE  search  routine  in  STABL  was  used.   How- 
ever, only  the  design  cross-section  corresponding  to  a  30  feet  embanhment 
and  2:1  slopes  was  analysed   (Figures  5.1  and  5-'0.   Data  and  results 
are  summarized  in  Figures  5»'+  and  5.5.   Results  show  STABL  to  give  a 
slightly  lower  factor  of  safety  than  the  simplified  Bishop  method  in  the 
undrained  analysis  and  about  5%  lower  values  in  the  drained  effective 
stress  analyses. 

The  Morgenstern-Price  method  with  SHANSEP  anisotropic  strength  para- 
meters was  also  applied  to  the  undrained,  end-of-construction  stage,  and 
showed  very  good  agreement  with  the  result  obtained  by  the  simplified 
Bishop  method  with  corrected  field  vane  strengths.   The  anisotropic 
strengths  are  given  in  Table  5-1  and  correspond  to  the  relationship 
shown  in  Figure  5.6. 

The  Morgenstern-Price  method  was  modeled  by  the  RANDOM  and  BLOCK 
procedures  in  the  STABL  program,  where  the  RANDOM  routine  randomly  gener- 
ates irregular  failure  surfaces  and  analyses  their  stability,  and  the 
BLOCK  routines  randomly,  but  more  restrictedly ,  generate  sliding  block  (or 
wedge)  modes  of  failure  and  analyses  their  stability.   (For  a  detailed 
description  of  the  surface  generation  routines,  see  Siegel,  1975a). 

The  results  of  Morgenstern-Price  (MP)  and  STABL,  RANDOM  and  BLOCK 
are  given  in  Figure  5.7.  The  factors  of  safety  for  the  STABL  procedures 
are  equal  to  or  slightly  lower  than  the  one  obtained  by  the  Morgenstern- 
Price  procedure.  All  failure  surfaces  lie  within  a  band  of  20  feet 
width. 
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FIGURE      5.6        ANISOTROPIC    UNDRAINED   STRENGTH 
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ANALYSES    (After  Ladd  ,  1975  ) 
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Figure  5.8  shows  the  most  critical  failure  surfaces  and  corres- 
pond ing  factors  of  safety  found  by  STABL  for  the  four  different  types 
of  surfaces. 
5.5.1.1  Frequency  Distribution  Curves  for  the  Factor  of  Safety. 

The  STABL  program  normally  prints  out  only  the  factors  of  safety 
for  the  ten  most  critical  surfaces.   In  this  study,  however,  a  modifi- 
cation provided  output  of  the  factors  of  safety  for  all  surfaces  inves- 
tigated.  A  histogram  was  made  with  intervals  of  AF  equal  to  0.1  (F  = 
1.0  includes  values  of  F  from  0.950  to  I.0H9,  etc.)  for  each  of  the  four 
procedures  used  in  the  undrained,  total  stress  analyses,  i.e.,  CIRCLE 
(Case  1A),  RANDOM  (Case  2B) ,  BL0CK1  (Case  3B) ,  and  BL0CK2  (Case  Ub). 
From  the  histograms,  frequency  distribution  curves  were  drawn  as  shown 
in  Figure  5-9.   It  is  seen  that  each  of  the  curves  has  a  clear  peak 
varying  from  F  =  0.9  to  1.3  for  the  different  cases. 

From  a  statistical  point  of  view,  the  peak  value  is  the  value  that 
is  most  likely  to  occur,  whereas  the  variance  gives  an  indication  of 
the  scattering  in  the  results.   It  is  seen  that  Cases  1A  and  UB  have 
less  variance  than  cases  2B  and  3B.  This  may  indicate  that  1A  and  UB 
surfaces  are  more  likely  than  2B  and  3B  surfaces.   For  a  truer  compar- 
ison, only  B  cases  should  be  compared  with  each  other,  as  soil  parameters 
are  changed  from  Case  A  to  Case  B.   However,  the  curves  are  too  close  to 
each  other  to  justify  specific  conclusions  with  respect  to  individual 
behavior. 

The  frequency  distribution  curve  can  also  be  interpreted  in  terms 
of  likelihood  of  failure.  For  the  CIRCLE  case,  9  out  of  100  surfaces 
have  a  factor  of  safety  less  than  1.0,  which  means  that  the  probability 
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FIGURE     5.9  FREQUENCY    DISTRIBUTION  CURVES 

FOR   THE    FACTOR  OF  SAFETY, 
CASES    IA,  2B,3B,4B 
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of  failure  is  9  out  of  a  100  or  9%.      When  dealing  with  probability, 
however,  many  more  things  have  to  be  taken  into  consideration.   The 
soil  parameters  should  not  be  handled  as  exact  values  ,  but  distributed 
around  mean  values  with  a  certain  variance  that  can  be  estimated  from 
laboratory  tests.   Gee  for  example  the  work  by  Athanasiou-Grivas  (1976). 
This  aspect  is  not  considered  here. 

5- 5-2   Final  Design  Stability  Analyses 

From  preliminary  stability  analyses  it  was  found  necessary  to  im- 
prove the  stability  of  the  design  cross-section  by  placing  berms  15 
feet  high  and  approximately  60  feet  wide  on  both  sides  of  the  embank- 
ment.  Settlements  were  then  calculated,  and  it  was  found  that  to  com- 
pensate for  settlements,  the  height  of  the  fill  should  be  increased  by 
8  feet.   To  accelerate  the  rate  of  settlements,  sand  drains  should  be 
installed  under  the  embankment  and  additional  load  applied  by  a  tempo- 
rary surcharge.   It  was  found,  however,  that  for  a  minimum  F  of  1.2, 
the  maximum  safe  height  is  only  k2   feet.   This  height  would  still  lead 
to  excessive  post  construction  settlements  if  the  required  time  schedule 
were  followed.   It  was  suggested  that  more  fill  could  be  placed  if  the 
construction  was  carefully  monitored  with  a  fairly  extensive  program  of 
field  measurements.   For  this  example,  however,  the  stability  analyses 
of  the  final  design  is  based  unon  the  \2   feet  of  fill,  with  stability 
berms  as  mentioned  above.   Two  situations  are  investigated  for  the 
final  design:   the  partially  drained  condition  and  the  long  term, 
totally  drained  condition. 
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5.5-2.1   Partially  Drained  Stability  Analyses 

The  1*2  feet  high  embankment  with  stability  berms  was  investigated 
for  a  partially  drained  condition  corresponding  to  95%   consolidation. 
At  this  stage  almost  8  feet  of  settlement  has  taken  place,  see  Figure 
5.10.   The  stability  analysis  was  carried  out  as  a  total  stress  analysis 
with  Morgenstern-Price  wedges  and  SHMSEP  anisotropic  strength  values 
increased  by  consolidation.   See  plots  of  strengths  in  Figure  5-11  and 
Table  5.2.   STABL  was  applied  with  all  four  types  of  analyses  and  aniso- 
tropic soil  parameters. 

Results  of  MP  and  STABL  analyses  are  given  in  Figure  5.12.   STABL ' s 
four  procedures  (Cases  1E.1  through  i*E.l)  give  results  from  +  6%   to  -  5$ 
of  that  obtained  by  the  MP  analysis. 
5.5.2.2  Long  Term  Stability  Analyses 

For  the  long  term  stability  the  embankment  is  analysed  for  a  final 
design  soil  profile  corresponding  to  0  feet  of  settlement  and  a  final 
height  of  30  feet  above  the  ground,  which  means  that  h    feet  of  material 
is  removed  compared  to  the  partially  drained  stage. 

Conventional  drained  effective  stress  analyses  were  performed  with 
Bishop  and  STABL  circular  arc  failure  surfaces.   Data  and  results  are 
given  in  Figure  5-13.   Results  obtained  by  STABL  are  approximately  13$ 
less  than  results  by  the  Bishop  method. 

Total  stress  analyses  were  also  applied  with  SHAH SEP  anisotropic 
strengths  corresponding  to  stresses  under  the  embankment  at  100%  consol- 
idation, see  Figure  5-1^  and  Table  5-3.   Again  MP  and  STABL 's  four  types 
of  analysis  were  applied.   Figure  5-15  shows  the  results.   The  factors 
of  safety  obtained  with  the  four  procedures  of  STABL  (Cases  IE. 2  through 
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UE.2)  varied  from  +  k%   to  -  J%   of  that  obtained  by  the  MP  analysis. 
Total  stress  analysis  is  not  normally  used  to  investigate  long 
term  stability,  where  it  is  assumed  that  the  embankment  foundation  has 
consolidated  100/*.   A  reason  for  using  total  stress  analysis  in  this 
case,  is  the  large  difference  between  the  factor  of  safety  by  the  total 
stress  analysis  and  the  effective  stress  analysis  (1.31*  vs.  2.25). 
Also,  it  is  believed  that  if  a  failure  occurred,  it  would  occur  rapidly 
and  develop  excess  pore  pressures.   Thus  even  though  the  embankment  was 
almost  100/u  consolidated,  an  undrained  shear  would  occur,  and  the  un- 
drained  total  stress  analysis  is  considered  most  appropriate. 

5.6  SUMMARY  Aim  CONCLUSION 

Results  of  stability  analyses  obtained  by  STARL  for  a  design 
example  problem  (Figure  5.1)  were  compared  with  results  obtained  by 
programs  based  on  the  simpli  fied  Bishop  method  and  Morgenstern-Price 
procedures.   The  results  are  summarized  in  Tables  $.h   and  5- 5-   Most  of 
the  cases  investigated  by  GTABL  are  documented  with  graphical  computer 
plots  in  Appendix  3. 

For  the  same  type  of  failure  surfaces,  GTABL  gave  lower  values  of 
the  factor  of  safety  than  the  other  programs.   The  difference  was  in 
most  cases  less  than  10$,  with  a  maximum  of  l6%    (0.83  vs.  1.05  for  the 
preliminary  design,  undrained  total  stress  analysis).   Considering  the 
uncertainty  in  the  input  data,  this  difference  is  tolerable. 

For  problems  with  complicated  soil  conditions,  GTABL  is  recommended 
compared  to  the  other  programs,  because  STABL  utilizes  the  computer 
more  in  the  search  for  the  more  critical  failure  surfaces.   STABL 
usually  gives  a  factor  of  safety  on  the  conservative  side  (see  Section 
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2.3).   If  for  any  reason  it  is  desirable  to  determine  a  "more  precise" 
value  of  the  factor  of  safety,  other  methods  of  slices  may  be  incorpo- 
rated in  the  CTABL  program  (see  Chapter  VI;  the  modified  (or  simplified] 
Bishop  method  was  later  introduced  as  a  permanent  option  in  the  CTABL 
program,  applicable  to  circular  failure  surfaces,  see  Chapter  IV). 

Differences  between  factors  of  safety  for  the  different  kinds  of 
failure  surfaces  generated  by  CTABL  is  less  significant  in  this  example 
problem  than  in  the  example  with  a  weak  layer  of  soil  investigated  in 
Chapter  III.   The  largest  difference  (12$)  is  between  CIRCLE  and  BL0CK2 
for  the  total  stress  analysis. 
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VI   COMPARISON  OF  STABL  WITH  EXISTING  PROCEDURES  FOR  SLOPE  STABILITY 
ANALYSIS  USED  BY  INDIANA  STATE  HIGHWAY  COMMISSION  (ISHC) 

Four  problems  from  the  Division  of  Materials  and  Tests,  ISHC,  have 
been  analysed  for  stability  by  existing  ISHC  procedures  and  by  STABL. 
The  first  problem  involves  analyses  by 

a)  the  computer  program  NYSTAB*  which  is  based  on  circular  arc 
failure  surfaces  and  the  modified  Bishop  method, 

b)  hand  calculated  sliding  wedge  analysis,  and 

c)  STABL. 

The  second  and  third  problem  involves  analyses  by 

a)  a  computer  program  used  by  ATEC**  in  which  circular  arc 
failure  surfaces  are  analysed  by  the  modified  Bishop 
method,  and 

b)  STABL  and  modifications  of  STABL. 
The  last  problem  is  solved  by 

a)  NYSTAB, 

b)  STABL  and  modifications  of  STA3L,  and 

c)  hand  calculation. 


*  NYSTAB  -  New  York  State,  Department  of  Transportation,  Slope  Stability 

Analysis . 
**  ATEC  =  American  Testing  and  Engineering  Corporation,  Indianapolis, 

Indiana. 
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6.1  PROBLEM  1 

The  problem  is  presented  in  Figure  6.1.  An  embankment  of  granular 
fill  is  resting  on  a  foundation  of  28  feet  of  granular  material  under- 
lain by  clay  of  medium  to  low  strength. 

The  factor  of  safety  for  the  embankment  slope  alone  is  given  by 

_        tan   <t>' 

=   tli — 6      '  where  *     is   the  effective   strength   angle,   and  8  is   the 

angle  of  the  slope.  In  this  case  <J>'  =  35°,  tan  8  =  32/56,  hence 

F  =  1.23.  This  result  was  also  found  by  STA3L,  circular  search  with  no 

restrictions  on  initial  angle,  the  critical  failure  surface  being  a 

shallow  circle. 

Stability  with  respect  to  the  underlying  clay  soil  was  analysed 

with  forced  deep  failure  surfaces,  both  circle  and  sliding  block  shapes. 

The  result  from  NYSTAB  is  F  .   =1.32  for  a  deep  circle.  Hand  calculated 

mm 

sliding  block  analvsis  gave  F  .   =1.26  for  the  deepest  surface  analysed 

mm  J 

by  this  method  (medium  depth),  and  3TABL  gave  F  .   =  l.l6  for  a  deep 
sliding  block.  All  sliding  blocks  have  approximately  horizontal  bases 
and  Rankine  active  and  passive  failure  zones. 

Figures  6.2  and  6.3  show  the  failure  surfaces  investigated  and 
corresponding  factors  of  safety.  All  results  are  summarized  in  Table 
6.1.  3TABL,  sliding  block  search,  gives  a  factor  of  safety  6%   lower 
than  the  one  found  by  hand  calculation  for  the  same  depth  (l.l8  vs. 
1.26),  whereas  the  deeper  sliding  block  found  by  STABL  is  more  critical 
{Q%   lower,  1.16  vs.  1.26). 

6.2  PROBLEM  2 

The  soil  profile  and  strength  data  for  problem  2  are  shown  in 
Figure  6.U.   An  embankment  of  granular  material  is  constructed  on  a 
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foundation  of  somewhat  clayey  material,  introducing  an  excess  pore 
water  pressure  in  this  material  corresponding  to  piezometric  surface 

#  2.  The  free  groundwater  table  is  represented  by  piezometric  surface 

#  1.  The  distance  between  piezometric  surface  #  1  and  #  2  represents 
the  excess  pore  water  pressure  induced  by  the  embankment  in  the  clayey 
foundation. 

The  ATEC  program  with  modified  Bishop  circular  arc  analysis  gave 
for  the  three  circles  (l,  2  and  3)  values  of  F  from  1.86  to  2.20,  see 
Figure  6.5- 

STABL  was  modified  to  handle  more  than  one  water  surface  (see 
Chapter  IV"),  and  factors  of  safety  were  obtained  for  the  three  circles 
from  the  ATEC  analyses  as  well  as  for  STABL  circular  search  and  sliding 
block  searches. 

The  circular  arc  failure  surface  proved  to  be  reasonable  for  this 
problem  (a  sliding  block  search  gave  a  factor  of  safety  only  2%   lower 
than  the  circle  search).  To  extend  the  comparison  between  different 
methods  of  analysis,  the  STABL  program  was  modified  to  calculate  the 
factor  of  safety  by 

1)  the  modified  Bishop  method,  and 

2)  Spencer's  method. 

The  Bishop  method  required  only  a  minor  change  in  the  subroutine 
FACTR  of  the  program,  whereas  the  Spencer  method  required  a  complete 
change  of  this  subroutine.  The  source  for  the  Spencer  method  was  Wright 
(1969).  See  Appendix  C  for  a  listing  of  subroutine  FACTR  by  both  methods, 
Later,  the  modified  Bishop  method  was  introduced  as  a  permanent  option 
in  the  STABL  program,  applicable  to  circular  failure  surfaces  (see 
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Chapter  IV).  The  FACTR  subroutine  used  in  the  present  STABL  program  is 
included  in  the  listing  of  the  program,  Appendix  E. 

The  results  of  the  different  methods  of  analysis  are  summarized 
in  Table  6.2.  STABL  with  the  normal  procedure  of  factor  of  safety  cal- 
culation (Janbu's  simplified  procedure)  gives  the  most  conservative 
values  of  F,  with  a  minimum  for  a  circular  arc  11$  lower  than  the  ATEC 
minimum  value  (1.66  vs.  1.86).  STABL  with  the  modified  Bishop  procedure 

gives  the  same  value  of  F  .   as  the  ATEC  method,  but  for  a  different 

min 

circle  (see  Figure  6.5).  STABL  with  Spencer's  procedure  gives  a  value 

slightly  below  the  ATEC  F  .   value,  also  for  a  different  circle.  The 
*  nun 

maximum  difference  in  results  occurs  between  STABL,  simplified  Janbu 
procedure,  sliding  block  search,  and  the  ATEC  solution  for  a  circular 
arc,  namely  12$  (1.63  vs.  1.86).  It  should  be  noted,  however,  that  the 
result  for  the  sliding  block  analysis  with  Spencer's  procedure  is  16$ 
larger  than  the  result  for  the  normal  STABL  procedure  (1.95  vs.  1.63). 
This  may  indicate  that  the  factor  of  safety  by  STABL  is  too  low  for 
the  sliding  block  analysis.  The  sliding  block  failure  surfaces  are 
shown  in  Figure  6.6. 

For  the  Spencer  result,  circle  5,  the  lines  of  total  and  effective 
thrust  are  shown  in  Figure  6.7-  They  look  reasonable  except  for  the 
right  termination.  Also,  interslice  factors  of  safety  with  respect  to 
shear  have  been  calculated.  The  minimum  value  was  2.12,  compared  to  an 
overall  factor  of  safety  F  =  1.83  for  the  sliding  surface,  indicating 
that  the  result  is  reasonable.  See  Appendix  C  for  calculations  of 
effective  line  of  thrust  and  interslice  factors  of  safety,  and  for  more 
information  about  the  influence  of  the  line  of  thrust  on  the  factor  of 
safety. 
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Table  6.2 
Comparative  Study,  Problem  2 
Factor  of  Safety  for  Different  Methods  of  Analysis 


Surface 

Method  of  Analysis 

1 

ATEC 

mod. Bishop 

Circ.  arc 

2 

STABL 

3 

STABL 

mod. Bishop 

Circ.  arc 

U 

STABL 

Spencer 

Circle  1 
2 
3 
h 
5 

BL0CK1-1 
"   1-2 

BL0CK2-1 
"  2-2 

1.86 
2.13 
2.20 

1.T5 
1.80 
1.T5 
1.66 
1.68 

2.01 
1.63 
2.09 
1.61+ 

1.95 
2.16 
2.01 
1.88 
1.86 

1.95 
2.15 
1.98 
1.85 
1.83 

1.95 
2.01 

Max.  Diff. 

0% 

-12% 

0% 

-  2% 

Circles  1,  2  and  3  from  ATEC  analysis, 
Circle  It  from  STABL,  circular  search,  case  2, 
Circle  5  from  STABL,  circular  search,  cases  3  &  *+. 
BL0CK1-1  and  2-1  from  STABL,  sliding  block  search  in 

soil  layer  no.  U,    case  2, 
BL0CK1-2  and  2-2  from  STABL,  sliding  block  search  in 

soil  layer  no.  5,  case  2. 


Between  minimum  values,  compared  to  method  1  (ATEC-solution) 
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6.3  PROBLEM  3 

Figure  6.8  shows  the  problem.   An  embankment  supported  on  one  side 
by  a  retaining  structure  is  investigated  for  a  total  failure. 

The  same  procedures  as  for  problem  2  were  applied  to  this  problem, 
however,  the  Spencer  method  did  not  converge  for  most  of  the  failure 
surfaces  investigated.   Even  with  convergence,  a  further  check  showed 
the  result  to  be  incorrect.   See  Appendix  C,  Tables  C.ll  and  C.12,  and 
Figure  C.13. 

The  results  for  circular  arc  analysis  by  ATEC  and  STABL  are  shown 
in  Figure  6.9.   For  the  ATEC  failure  surface,  STABL  gives  a  factor  of 
safety  1%   lower  than  the  ATEC  value  (2.15  vs.  2.32),  whereas  a  circular 
search  by  STABL  results  in  a  failure  surface  having  a  STABL  factor  of 
safety  16%   lower  than  the  ATEC  value.   If  the  safety  factor  is  calcu- 
lated by  the  modified  Bishop  method,  the  same  two  surfaces  result  in 
factors  of  safety  1%   to  3%   higher  than  the  ATEC  value.   Hence  the  values 
found  by  the  normal  STABL  procedure  are  somewhat  low.   Even  lower  values 
are  found  by  a  STABL  sliding  block  search,  namely  as  low  as  F  =  1.5k, 
which  is  3k%   lower  than  the  ATEC  minimum  (see  Figure  6.10). 

In  applying  the  Spencer  procedure  to  the  same  failure  surfaces ,  all 
cases  except  one  did  not  converge  in  the  iteration  process  for  the  fac- 
tor of  safety  calculation.   The  one  case  that  did  converge  gave  a  safety 
factor  of  2.30,  where  a  value  of  1.65  was  obtained  by  STABL.   The 
Spencer  result  was  checked  by  calculation  of  interslice  factors  of 
safety  with  respect  to  shear.   These  were  found  to  be  less  than  or 
equal  to  l.k    (see  Appendix  C,  Table  C.12),  hence  the  value  of  2.80  is 
not  correct.   The  results  of  the  analyses  are  summarized  in  Table  6.3. 
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Table  6.3 
Comparative  Study,  Problem  3 
Factor  of  Safety  for  Different  Methods  of  Analysis 


Surface 

Method  of  Analysis 

1 
ATEC 
mod. Bishop 
Circ.  arc 

2 

STABL 

3 

STABL 

mod. Bishop 

Circ.  arc 

k 

STABL 

Spencer 

Circle  1 

2.32 

2.15 

2.38 

- 

2 

1.95 

2.3<+ 

- 

BL0CK2-1 

1.72 

- 

"   _2 

1.63 

- 

"   -3 

1.80 

- 

"   -1+ 

1.65 

2.80** 

"   _2* 

1.5U 

- 

Max.  Diff. 

0% 

-3k% 

1% 

*  Result  for  search  with  50  surfaces  instead  of  only  25. 
**  Calculation  of  interslice  factor  of  safety  gave  Fi  l.k, 
so  2.80  is  not  correct. 
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6.k      PROBLEM  k 

The  problem  with  soil  profile  and  strength  data  is  shown  in  Figure 
6.11.   An  embankment  of  cohesive  material,  9  feet  high,  is  built  on  a 
relatively  soft  foundation.   The  strength  of  the  foundation  decreases 
with  depth  until  a  firm  strata  is  reached  at  approximately  13  feet 
below  the  Ground  surface. 

A  search  with  the  NYSTAB  program  gave  a  factor  of  safety  F  =  1.1*1 
for  a  deep  circular  failure  surface  (circle  1,  Figure  6.12).   This  sur- 
face was  analysed  by  the  STABL  program,  resulting  in  a  factor  of  safety 
F  =  2.21,  51%   higher  than  the  result  by  NYSTAB.  A  circular  search  with 
STABL  gave  F  =  1.7U,  while  specifying  the  critical  circle  according  to 
Taylor  (from  Figure  U.7,  page  120),  for  a  similar  slope  geometry,  but 
homogeneous  soil  with  <J>  =  0  (circle  2,  Figure  6.12),  resulted  in  a  STABL 
factor  of  safety  F  =  1.70.   This  value  is  still  21$  larger  than  the 
result  obtained  with  the  NYSTAB  program.   Because  of  the  large  differ- 
ence in  results,  and  the  fact  that  the  STABL  program  gave  an  unconser- 
vative  safety  factor,  the  case  was  further  investigated. 

Because  a  total  strength  analysis  is  applied  in  this  problem,  it 
could  easily  be  hand  calculated  for  a  specific  circular  failure  surface 
by  the  <f>  =  0  method.   This  method  gave  a  factor  of  safety  F  =  1.1+5  for 
circle  1  (Figure  6.12)  confirming  the  NYSTAB  result.   Finally,  STABL 
with  the  modified  Bishop  method  was  applied  to  the  same  failure  surface, 
resulting  in  a  factor  of  safety  F  =  1.55,  again  much  closer  to  the 
NYSTAB  result  than  STABL  with  the  normal  (simplified  Janbu)  procedure. 

The  fact  that  the  simplified  Janbu  procedure,  normally  used  in 
STABL,  in  this  case  gives  a  higher  factor  of  safety  than  the  modified 
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Bishop  procedure,  contradicts  Carter's  conclusion  (1971)  that  the  fac- 
tor of  safety  reaches  a  minimum  value  when  the  moment  center  (y)  goes 
towards  infinity. 

Taking  Carter's   equation   for  the   factor  of  safety  as   a   function 
of  the  y-coordinate  of  the  moment   center,   and  applying  <j>  =   0  and  c  = 
a  constant,   the   factor  of  safety  F     can  be  expressed  explicitly  as   a 
function  of  strength   intercept   c,   unit  weight  y,   geometry  of  slope  and 
failure  surface,   and  moment   center   (y).      Replacing  the  summation  over 

slices  with  integration  over  a  circular  failure  surface,  a  mathematical 

3F 

Q 

expression  is  obtained  for  F  ,  and  ■= —  can  be  found.   Derivation  of 

c         3y  3f 

these  expressions  are  included  in  Appendix  C.II.   It  is  found  that  -r— 

goes  to  0  (zero)  for  y  -*■  °°,  hence  F  reaches  a  minimum  or_  a  maximum  for 

9F  C 
c 
y  -*■  °°,  depending  on  the  sign  of  ~ —  . 

3F  y 

F     and  -r— ^     were   calculated  for  different   slopes   and  different 
c  <3y 

circular  failure  surfaces.   The  results  are  illustrated  in  Figures  6.13 

9F 
c 
through  6.19.   Defining  F  =  lim  F  ,  it  is  found  that  ~ —  is  positive, 

and  F     =  FmaX   for  circular   failure  surfaces   that   intersect  the  ground 
00         c 

surface  with  a  steep  angle  at  the  top  of  the  slope  (a_  >  60°  to  80° , 

U3F 
c 

see  Figure  6.12).   For  small  angles  of  intersection,  -r—     is  negative 

,  _    „min 
and  F  =  F^ 

OO  Q 

Since  F  ,  the  factor  of  safety  for  moment  taken  around  the  center 
o 

of  the  circular  failure  surface  (corresponding  to  y  =  y  ) ,  in  the  case 
of  <j>  =  0,  represents  a  mathematical  correct  solution  to  the  problem, 
comparison  was  made  between  the  values  of  Fot  and  F  .   Figure  6.17  shows 
F  and  F  for  homogeneous  soil  conditions  (strength  of  embankment  soil, 
c  ,  equals  strength  of  foundation  soil,  c  ).   It  is  seen  that  whereas 
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FcXh 


FIGURE     6.13 


FACTOR  OF  SAFETY  Fc  vs.  MOMENT 

CENTER  y  FOR  VARIOUS  0.,r«=2.5h, 
cE-cF 
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FIGURE  6.14  FACTOR  OF  SAFETY  Fc  vs.  MOMENT 
CENTER  y  FOR  VARIOUS  0O  ,  r=2.5h, 
cE=2cF 
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FIGURE    6.15        FACTOR  OF  SAFETY   Fc  vs.  MOMENT 

CENTER  y  FOR  VARIOUS  0.,r-5h, ce»2c. 
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FIGURE  6.16  FACTOR  OF  SAFETY  Fc  VS.  MOMENT 
CENTER  y  FOR  VARIOUS  0. ,  r-2.5h, 
Cj-O.Scp 
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FIGURE    6.I7        STABILITY    NUMBER 
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FOR  VARIOUS   RADII,  CIRCULAR 
FAILURE  SURFACES.  Cj-c,,  *-0 
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Fc     ~     moment  around  center  of  circle 
(modified  Bishop) 

F^    ~     horizontal  force  equilibrium 
(simplified  Janbu) 
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FIGURE     6.18         STABILITY  NUMBER    Fcyh       VS.0OFOR 

VARIOUS  RADII,  CIRCULAR  FAILURE 
SURFACES.  ce«2c„*«0 
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F.   rsj    moment  around  center  of  circle 
(modified  Bishop) 

Fee  ~     horizontal  force  equilibrium 
(simplified  Janbu) 
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FIGURE    6.19        STABILITY  NUMBER    Fc^h    VS .  0.    FOR 

VARIOUS   RADII,  CIRCULAR  FAILURE 
SURFACES.cE=0.5cF,*-0 
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F  =  F    is  only  slightly  less  than  the  corresponding  F  ,  F  =  F 
co     c  o    oo     c 

may  be  considerably  larger  than  the  corresponding  F  ,  and  actually  goes 

towards  infinity  for  a  approaching  90°  (0  =  0).   It  is  seen  that 

F     =  Fmax  for  8„  <  10°   corresponding  to  or,  >  80°    (a.  =  90°   -  9  J . 
co  c  0  0  0  0 

If  the  soil  is  inhomogeneous ,  with  the  slope  material  being  stron- 
ger than  the  foundation  material,  e.g.,  an  embankment  on  a  soft  founda- 
tion, as  in  this  problem,  the  result  by  the  simplified  Janbu  method 
(F  )  deviates  even  more  from  the  mathematical  correct  solution  (F  ). 

co  O 

For  the  embankment  soil  being  twice  as  strong  as  the  foundation  soil 

3F 

(o  =  2c-,),  t-5-  is  positive  and  Fm  =  FmaX  for  a_  >  60°  ,  and  the  dif- 
r>      r    dy  C         U 

ference  between  F  and  F  increases  significantly  as  shown  in  Figure 

°o        O 

6.18. 

If,  on  the  other  hand,  the  foundation  soil  is  stronger  than  the 
embankment  soil,  F  and  F  are  almost  identically,  except  for  very  small 

co        o 

values  of  8.  (<  5°  for  the  foundation  soil  being  twice  as  strong  as  the 
embankment  soil,  see  Figure  6.19)-   However,  normally  the  embankment 
soil  will  be  the  strongest. 

The  conclusion  is  that  the  simplified  Janbu  method  cannot  be  relied 
upon  where  deep  failure  surfaces  intersect  the  ground  surface  at  the  top 
of  the  slope  at  angles  steeper  than  a  =  60° ,  coincident  with  a  relative 
large  total  strength*  in  the  zone  of  intersection. 

In  the  cases  of  circular  shaped  failure  surfaces,  Fq  provides  a 
reliable  solution.   For  this  reason  the  STABL  program  was  modified  to 
include  the  option  of  using  the  modified  Bishop  method  (F  )  for  circular 
failure  surfaces,  see  Chapter  IV.   Since  F  is  much  less  sensitive  to 


*  In  case  of  effective  stress  analysis,  a  relatively  large  strength 
intercept  c' . 


POO 

variations  in  a^ ,  it  is  recommended  that  the  modified  Bishop  method  is 
used  for  circular  failure  surfaces  in  general. 

A  sliding  block  search  with  STABL  gave  a  factor  of  safety  F  =  I.I42, 
almost  identical  with  the  NYSTAB  result.   This  confirms  that  the  Janbu 
procedure  normally  used  in  the  STABL  program  is  still  reliable  for 
failure  surfaces  where  a  is  smaller  than  approximately  60° .   The  slid- 
ing block  has  ctn  =  h^° . 

However,  STABL  with  the  simplified  Janbu  procedure  should  not  be 
exclusively  relied  upon  for  analyses  of  non-circular  failure  surfaces, 
intersecting  the  ground  surface  at  the  top  of  the  slope  with  a  steep 
angle,  coincident  with  a  relatively  large  total  strength  in  the  zone  of 
intersection. 

Results  of  the  different  analyses  are  summarized  in  Table  6.h. 

6.5  CONCLUSION  OF  COMPARATIVE  ANALYSIS 

The  comparisons  between  STABL  results  and  results  obtained  by 
other  methods  of  stability  analysis  show  that  STABL  generally  provides 
solutions  that  are  somewhat  conservative.   However,  the  difference  is 
usually  less  than  10/5.   This  is  acceptable  when  compared  to  the  much 
larger  variation  in  the  factor  of  safety  that  may  result  from  variable 
soil  strength  data,  see  Chapter  III. 

In  one  case,  however,  a  STABL  sliding  block  surface  gave  a  result 
3h%   lower  than  a  circular  surface  analysed  with  the  modified  Bishop 
method.   This  result  could  not  be  verified,  as  the  Spencer  procedure 
did  not  converge  for  this  failure  surface,  and  the  Bishop  method  only 
applies  to  circular  surfaces. 
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An  exception  to  the  general  conservative  tendency  are  the  cases  of 
deep  failure  surfaces  in  soils  with  strength  defined  mainly  in  terms  of 
a  strength  intercept  c  (or  c1).   In  such  cases  the  simplified  Janbu 
method  used  in  GTABL  may  give  erroneous  and  un conservative  results , 
and  should  not  be  used. 

For  the  general  conservative  solutions,  Janbu  (1956)  has  intro- 
duced a  correction  factor  for  the  factor  of  safety  obtained  by  the 
simplified  Janbu  method,  see  Figure  6.20.   These  correction  factors 
were  applied  to  the  results  presented  in  Chapter  II,  Table  2.3.   They 
showed  reasonable  good  agreement ,  although  the  corrected  factor  of 
safety  seemed  to  be  slightly  on  the  unconservative  side  compared  to 
the  Bishop  solution.   However,  for  the  non  homogeneous  soil  conditions 
analysed  in  Chapter  V  and  in  this  chapter,  the  correction  factors  did 
not  seem  to  work. 

A  more  thorough  analysis  is  required  to  determine  the  relationship 
between  type  of  failure  surface,  soil  profile  and  strength  properties, 
and  the  conservatism  of  the  simplified  Janbu  method  of  slices. 
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F-f0F. 
F«  ■  Factor  of  safety  by  Janbu's  procedure 
F  ■  Corrected  factor  of  safety 


shear  surface 


t 


FIGURE   6.20         CORRECTION  FACTORS  FOR  THE 
SIMPLIFIED  JANBU  PROCEDURE 
(AFTER  JANBU  et  al.,  1956) 
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VII   SUMMARY  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

7.1  WORK  AND  RESULTS 

The  present  study  is  a  continuation  of  a  previous  study  by  Siegel 
(1975),  dealing  with  computerized  slope  stability  analysis.   All  inves- 
tigations are  related  to  a  computer  program,  STABL,  developed  by  Siegel. 
STABL  performs  analyses  of  general  slope  stability  problems. 

7.1.1  Parametric  Study 

A  parametric  study  (sensitivity  analysis)  was  performed  on  the 
STABL  program  to  assess  the  precision  with  which  input  data  need  to  be 
defined.   The  results  of  this  study  show  that  variations  in  the  soil 
parameters  c  (c')  and  (J)  ((J)1)  have  the  largest  influence  on  the  stability 
of  a  slope,  in  homogeneous  soils,  as  well  as  for  more  variable  soil 
conditions. 

7.1.2  Modifications  of  STABL 

In  addition  to  the  three  original  surface  generators,  which  are 
the  circular  arc,  the  irregular,  and  the  sliding  block  type  surface,  a 
new  sliding  block  surface  was  introduced,  which  generates  the  active 
and  passive  portions  of  the  surface  according  to  Rankine's  theory, 
(1*5°  +  <f>'/2)  with  respect  to  the  horizontal  for  the  active,  and  (U50  - 
<J>'/2)  with  respect  to  the  horizontal  for  the  passive  portion.   The  new 
sliding  block  is  recommended  for  slopes  with  relatively  simple  soil 
profiles. 
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Also,  the  water  surface  option  is  extended,  so  that  it  is  now 
possible  to  define  different  piezometric  surfaces  for  different  soil 
layers.   This  is  convenient  in  the  analysis  of  problems  encountering 
artesian  water  conditions  or  a  perched  water  table.   It  may  also  be 
used  to  approximate  excess  pore  water  pressures  in  low  permeable  soil 
layers  due  to  embankment  construction,  and  to  analyse  reservoir  embank- 
ments for  rapid  drawdown. 

The  modified  Bishop  factor  of  safety  was  introduced  for  circular 
failure  surfaces,  in  addition  to  the  simplified  Janbu  method.   This  was 
done  because  of  erroneous  results  obtained  with  the  Janbu  method  for 
deep  circular  failure  surfaces  in  soils  with  <|>  =  0.   The  Bishop  method 
satisfies  overall  moment  equilibrium,  which  in  the  case  of  a  soil  with 
<f>  =  0 ,  provides  a  mathematically  correct  solution,  whereas  the  Janbu 
method  with  overall  horizontal  force  equilibrium  represents  an  approxi- 
mate solution.   The  Bishop  method  can  be  applied  only  to  circular 
shaped  failure  surfaces. 

7.1.3  Comparative  Study 

Results  of  analyses  performed  with  STABL  (simplified  Janbu  method) 
were  compared  with  results  obtained  by  other  methods  of  slope  stability 
analysis.   Results  by  STABL  are  usually  conservative,  an  exception  being 
the  deep  circular  failure  surfaces  mentioned  above.   The  differences 
are  normally  less  than  10$,  but  a  value  as  high  as  3h%   was  obtained. 
The  conservatism  seems  to  be  larger  for  sliding  block  type  failure 
surfaces,  than  for  failure  surfaces  of  circular  shape. 
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7.?     SUGGESTIONS  FOR  FUTURE  WORK 

GTARIi  has  been  used  routinely  by  the  Indiana  State  Highway  Commis- 
sion for  the  last  two  years,  during  which  period  (1976-77)  more  options 
have  been  added  to  the  program,  and  comparisons  have  been  made  with 
other  methods  of  slope  stability  analysis,  as  described  above. 

The  result  is  a  program,  even  more  versatile  than  the  original 
STABL  program.   However,  the  program  is  also  rather  bulky  and  for  many 
problems,  only  a  fraction  of  the  program  is  actually  used,  even  though 
the  total  program  has  to  be  loaded  into  the  computer  memory.   This  is 
inefficient  use  of  computer  space,  and  could  be  avoided  by  means  of  an 
overlay  structure.   Such  a  structure  divides  a  program  into  overlays, 
and  only  those  overlays  that  are  needed  for  a  particular  problem  are 
loaded  into  the  central  core  memory  of  the  comt>uter.   This  improvement 
is  recommended  for  the  STABL  program. 

A  more  thorough  analysis  is  required  to  determine  the  relationship 
between  type  of  failure  surface,  soil  profile  and  strength  properties, 
and  the  conservatism  of  the  simplified  Janbu  method  of  slices.   Janbu 
did  provide  correction  factors  to  adjust  the  simplified  solutions  to 
more  accurate  results.   However,  these  correction  factors  seem  to  apply 
only  to  slopes  in  homogeneous  soil. 

The  task  still  remains  to  develop  an  analytical  model  for  slope 
stability  analysis  that  combines  the  simplicity  of  the  Janbu  method 
with  the  higher  accuracy  of,  for  example,  the  Morgenstern  and  Price 
procedure,  and  does  not  have  the  problems  related  to  misleading  factors 
of  safety. 
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